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SECTION 1 
INTRODUCTION 



Tymshare FORTRAN IV is composed of: 1) FORTRAN IV source language statements with which the 
actual programs are written, and which are upward compatible with most FORTRAN IV compilers; and 2) 
The Conversational Compiler System CCCS) command mode which controls the operating system under which 
FORTRAN IV language statements are used on the Tymshare system. 

The FORTRAN IV Language statements are compiled and analyzed for errors one at a time as they 
are being entered. If the statement is incorrect, a diagnostic is given immediately. The statement 
then may be corrected. After modification, the program can be listed on the terminal, saved on a disk 
file, and/ or executed. The user can specify a single statement or a range of statements to be executed 
as well as execute the entire program. Program execution can be interrupted at any time, and direct 
statements can be entered for immediate execution. Program execution can be resumed at the point of 
interruption. 

Source programs in FORTRAN IV can be entered directly from the terminal, from a paper tape, or 
from a file. The user with a source program punched on cards can request that the cards be put on a 
disk file at the computer center, and then load the program from the file. 

Data can be read in from the terminal, or from a file. Data on paper tape, magnetic tape, and 
cards may be read onto a disk file and then used. Similarly, results can be printed on the terminal 
or saved on a file. Large amounts of output can be printed on the high speed printer at the user's 
request at the computer center. 

In addition to being conversational, Tymshare FORTRAN IV is compatible with other FORTRAN IV 
systems. Other advantages of Tymshare FORTRAN IV include: 

Editing and debugging features provided by CCS which make it easy to write and 
debug a program on-line. 

Large programs handled through program segmentation. 
Recursive calls to subprograms (functions, subroutines, segments) 
Dynamic storage allocation (subroutines, segments) 
Dynamic parameters in subroutines 
FORTRAN IV should be used when any of the following characteristics are desired: 
When the user wishes to do complex arithmetic 
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When the user wishes to have double precision 

When the user wants to use an existing FORTRAN IV program from other computers; or 
to debug a FORTRAN IV program to be run on other computers by taking advantage of 
the fact that FORTRAN IV is a common language. 

When the user wishes to use local names to isolate subroutines. 

When the user wishes to use long names for greater readability. 

Purpose Of The Manual 

This manual can be used both as a Tymshare tutorial guide and as a reference manual for 
FORTRAN IV. One of its aims is to provide sufficient instruction for a beginner in computer program- 
ming to write complete programs in FORTRAN IV and run them successfully on the Tymshare system. 

This manual is organized so that an experienced programmer can skip Section 2 - An Introduc- 
tion To FORTRAN! IV Language Programming - and proceed directly to Section 3 which describes the com- 
plete FORTRAN IV language. Section 4 describes the commands in CCS; that is, how programs are entered, 
stored, and executed in the Tymshare system. It includes also the CCS commands for control of running 
programs and for debugging. Section 5 contains some user aids. 

Section 6 provides a handy reference of Tymshare FORTRAN IV commands and statements. 
Section 7 gives some sample programs written in FORTRAN IV and executed on the system. 
Appendix A contains instructions for planning for efficient use of computer storage, using 
memory allocation estimation; and Appendix B contains the internal representation of ASCII codes. 

Procedures For Entering And Leaving The System And 
A Sample FORTRAN Program Entered And Executed 
On-Line 

The following is a smtmary of the log in and log out pvoaeduves for the Tymshare system. The 
complete description of the keyboard and terminal may he found in "The Tymshare Terminal" folder. A 
complete description of the EXECUTIVE user commands is found in the Tymshare EXECUTIVE Manual, Refer- 
enoe Series. 

To gain access to the Tymshare time sharing system, you must first log in. As soon as the con- 
nection to the Tymshare computer is made, the system will type: 
PLEASE LOG IN:p 

Type a Carriage Return. The system replies with: 
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ACCOUNT: A3 5 

Type your account niamber (AS in this case) followed by a Carriage Return. The system then 

types : 
PASSWORD: ^ 

Type your password followed by a Carriage Return. The letters in the password do not print. 
The system next types: 
USER NAME: JONES f 

The user name JONES is followed by a Carriage Return. The system next asks for a project code. 
PROJ C0DE:K-123-Kp 

K-123-K is a project code. NOTE: A project code is optional. If no project code is wanted, 
simply type a Carriage Return in response to the system's request. 

After you have entered the requested information correctly, the system will type: 
READY 4/8 11:20 

You are now in the EXECUTIVE system and can call FORTRAN by typing FORTRAN followed by a Car- 
riage Return. FORTRAN will reply with a > when it is ready to accept a command. 

LOGGING OUT 

To exit from the Tymshare system, you first must be in the EXECUTIVE, characterized by a dash 
in the left hand margin. To return to the EXECUTIVE system from FORTRAN type: 
QUITp 

The EXECUTIVE dash will appear in the left-hand margin. Now type: 
- LOGOUT p 

followed by a Carriage Return. The system then will type: 

TIME USED 0:37:12 
PLEASE LOG IN: 

You now may disconnect the line, or let another user log in. 

The following is a sample FORTRAN IV program entered from the terminal and executed. 
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PLEASE LOG IN: 
ACCOUNT: A3 
PASSWORD: 
USER NAME: MM 
PROJ CODE: KL-456 
READY 10/23 12:54 
-FORTRAN 

>10 *THIS PROGRAM COMPUTES BALANCE ON A SERIES OF LOANS 
>20 10 ACCEPT "LOAN NUMBER = "*NUM 

>30 ACCEPT ••PRINCIPAL= PRIN, "RATE^ R* "PAYMENT^ S'SPAY 
>40 DISPLAY "MONTH BALANCE" 
>50 MONTHS 1 

>70 20 XINTEREST=PRIN*R/12 
>80 PRINapRlN+XINTEREST-PAY 
>85. IF (PRIN.LE.O) GO TO 10 
>90 WRITE Cl*7> MONTH, PRIN 
>95. 7 FORMAT ( I 5* 1 OX, "$"* F8 . 2) 
>I00 MONTH=MONTH+ 1 
>110 GO TO 20 
>120 END 
>RUN 

LOAN NUMBER ~ 13556 
PRINCIPALS $1500.00 
RATE= .08 
PAYMENTS $125.00 



MONTH 




BALANCE 


1 


$ 


1385.00 


2 


$ 


1269.23 


3 


$ 


1152.69 


4 


$ 


1035.37 


5 


$ 


917.28 


6 


$ 


798.39 


7 


$ 


678.71 


8 


$ 


558.24 


9 


$ 


436.96 


10 


S 


314.87 


11 


$ 


191.97 


12 


$ 


68.25 



LOAN NUMBER = 
INTERRUPTED BEFORE: 

20. 10 ACCEPT "LOAN NUMBER - "*NUM 

20. >QUIT 

-LOG 

TIME USED 0: 11:56 
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SECTION 2 

INTRODUCTION TO CCS FORTRAN PROGRAMMING 

A. INTRODUCTION TO FORTRAN IV LANGUAGE ELEMENTS 

A program is made up of the statements necessary to tell the computer how to solve a problem - 
how to accept data, how to process it, and how to return the results to the user. 

To analyze a problem to be solved, we must first examine the data and decide what calculations 
are necessary to produce the desired results using the data available. 

Let us compute the interest on a loan. The basic computation is to multiply the principal by 
the interest rate. It may be written as: 
PRINCIPAL * RATE 

This is a FORTRAN expression and its value could be assigned to a variable written in an assignment 
statement . 

XINTEREST = PRINCIPAL * RATE (* is the FORTRAN IV multiplication symbol.) 

Assume that the interest is compounded monthly. The borrower makes a monthly payment to cover 
the interest charge and to reduce the principal. Only two assignment statements are necessary to make 
the necessary computation: 

XINTEREST = PRINCIPAL * RATE/ 12 (/ is the FORTRAN IV division symbol) 

BALANCE = PRINCIPAL + XINTEREST - PAYMENT 

The first assignment statement computes the monthly interest charge; the second statement com- 
pounds the interest and deducts the payment to determine the balance due. 

Before the computations can be made, however, the computer, must accept the data; that is, it 
must have the values for PRINCIPAL, RATE, and PAYMENT. An ACCEPT statement is required to obtain the 
data. 

ACCEPT PRINCIPAL, RATE, PAYMENT 
XINTEREST = PRINCIPAL * RATE/ 12 
BALANCE = PRINCIPAL + XINTEREST - PAYMENT 

The ACCEPT statement is an input statement for reading data into a program from the terminal . 

After the data is read, the assignment statements are executed. The results of the processing 
are then ready to be returned. 

A DISPLAY statement is added to tell the computer to print the results. 

ACCEPT PRINCIPAL, RATE, PAYMENT 
XINTEREST = PRINCIPAL * RATE/ 12 
BALANCE = PRINCIPAL + XINTEREST - PAYMENT 
DISPLAY PRINCIPAL, RATE, PAYMENT, BALANCE 
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The DISPLAY statement is an output statement which directs the computer to write on the termi- 
nal the current values o£ the following list of variables: 
PRINCIPAL, XINTEREST, PAYMENT and BALANCE. 

At this point the statements provide all the instructions necessary to read, compute, and 
write the information for a single loan record. 

Suppose that many loans must be processed. The program must loop; control must be transferred 
from the last of the sequence of statements back to the first. A transfer of this type can be made 
only to a statement having a statement number , which is a number chosen by the programmer and prefixed 
to the statement to which the transfer is to be made. Any number may be used. 

Now a GO TO statement is used to transfer control. 

100 ACCEPT LOAN-NUMBER, PRINCIPAL, RATE, PAYMENT 
XINTEREST = PRINCIPAL * RATE/12 
BALANCE = PRINCIPAL + XINTEREST - PAYMENT 
DISPLAY LOAN-NUMBER, PRINCIPAL, RATE, PAYMENT, BALANCE 
GO TO 100 

The GO TO statement is one of the control statements in FORTRAN IV. 

We have added LOAN-NUMBER to the list of variables to be read and written. Now with this pro- 
gram many loans with varying rates of interest can be processed. 

To amplify the problem above, suppose that instead of processing many loans, the balance on a 
single loan is to be processed for all of the months for which payment is made and the results printed 
in the form of a chart. 

Again, the program must loop, this time for each month's payment. The number of the month is 
given an initial value of one and then increased by one each time that the loop is executed. 

In this program an additional check must be made to determine when the entire principal has 
been paid, that is, when the balance is less than or equal to zero (BALANCE .LE. 0). This check is 
made using a logical IF statement. If the logical expression (BALANCE .LE. 0) is true, the statement 
following it (GO TO 90) is executed and control is transferred out of the loop to the statement 
labeled 90. If the logical expression is false, that is, the balance is not less than or equal to 
zero, the statement GO TO 90 is ignored and the next statement in sequence (DISPLAY MONTH, BALANCE) 
is executed. 

The fourth statement in the program (BALANCE = PRINCIPAL) initializes the variable BALANCE 
to the value of PRINCIPAL that was accepted as input. Each time through the subsequent loop, the 
new value of BALANCE is used to calculate XINTEREST. 

The text enclosed in the quote marks will be printed on the terminal when the program is 
executed. 
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ACCEPT "PRINCIPAL = PRINCIPAL, "RATE = RATE, "PAYMENT = ", PAYMENT 
DISPLAY "MONTH BALANCE" 
MONTH =1 

BALANCE = PRINCIPAL 
20 XINTEREST = BALANCE * RATE/ 12 

BALANCE = BALANCE + XINTEREST - PAYMENT 

IF (BALANCE .LE. 0) GO TO 90 

DISPLAY MONTH, BALANCE 

MONTH = MONTH + 1 

GO TO 20 
90 STOP 

END 



The program will now produce a chart which shews the monthly balance on the loan. This pro- 
gram is essentially the one first presented on page 4. 
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B. KEY STEPS IN PROGRAMMING 

As shown above, a program can be written in a computer programming language such as FORTRAN IV 
to solve a problem. In this section we will outline the key steps in writing a program. 

To write a program, the key steps are: 

1. Defining the problem 

2. Selecting a method for solution 

3. Analyzing the problem 

4. Writing the instructions in the language 

5. Debugging and checking the program 

6. Documenting the program 

A brief description of each step is as follows. 

lij. Defining The Problem 

Before you can write a program to solve a problem, you must know what the problem is. To de- 
fine a problem you should consider: 

a. What information is given: Cl^put) 

b. What answers are required:, (Output) 

c. In what manner is the data supplied? 

d. How are answers to be given and what accuracy is desired? 

e. Are there any options? 

2_. Selecting The Method 

Usually there are several ways to solve a problem. Decide on one that is best for your par- 
ticular problem. You may begin by preparing a brief step-by-step numerical solution to your problem 
using as many methods as you can. After some experience, it will be easier to see which method is 
better for computer solution. 

3_: Analyzing The Problem 

Once the method for solution is decided, you may decide on the steps for solving the problem. 
A program may be solved primarily by evaluating a sequence of formulas, or may require more involved 
steps. In the latter case, a detailed list of steps in mathematical form or by means of a flowchart 
is necessary. A flowchart is a pictorial representation of the problem and is used to illustrate the 
logical flow of steps necessary to solve the problem. 

Conventionally: 

Rectangular boxes are used for input, output 
and computations 
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Diamond- shaped boxes are used for decisions. 

Circles are used as connectors when flowcharts 
become more complicated. 
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Figure 1 is an example of a simple flowchart for a program to find the sum of the reciprocals 
of n numbers. 

Usually it is worthwhile to prepare test cases that exercise all options; that is, follow all 
paths of the flowcharts... A test case should contain input data for which correct answers are known. 
±1 Writing The Instructions (Coding) 

Writing the instructions usually is called coding, and is the actual writing of the program 
in a language "understood" by the computer. FORTRAN IV is a procedure-oriented rather than a machine- 
oriented language, and is made up of a small group of statement types. Each step in the problem solu- 
tion frequently corresponds to a single FORTRAN IV statement. In forming these statements, the correct 
order and proper syntax of the language must be used. You may write down all instructions before typ- 
ing them into the machine from the terminal. This allows you to do much of your initial thinking while 
not connected to the computer. 

5^ Debugging 

Debugging is testing and checking a program. When you are entering the instructions from the 
keyboard, the interactive features of the Tymshare system can be a great help to you in finding and 
correcting errors in your program. Any faulty statement will be detected and printed with an error 
diagnostic. The first step in debugging is to make the program run on the machine. The next step is 
to make sure that the answers given are correct. After FORTRAN syntax errors are corrected and the 
program runs, it is advisable to test it using data for which the correct answer is known. The final 
checkout is of course continued satisfactory use of the program. 
6. Documenting The Program 

It is a good idea to include comments in the program to remind yourself of what the program 
does and what you have done in the different steps. After the program is debugged, document the pro- 
gram if you want to use it again. 

Now that we have introduced you to programming, we will proceed to a simple course in FORTRAN 
IV, which will enable you to write and use programs of your own on the Tymshare system. We suggest 
that you try to run some simple programs on the terminal to see how it works. FORTRAN IV used on the 
Tymshare conversational system is very easy to use. You will find that the actual use of the computer 
is your most valuable learning experience, and the fastest way to get acquainted with the system. 
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INITIALIZE: 
COUNTER = 
SUM = 



INPUT: N 
(HOW MANY NUMBERS) 



INPUT: 
NEXT NUMBER 



X 

CALCULATE: 
RECIPROCAL 
OF NUMBER 



ADD: 

RECIPROCAL TO SUM, 
INCREASE COUNTER 
BY 1 




Figure 1 - Example Of A Flowchart 
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C. A SAMPLE FORTRAN PROGRAM EXECUTED IN CCS 

We are going to take a sample problem and go through, the programming steps to arrive at a 
solution. We will then show how the program is entered and executed on the terminal. 

Problem 

Write a general FORTRAN IV program to solve two simultaneous linear equations in two unknowns. 
In particular, give answers to X and Y in the following set o£ equations: 

X + 2Y = -7 
4X + 2Y = 5 

The problem can be generalized as: 
X . Y = 

A3 X . A4 Y = 

Different sets of two simultaneous linear equations could be solved by applying different 
values of B^ and or different values of A^^, A2, A^, and A^. 

Solving for X and Y, the solution appears as follows: 



When 



(A^A^ - '^3^2) CAi\ - ^3^2) 

the denominator A^A^ - A^A^ is equal to zero, there is no unique solution. 



Coding 



Here is the coding of the program: 



10 ACCEPT Al, A2, A3, A4 
15 DIFF = Al * A4 - A3 * A2 
20 IF CDIFF .EQ. 0) GO TO 100 
25 ACCEPT Bl, B2 
30 X = (Bl * A4 - B2 * A2)/DIFF 
35 Y = (B2 * Al - Bl * A3) /DIFF 
40 DISPLAY X, Y 
45 GO TO 200 

50 100 DISPLAY " NO UNIQUE SOLUTION" 
55 200 STOP 
60 END 

Each statement corresponds to a step in the flowchart. Every statement has a line number (10, 
15, 20, ...). The statements are executed in the order in which they are numbered unless a transfer 
occurs. Each line number identifies a line. Line numbers keep the program statements in order. As 
we shall see later, line numbers are used to refer to the program statements when modifying and insert- 
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^ START ^ 







INPUT: 
Ai, A2, A3, A4 






CALCULATE 
THE DIFFERENCE: 
A1A4-A3A2 








YES 



INPUT: 
Bi, B2 



CALCULATE: 
X, Y 



OUTPUT: 
X, Y 



OUTPUT: 
A MESSAGE 
'NO UNIOUE SOLUTION 
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ing statements. The choice of line nximbers is arbitrary; they may be any number from .001 to 999.999, 
assigned either explicitly by the programmer or implicitly by the system. 

There are two lines in the program with statement niombers in addition to line numbers (lines 
50 and 55) . These are statement numbers . Statement numbers are part of the FORTRAN IV language and 
are used as reference points by the program. Statement numbers can be any integer number and need not 
be in order. 

Now we will go through the program step by step. 
10 ACCEPT Al, A2, A3, A4 

The first statement (10) is an ACCEPT statement for reading input data. The ACCEPT statement 
is unique to Tymshare FORTRAN IV and allows data to be read in a format-free form from the terminal. 
By format-free form we mean that the user simply types in the data values separated by a comma, a 
space, or a Carriage Return. The data is not input according to a specific format. When the program 
is executed, the system will ring a bell and wait for the user to type in the data values specified 
in the ACCEPT statement. Al, A2, A3, and A4 are variable names in which the values are to be stored. 

15 DIFF = Al * A4 - A3 * A2 

The next statement (15) is an assignment statement . In this statement the value of the expres- 
sion Aj^A^ - k^k^ is assigned to the variable DIFF. DIFF is another variable name. A variable name 
can be one or more alphabetic characters or digits, but the first character must be a letter of the 
alphabet. It is good practice in programming to choose a variable name that has some meaningful rela- 
tionship to the program so it can be read easily. Notice two things: 1) multiplication is stated 
explicitly with a '*' sign. This is one of the arithmetic operators in FORTRAN IV. The others are 
'+' for addition, '-' for subtraction, •/' for division, and for exponentiation (raising a number 
to a power). 2) The correct form of an assignment statement is one in which the left side is the name 
of a single variable. The equal sign means "is to be replaced by" rather than "is equivalent to". 
This is an important distinction. For example, X = X + 1 is not an equation but a valid assignment 
statement. 

20 IF (DIFF .EQ. 0) GO TO 100 

In this statement we ask the question: "Is the difference equal to zero?" This type of state- 
ment is one of the control statements for making decisions. This statement allows you to branch to a 
different part of the program if the answer to the above question is "yes". This type of statement 
is an example of a logical IF statement . If the condition (DIFF equal to zero) is true, the statement 
following (GO TO 100) will be executed; otherwise GO TO 100 is not executed and the program goes on 
to the next statement. The statement GO TO 100 here is called the conditional GO TO statement. It 
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consists of the word GO TO and a statement number. It transfers operation to the part of the program 

specified by the statement number Cin this case, statement 100), if and only if the condition speci- 
fied is true. 

.EQ. stands for "equal to" in FORTRAN IV, and must be surrounded by dots. .EQ. is one of the 

logical operators in the language. The others are not equal to C-NE.), less than (.LT.), less than 

or equal to (.LE.), greater than or equal to C«GE.), and greater than C-GT.). We call the expression 
DIFF .EQ. a logical expression . 

25 ACCEPT Bl, B2 

Line 25 is another ACCEPT statement to read in the values of Bl and B2. 

30 X = (Bl * A4 - B2 * A2)/DIFF 

35 Y = CB2 * Al - Bl * A3)/DIFF 

Statements numbered 30 and 35 complete the computation of the solution. Notice that the de- 
nominator is evaluated previously as the variable DIFF and is used in these two statements. Notice 

also that parentheses are used to specify that the numerator is the entire quantity (Bl * A4 - B2 * A2) . 

B A B 

FORTRAN distinguishes between the quantities A - ^ and — ^ — . The first is written as A - B/C, the 
second as (A - B)/C. Parentheses may be used to specify the order of operations in expressions. If 
parentheses are omitted, the order of arithmetic operations is -b, * and /, + and -. 

40 DISPLAY X, Y 

The next statement (40), a DISPLAY statement , is to output answers. When the statement is en- 
countered, the values of variables X and Y are printed on the terminal. The DISPLAY statement is for 
format-free output just as the ACCEPT statement is for format-free input. This means that the system 
supplies the format that is used for output. Format-free input and output statements make it very 
easy to input and output data. 

45 GO TO 200 

Line 45 is another control statement. It is an unconditional GO TO statement which transfers 
the program to statement 200. We put in this transfer because, after displaying the solutions to 
X,Y, we simply want to end the program. This transfer allows us to bypass execution of the next state- 
ment and go on to the end of the program. 

50 100 DISPLAY "NO UNIQUE SOLUTION" 

The next statement has a statement number (in addition to a line number). Earlier in line 20 
we have told the program to transfer to this statement if the difference (Al * A4 - A3 * A2) is zero. 
Any number can be used as a statement number. This statement prints a message NO UNIQUE SOLUTION. 
Any message (text) inside a pair of double quote marks can be printed with a DISPLAY statement. 
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55 200 STOP 



The STOP statement stops execution o£ the program. Notice that it also has a statement number. 
Earlier in the program we have asked for a transfer to here after displaying the results. 



60 END 



The last statement is the END statement . Every Tymshare FORTRAN IV program must end with an 
END statement. The END statement may not be labelled with a statement number; that is, we cannot 
transfer to this statement from other parts of the program. 

Before we go on to show you how to run this program on the computer, we will summarize what 
we have learned about FORTRAN IV from this example. 

1. ACCEPT statements are used to input data values used for computation. 

2. DISPLAY statements are used to output the results of computation, as well as any literal 
text. 

3. Variables are used to identify the values of input data and to store computational results. 

4. The computer evaluates expressions whicTv are made up of variables and operators. Paren- 
theses may be needed to specify the order of computation. 

5. An assignment statement stores the value of an expression in a single variable. 

6. Control statements allow the user to override the normal sequential order of execution in 
the program. Here we have seen two control statements; namely, the logical IF with a con- 
ditional GO TO, and the unconditional GO TO statement. Logical operators and logical ex- 
pressions are used with logical IF statements. 

7. Statement numbers such as 100, 200 are used to label statements in the program. 

8. A STOP statement stops execution of the program. 

9. The last statement of a FORTRAN IV program must be an END statement. 

To execute the program on the computer, log in and call FORTRAN. When FORTRAN is ready to 
receive instructions, it prints the > symbol. Then you may start typing in the statements of your 
program with line numbers, ending each line with a Carriage Return. If you type an incorrect state- 
ment, the computer will ring a bell and print out the statement in error with an f pointing to the 
mistake. Just retype the statement correctly. When you have finished typing all the statements, the 
computer will type a > and the program is ready to be executed. 

To start execution of a program, type RUN or EXECUTE followed by a Carriage Return. 
The following is the program typed on line and executed. 



12 



PLEASE LOG liit ASIMMI 

READY 9/2 5 9:36 
-FORTRAN 

>10 ACCEPT A1,A2*A3#A4 

>15 DIFF=Al*A4-A3*Aa 

>20 IF (DIFF .EQ.O> GO TO 100 

>2 5 ACCEPT B1#B2 

>30 X=(B1*A4-B2*A2)/DIFF 

>35 Ys(B2*Al -B1*A3>/DIFF 

>40 DISPLAY X*Y 

>45 GO TO 200 

>50 100 DISPLAY "NO UNIQUE SOLUTION" 
>55 200 STOP 
>60 END 

>RUN 

2> 5>9> 4 

7*8 

0.3243243243 1.27027027 
STOPPED AT: 55 

When execution o£ the program is completed, the message STOPPED AT 55 is printed, 55 being 
the last statement executed before END was encountered. 
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To save a debugged program for future use, type SAVE/File Name/^ . The file name can be any 
name you want to give the program. After we have executed the above program, we could do the followr 
ing : 

>SAVE. /EQUATION/;) 
NEW FILE 

The computer tells you whether it is a NEW FILE or an OLD FILE. Saving the program on an old 
file will erase the contents of that file and replace them with the current program. 

Type a Carriage Return after NEW FILE. Type a Carriage Return after OLD FILE if you do not 
want to keep the file. If you want to keep the old file, type ALT MODE or ESCAPE and choose another 
file name. The entire program is now saved with line numbers on the file named /EQUATION/, and a > is 
given . 

To use this program in the future type: 
>LOAD /EQUATION/ 

The computer will say O.K. and start loading. When it has finished, it will type a >. Now 
you can execute your program again. 

Suppose you want to have a listing of the program in the file /EQUATION/. After you have 
loaded the file, type LIST ^ as shown in the following example. 

>LOAD /EQUATION/ 
OK. 

>LIST 

10* ACCEPT A1»A8«A3«A4 

IS* DIFP«A1*A4-A3*A2 

80. IF CDirr .EQ.O) (30 TO 100 

as. ACCEPT B1«B2 

30. XMBt*A4«B8*A8>/Dirr 

35. Y«€Ba4iAI-Bl4iA3>/DirF 

40. DISPLAY X,Y 

45. aO TO 800 

50. 100 DISPLAY ••MO UNIQUE SOLUTION" 

55. 200 STOP 
60. END 



The program is printed neatly on the terminal with line numbers. The statement numbers are 
aligned and the statements indented. 
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The commands RUN, SAVE, LOAD, LIST are part o£ Tymshare FORTRAN IV commands. They are not 
parts of the FORTRAN language but are operating features of the Tymshare Conversational Compiler Sys- 
tem (CCS). You will learn about the other CCS commands later in the manual. 
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D. A SECOND EXAMPLE TO SHOW THE USE OF OTHER FORTRAN FEATURES 

AND CONVERSATIONAL COMMANDS 

In the following example we will introduce you to some other FORTRAN language elements as 
well as other Tymshare FORTRAN CCS commands. 

Problem 

Write a program to calculate the mean and standard deviation of N numbers, where N is greater 
than 1. Print also the normalized data. 

Input : Required are: 

N - the number of observations 
Ai through An - the observations 

Compute mean and standard deviation: 



Mean = 



n 

I 

i=l Ai 
N 



Standard Deviation = 



n 

E Ai 



2 i=l 



N - 1 



Normalized data = „^ ~ f?^^" ^. i = l,2,...,n 

Standard Deviation 

Output : 

Original data, normalized data, mean, and standard deviation. 
A(l) - ACN) NA(1) - NACN) MEAN STD 

Program Coding 
10 DIMENSION A(15) 
20 REAL N, MEAN, NA (15) 
30 ACCEPT "NUMBER OF OBSERVATIONS", N 
40 ACCEPT (A(I),I=1,N) 
50 SUM=0. 

60 SUMSQ=0. 

70 DO 20 1=1, N 

80 SUM=SUM + ACI) 

90 SUMSQ=SUMSQ + A(I)f2 

100 20 CONTINUE 
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^ START ^ 







INPUT: 
N, 

Aj, i = 1 TO N 


> 




CALCULATE: 
SUM OF N NUMBERS 
SUM OF SQUARE 
OF N NUMBERS 






CALCULATE: 
MEAN, 
STANDARD DEVIATION 






CALCULATE: 
NORMALIZED DATA 






OUTPUT: 
ORIGINAL DATA. 
NORMALIZED DATA 
MEAN, 
STANDARD DEVIATION 







STOP 
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110 MEAN = SUM/N 

120 STD = SQRT[(SUMSQ-SUM+2/N)/CN-l)] 

130 DO 50 1=1, N 

140 NA(I) = (A(I)-MEAN/STD 

150 50 CONTINUE 

160 DISPLAY "INPUT DATA", A 

170 DISPLAY "NORMALIZED DATA",NA 

180 DISPLAY "MEAN=", MEAN, "STAND. DEV. =", STD 

190 END 

In this program we are reading the N data values into a subscripted variable or an array, A. 
Thus A(l) is the first data observation and A(N) the Nth data observation. Storing data into an array- 
makes it possible to refer to the list of observations by a single name. 

To ask the computer to reserve space for the array A we use a DIMENSION statement as is done 
in line 10. Fifteen spaces are reserved for array elements AQ) to AC15) in the array A (for input 
data) . 

Line 20 is a type declaration statement , specifying that the variables N and MEAN are real 
numbers and the array NA contains real data Cnumbers with decimal points) . This is necessary because 
FORTRAN distinguishes between integer numbers and real numbers. Thus 3/2=1, but 3./2.=1.5 in FORTRAN. 
Any variable name starting with I through N inclusive is treated as an integer variable unless speci- 
fied to be a real number as we have done here. If this statement is left out, integer divisions will 
be performed and integer results will be given. The array NA is dimensioned using implied dimension- 
ing in the type declaration statement. 

Line 30 is an ACCEPT statement with descriptive text in quote marks. The text will be printed 
on the terminal when the program is executed. Line 40 reads in the N data values and stores them in 
the array A. 

Lines 50 to 100 correspond to box 2 of the flowchart. These statements sum the N numbers and 
the square of the N numbers. Lines 50 and 60 initialize the variables SUM and SUMSQ to zero] The 
statement, DO 20 I=1,N controls the repetition of all succeeding statements down through and including 
the statement labeled 20. Repetition of these statements is controlled by varying an index called I, 
from an initial value of 1 to a terminal value of N in increments of 1. Lines 70 to 100 create a 
DO_loo£ which sums the N numbers and the squares of the N numbers; each time adding the present value 
or square to the variables SUM and SUMSQ, the previously cumulated values. The CONTINUE statement on 
line 100 indicates that the process is to be continued. 
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After we calculated the sum (^I^ A^) and the sum of the squares (.S^ ^i^} > we can proceed to 
calculate the mean and standard deviation. Line 110 calculates the mean. In calculating the standard 
deviation on line 120, we use the library function SQRT to find the square root. The expression 
(argument) is surrounded by square brackets. SQRT is one of the many built-in functions provided by 
the FORTRAN library. 

Lines 130 through 150 correspond to box 4 in the flowchart. Another DO loop is used to normal- 
ize each of the N data values. 

The last three statements of the program are output statements. Literals within quote marks 
are displayed describing what the output is in each case. 

Before we run this program on the computer we will summarize the FORTRAN IV features in this 

example. 

1. Arrays are dimensioned in DIMENSION statements. 

2. Specific types of variables can be declared with a type declaration statement. Arrays 
may be dimensioned in a type declaration statement. 

3. ACCEPT and DISPLAY statements can contain descriptive text as well as variables. 

4. All variables must be defined when used on the right side of a replacement statement; for 
example, SUM (line 50) must be initialized (to zero in this case) before it can be used in 
line 80. 

5. A DO loop is used for a repetitive process. 

6. The FORTRAN library of functions can be used to simplify computation. The argument must 
be enclosed in square brackets. 

Tymshare FORTRAN IV provides commands which make it easy to create programs on the terminal. 
These commands include line prompting, debugging, and editing features. 

If you do not want to type in line numbers when you are first entering a program, ask the com- 
puter to prompt you with line numbers. 

To do this, type in a line number followed by an increment in parentheses. (An implied ENTER 
command; the word ENTER is optional.) The computer will prompt you with line numbers starting with 
the number you specify increasing each time at the increment specified. A Carriage Return is neces- 
sary at the end of every statement. When the last statement of the program has been entered, you must 
type a Control D (D^) to terminate the ENTER command. 

Many editing features are available to you when entering a statement. For example, use a Con- 
trol A (AC) to delete the last character typed. 
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If you have entered a faulty statement, the computer will ring a bell when you hit the Carriage 
Return and print out the statement with an + pointing to the first character that caused the error. 
This statement becomes the old line for editing purposes and any of the control characters described 
in this manual may be used. 

In creating the following program, we will show you how to use some of these features. 



>10C5> 

10. DIMEMSION ACS) 

15. ACCEPT ••NUMBER OF OBSERVATIONS ' 

20. ACCEPT (A(I)#I=1#N) 

25. SUM=0. 

30. S(JMM-SQ=0. 

35. DO 20 I=l#N 

40. SUM=SUM+A<I) 

45. SUMSQ=SUMSQ+A(I)t2 

50. 20 CONTINUE 

55. MEAN=SUM/N 

60 . STD=SQRTC ( SUMSQ-SUM*SUM/N/(N-n C 



N 



This command caused the computer to 
prompt with line numbers beginning 
with line 10 in increments of 5. 



The user typed another M instead of 
S. He typed a Control ACA^) to erase 
the M, then typed the rest of the 
line. 



60. STD=SQRTC(SUMSQ-SUM*SUM/N/CN-1 ) C 

t 

60. STD=SQRTC(SUMSQ-SUM*SUM/N)/CN-1 )3 

65. DO 50 I=1>N 

70. NACI)~(A<I)-MEAN)/STD 

75. 50 CONTINUE 

80. DISPLAY ••INPUT DATA'»>A 

85. DISPLAY ••NORMALIZED DATA"*NA 

90. DISPLAY ''MEANS "^MEAN* ''STAND.DEV. 

95. END 

100. 



>12 REAL N*MEAN 



STD 



This statement is in error. The com- 
puter rang a bell and typed out the 
line with an f pointing to the error. 
The user typed a Control Z(ZC) followed 
by N. Line 60 up to the first N is 
copied. He typed in the rest of the 
line. 



When the computer prompted line 100, 
the user typed a Control DCDC) to end 
the line prompt. The user had forgotten 
to declare N, MEAN to be real variables, 
so he simply typed in the statement as 
line 12. 

This statement will be inserted between 
lines 10 and 15. 



So far the program has been entered. Now it is ready for execution. When an error is de- 
tected during execution, an error message is printed together with the statement that caused the error; 
execution will be terminated. At this point you may enter direct statements for execution to find out 
what caused the error. A direct statement is preceded by an @ sign. The statement can be entered 
when a > is given. A Direct statement, once executed, is discarded and is not part of the program. 
Once the error is detected and corrected, the program may be executed again. 

In the following examples, the program just entered is to be executed. 



19 



>RUN 

NUMBER OF OBSERVATIONS 
1,2,3,4,5, 

SUBSCRIPT OUT OF RANGE 
20. ACCEPT (A(I),I=1,N) 



10 



The user tried to execute the program, 



The computer typed out an error message and the statement 
in error after 5 data values were entered, then it typed 
20. >. 



20. >@DISPLAY I 
6 



The user typed in @ for iiranediate execution of the state- 
ment: DISPLAY I 



20. >DEFINITIONS A Then he typed the CCS command DEFINITIONS A which listed 

out the statement in which A is dimensioned. He found out 
10. DIMENSION AC5) that he had not reserved enough space for the array. 

20. >10 DIMENSION A(10) The computer continued to prompt 20. >. The user typed 

line 10 again to re-dimension A. The old line 10 is re- 
placed by this line. 

^RUN The user tried to execute the program again. 

NUMBER OF OBSERVATIONS 10 

1,2,3,4,5,6,7,8,9,10 

SUBSCRIPTED VARIABLE NOT DIMENSIONED 

70. NA(I)=(A(I)-MEAN)/STD Another error is detected, the line that caused the error 

is listed. The user realized he had not dimensioned the 
>12 REAL N,MEAN,NA(10) array NA. He retyped line 12 including NA as a real array. 

^^UN He ran the program again, this time successfully. 

NUMBER OF OBSERVATIONS 10 

1,2,3,4,5,6,7,8,9,10 

INPUT DATA 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 

NORMALIZED DATA -1.486301082 -1.156011953 -0.8257228238 

-0.4954336943 -0.1651445647 0.1651445647 0.4954336943 

0.8257228238 1.156011953 1.486301082 
MEAN= 5.5STAND.DEV.= 3.027650354 

STOPPED AT: 90. 

Two things should be noticed. First, the command DEFINITIONS is used to type out the dimen- 
sion statement for A. It is one of the commands provided by the CCS system to help you debug your 
program. Secondly, when the user types in line 10 again, it replaces the old line 10. 

If you want the output to be in a neat form, you may use a formatted output. 

A formatted output for mean and standard deviation may be: 
WRITE (1,100) MEAN, STD 
100 FORMAT ("MEAN=",F10.3,"STD=", F6.3) 

The 1 in the WRITE statement is the file number and tells the computer that the values of the 
variables are to be written on the terminal (1 is reserved for terminal output). 100 is an arbitrary 
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number used to identify the format used. The values of the variables MEAN and STD are to be printed. 

The FORMAT statement has a statement number corresponding to the one specified in the WRITE 
statement. Descriptive text to be printed is to be enclosed in quotes. F10.3 specifies the format 
in which MEAN is to be printed. The F-format prints a real number. F10.3 causes a maximum of ten 
characters (including the decimal point) to be printed, with three decimal places. STD is printed 
with the next format F6.3. STD is printed with a maximtmi of six characters with three digits after 
the decimal. 

The formatted output for writing the data and the normalized data of the above may be 
WRITE (1,200) (A(I),NA(I),I=1,N) 
200 FORMAT (FIO . 3, 5X,F10. 3) 

Data to be printed is A^ and NA^ for N values of each. 

Each element of A and NA is to be printed as a real number with three places after the decimal 
point. 5X gives five spaces between the two numbers. 

The program may be changed as follows to use formatted output. The user typed in the implied 
ENTER command 80:90, which allowed him to enter the statements. The computer prompts with an @ sign 
at the beginning of each line. The statements entered will be numbered uniformly between 80 and 90 
when a terminating D^ is given. NOTE: The new lines will replace any old lines in that range. 

>80:90 

©DISPLAY " INPUT DATA NORMALIZED DATA" 
@WRITE(1,200) (A(I) ,NA(I),I=1,N) 
@200 FORMAT (FIO. 3, 5X, FIO. 3) 
@ WR ITE ( 1 , 30 ) , MEAN , STD 

@300 FORMAT ("MEAN=", F6.2,2X, "STAND. DEV .=",F6.3) 

He typed a Control D (DC) to tell the 
computer that he has finished entering 
statements. The computer prompted a >. 
>LIST 80:90 He listed out the statements in this 

range . 

80. DISPLAY " INPUT DATA NORMALIZED DATA" 

82. WRITE(1,200) (A(I) ,NA(I) , I=1,N) 

84. 200 FORMAT(F10.3,5X,F10.3) 

86. WRITE (1,300) MEAN, STD 

88. 300 FORMAT ("MEAN=", F6.2,2X, "STAND.DEV .=",F6.3) 

Now, when the program is re-run, the output is aligned in the form that the user wished it to 

be. 
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>RUN 

NUMBER OF OBSERVATIONS 10 
1,2,3,4,5,6,7,8,9,10 

INPUT DATA NORMALIZED DATA 



1.000 
2.000 
3.000 
4.000 
5.000 
6.000 
7.000 
8.000 
9.000 
10.000 



-1.486 
-1.156 
-0.825 
-0.495 
-0.165 
0.165 
0.495 
0.825 
1.156 
1.486 



MEAN= 5.50 STAND. DEV 



STOPPED AT: 88. 



3.027 



Formatted input also may be used for reading in the data. Suppose in this example that the 
data comes from a score sheet as follows: 



11 
8 
9 

10 



75 
64 
87 
59 



where the double line is the decimal place. You may write the following statements to read the input 
with a format: 

READ (0, 500) CA(I), 1=1, N) 
500 FORMAT CF5.2) 

To input data from the terminal, you merely type 
1175 864 987 1059 

allowing five columns for each number. You need not type in the decimal point. Each number would be 
read in with the format F5.2, that is, five digits for each nxamber, two of which are decimal digits. 

Now the program is modified to use formatted input as follows: 

>20 READ(0,500) (A(I),I=1,N) 



>22 500 FORMAT (4F5. 2) 

New line 20 replaces old line 20. Line 22 is inserted between lines 20 and 25 
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The program is executed: 



>RUN 

NUMBER OF OBSERVATIONS 4 
1175 864, 987 1059 
INPUT DATA NORMALIZED DATA 

11.750 1.179 
8.640 -1.206 
9.870 -0.262 
10.590 0.289 
MEAN= 10.21 STAND. DEV. = 1.303: 



STOPPED AT: 88. 



And finally, here is the listing of the entire program: 



>LIST 






10. 




DIMENSION AC 10) 


12. 




REAL N*MEAN#NA(10> 


15. 




ACCEPT "NUMBER OF OBSERVATIONS 'SN 


20. 




READ(0#500) (A(I)*Ial*N) 


22. 


500 


FORMAT <4F5. 2) 


25. 




SUM«0. 


30. 




SUMSQaO. 


35. 




DO 20 I=1>N 


40. 




SUMaSUM+ACI) 


45. 




SUMSQaSUMSQ+AC I ) t2 


50. 


20 


CONTINUE 


55. 




MEAN=SUM/N 


60. 




STD=SQRTC<SUMSQ-SUM*SUM/N)/(N-1) ] 


65. 




DO 50 I«1*N 


70. 




NA<I)aCA<I)-MEAN)/STD 


75. 


50 


CONTINUE 


80. 




DISPLAY INPUT DATA NORMALIZED DATA' 


82. 




WRITEC1*200) (A(I)#NA(I>* I«l*N) 


84. 


200 


F0RMAT(F10.3*5X#F10.3) 


86. 




WRITEC 1*300) MEAN#STD 


88. 


300 


FORMAT C "MEANS" , F6 . 2 * 2X# " STAND . DE V. =" * F 6 . 3 ) 


95. 




END 
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SECTION 3 
FORTRAN IV LANGUAGE 



A. STATEMENT ELEMENTS 

K General 

The FORTRAN IV language is composed primarily of combinations of constants, variable names, 
array names, together with the usual arithmetic, logical, and relational operators. There are certain 
rules for their formation and combination, which are presented here. 

The result of the combination of constants and variables with arithmetic operators is an 
arithmetic expression. An arithmetic expression has a numeric value. The result of the combination 
of arithmetic expressions with relational operators, or of logical variables with logical operators, 
is a logical expression. A logical expression always has the value of .TRUE, or .FALSE. 

2. Constants 

A constant is a quantity which appears in a FORTRAN statement in numerical form. 
Integer Constants 

An integer constant is a positive or negative whole number up to eleven digits such as: 


-1245 
3859437 

Real Constants 

There are two forms of real constants, either with or without exponents. 

Without exponents, the constant contains one or more digits and a decimal point (.) which must 
be written and may appear anywhere in the number. For example, 
0. 

3.14159265359 

-0.07 

.0000567 

The magnitude of the constant ranges from 10"^^ to 10^^ or zero. The accuracy of a real con- 
stant is eleven significant digits. 

With an exponent, the constant contains one or more digits with or without a decimal point, 
followed by the letter E, followed by the exponent. The exponent represents the power of ten by which 
the number to the left of the exponent is to be multiplied. This is called floating point or exponen- 
tial representation. 
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Examples 

5E-2 means .05 

1.E7 means 10^ 

1.973E3 means 1973. 

.0271828E+2 means 2.71828 
25E3 means 25000. 

The magnitude of the exponent is from E-75 to E+75, 
The value of zero is represented simply as 0. 
Complex Constants 

A complex constant is expressed as two single precision numbers (eleven digits of accuracy) 
separated by a comma and enclosed in parentheses . The first number represents the real part of the 
complex number; the second represents the imaginary part of the complex number. A complex number 
always is printed with a decimal point. 

Examples 

C3, 5.2) has the value of 3+5. 2i 

(-1.8, .16E2) has the value of -1.8+16i 

(2.4, 0) has the value of 2.4+Oi 

Logical Constants 

A logical constant can be either .TRUE, or .FALSE. . 
Hollerith Constants 

A Hollerith constant is expressed as nH followed by n characters, or as a group of characters 
enclosed by a pair of double quotes, a pair of single quotes, or a pair of $ signs. (The characters 
can be any character or digit except a semicolon and an up arrow.) 

Examples 

3HYES 

5HABC12 

7HAB CDE 

"XXX" 

'CODE' 

$YES$ 
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3j, Variables 

Var^iable Types 

A variable is a symbol or name which represents a quantity. The value of a variable may be 
changed throughout the program. 

In CCS FORTRAN IV the name of a variable may consist of any number of alphanumeric characters. 
The first character must be alphabetic CA through Z) . 

There are five types of variables in CCS FORTRAN IV: Integer, Real, Double Precision, Complex, 
and Logical. The variable type corresponds to the type of data the variable represents. Thus, an 
integer variable represents integer data, a real variable represents real data, and so on. A ntmiber 
stored in a Double Precision variable has at least seventeen digits of precision. 

The programmer must show the data type for each variable in his program. Data types other 
than integer and real must be declared explicitly with a "type declaration" statement. (See Section 
, Declaration Statements.) 

Integer and real data types may be declared implicitly. If the first letter of a variable 
name is I, J, K, L, M, or N, the variable is Integer. Any other variable name not appearing in any 
type declaration statement is Real. This is implicit type declaration. Explicit type declaration in 
a type declaration statement however, overrides implicit type declaration. 

Scalar Variables 

A scalar variable represents a single quantity such as: 
N 

INDEX 
ALPHA 
X12 

Array and Array Variables 

A group of variables which form or belong to a single class or collection may be related to 
one another by subscript notation. Such a collection is called an array, and the variables belonging 
to the array are called array elements. 

A string of numbers in a single row or column is thought of as a one-dimensional array. The 
elements are identified by a single subscript as shown below. 

Usual Notation FORTRAN Notation 

al ACl) 
a2 AC2) 
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ai 



ACD 



an 



ACN) 



If two subscripts are used to identify the elements of an array, the array is a two-dimensional 



array. For example, if there are three rows and four columns in a table, AC2,3) refers to the element 
in the second row and third colimn. 

In CCS FORTRAN IV, there is no limit to the number of dimensions of an array. 

The array type may be Integer, Real, Double Precision, Complex, or Logical. 

A subscripted variable represents a simple element of an array. A subscripted variable is de- 
noted by a variable name followed by an arithmetic expression in parentheses (see "Expressions" below) . 

Examples 

AC3) 

B(-5) 

CCO) 

MOVE (X, 5) 

Y(M,l,l,N+3) 

VOLTAGE C2*N+1,L,L+1) 

4j Expressions 

An expression in its simplest form may consist of a single constant, a variable, or a function 
Csee "Functions"). Also, an expression may denote a computation between two or more constants and/ or 
variables or functions. 

There are two kinds of expressions in FORTRAN IV; arithmetic and logical. The value of an 
arithmetic expression is always an integer, real, or complex number. The value of a logical expres- 
sion is either .TRUE, or .FALSE. . 

a. Arithmetic Expressions 

An arithmetic expression may consist of a single basic element; that is, a constant, variable, 

or function. 
Examples 



X 



3.14 



ACS) 

SQRT [ALPHA] 
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More complicated arithmetic expressions may be formed by using arithmetic operators which de- 
termine the computations to be performed. 

Arithmetic Operators 

+ Addition 

Subtraction 
* Multiplication 
I Division 
** or + Exponentiation 

Unary arithmetic operations; that is, operations involving only one argument or operand, also 
are possible. For example, 

-B means negative of B 

+A means A 

-C**2 means negative of 

Order of Operation 

1. Parentheses may be used to specify the order of operation in an expression. When sets of 
parentheses appear within other sets of parentheses, the innermost set is evaluated first, then the 
next set, and so on. 

2. When parentheses are not used, expressions are evaluated in the following order: 
Exponentiation C** or +) 

Unary minus (-) 

Multiplication and Division (* and 
Addition and Subtraction (+ and -) 

3. Arithmetic expressions are evaluated from left to right. Thus, when two operators of the 
same precedence appear, the left operation is evaluated first. 

Examples 

A+B*C means A+CB*C) 



A/B/C means 



A/B*C/D means 



. A+B/C**2 means A+ — - 

L C2 
B 

C 



B ^ 



(CA+B)/C)**2 means ^^j' 
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Modes of Expressions :Mixed Expressions 

The kind of arithmetic performed depends on the type or modes of the operands. For example, 
if both operands are integers, integer arithmetic is used. Thus, 3/4 causes an integer division 
and has the value zero. But 3./4. gives a result of .75 since both operands are real and division 
with real numbers is performed. 

Arithmetic expressions containing constants or variables of more than one type are called 
mixed expressions. In a mixed expression, as each operation is performed, the types of the two 
operands are compared; the smaller type is converted to the larger type and the result is of the 
larger type. 

Number 

Type Magnitude 
Complex Largest 
Double Precision 
Real 
Integer 

Logical Smallest 

Examples 

A/l The integer I is converted to real and real division is performed. 

N/CI+2.) The expression CI+2.) is first considered. Since 2. is real, I is converted to real and 

the result of the addition is real. N is then converted to real since it is to be divided 
by a real number. The result of the whole expression is real. If N=3 and 1=4, the re- 
sult is 0.5. 

A* (N/2) The result of the expression CN/2) is an integer since both operands are integers. (N/2) 

is then converted to a real number since it is to be multiplied by A which is real. If 
A=6 and N=3, the result is 6. . 

A* (N/2.) N is first converted to real since it is to be divided by a real number 2. . Real divi- 
sion is performed. Two real numbers are then multiplied together. If the values for A 
and N are the same as in the previous example, the result is 9. (3./2.=1.5 but 3/2=1). 

B+l If B is a complex variable with a value of (2.7,8.1) and 1=1, I will be converted to a 

real number (l.,0). The result will be (3.7,8.1). 

D+B If D is a double precision variable and B a complex variable, the value of D will be 

truncated to single precision, then treated as a complex number with zero as the imagi- 
nary part. 
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An exception to the conversion rule is with exponentiation. Integer expressions which appear 
as exponents are evaluated in integer mode regardless of what the other operand is. For example, a 
real operand may be raised to an integer power and the result is real. 

h. Logical Expressions 

A logical expression may consist of a single logical constant or a logical variable. The value 
of a logical expression is always a truth value, .TRUE, or .FALSE. . 

More complicated logical expressions may be formed by using logical and relational operators. 
These expressions may be one of the following forms: 

1. Relational operators combined with arithmetic expressions. 

2. Logical operators combined with logical constants or logical variables. 

3. Logical operators combined with either or both forms of the logical expressions described 
in 1. and 2. above. 

Relational Operators 

A relational operator makes a comparison between expressions. For example, the expression 
X .GT. Y means X is greater than Y. This expression has a value which at any time is either true 
or false. Comparisons may be made by means of any of the following relational operators. 
Symbol Customary Notation Meaning 

•EQ. = Equal to 

•NE. ^ Not equal to 

•LT. < Less than 

•LE. ^ Less than or equal to 

•GT. > Greater than 

•GE. > Greater than or equal to 

Examples 

A .EQ. B 
X-5 .LT. Y+4 
SQRT[BETA] .GE. ALPHAI'2 

Unlike most FORTRAN IV compilers, CCS FORTRAN IV can make comparisons between arithmetic ex- 
pressions of different types. For example, A .LE. I is valid. I will be converted to a real num- 
ber before the relation is evaluated. 

The relative magnitude of expression types and rules for comparing expressions are as follows: 
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Expression 

Type Magnitude 

Complex Largest 

Double 

Real 

Integer 

Logical Smallest 

Let El and E2 be arithmetic expressions o£ any type in the following: 

(1) Comparing with .EQ. and .NE. : 

El .EQ. E2 
El .NE. E2 

The expression type of smaller magnitude is converted to the type of larger magnitude. Com- 
parison is then made. In the case of complex expressions, both real and imaginary parts are compared, 
For example, (l.,2.) .EQ. 1 is false. The integer 1 is converted to a complex number (l.,0) and 
the imaginary part of this number is not equal to 2. . 

(2) Comparing with .LT., .GT., .LB., .GE. : 

El .LT. E2 
El .GT. E2 
El .LE. E2 
El .GE. E2 

The expression type of smaller magnitude is converted to the type of larger magnitude as 
above. Take the difference between the two expression types of the same magnitude. In the case of 
complex expressions, compare the real part of the difference between El and E2 to zero. If the dif- 
ference is positive, El is greater than E2. If the difference is equal to zero. El and E2 are equal. 
If the difference is negative, El is less than E2. For example, (l.,2.) .LT. (l.,3.) is false 
since the difference between the real parts is zero. (2., -100) .GT. (l.,4.) is true since the 
difference between the real parts is positive. 

Logical Operators 

A logical operator operates on logical expressions. The result of a logical operation is 
.TRUE, or .FALSE. . 

There are four logical operators. In the following table where X and Y are logical expres- 
sions, the results of the logical operations are shown. 
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OPERATOR 


X 
Y 


.TRUE. 
.TRUE. 


.TRUE. 
.FALSE. 


.FALSE. 
.TRUE. 


.FALSE. 
.FALSE. 


AND 


X .AND. Y 


.TRUE. 


.FALSE. 


.FALSE. 


.TRUE. 


Inclusive OR 


X .OR. Y 


.TRUE. 


.TRUE. 


.TRUE. 


.FALSE. 


Exclusive OR 


X .EOR. Y 


.FALSE. 


.TRUE. 


•TRUE. 


.FALSE. 


NOT 


.NOT. X 


I£ A is .TRUE., then .NOT. A is .FALSE.; 
i£ A is .FALSE., then .NOT, A is .TRUE. . 



Examples Of Logical Expressions With Logical Operators 
A .AND. B. where A,B are logical variables. 
(X .EQ. 3) .OR. Y .LE. 4 
(E*5 .GT. A-B) .OR. CE .LE. 100) 
(A/2 .LT. ,SQRT[5]) .AND. I .EQ. J 
.NOT. W .AND. .NOT. L 
(A .GT. 100) .EOR. CB -GT. 200) 

Order of Computation In Logical Expressions 

Where parentheses are omitted, the order in which operations are performed is as follows: 
Evaluation of Functions (for example, SQRT) 
Exponentiation (** or f) 
Unary minus (negative) 
Multiplication and Division (* and /) 
Addition and Subtraction (+ and -) 

Relational operators (.EQ., .NE., .LT., .LE., .GT., .GE.) 

.NOT. 

• AND. 

.OR. 

.EOR. 

For example, the expression SQRT[A] .GT. C+3 .AND. .NOT. R .OR. S is evaluated in the 
following order: 



1. 


SQRT [A] 


call 


the result 


V 


(function) 


2. 


C+3 


call 


the result 


W 


(exponent i at i on) 


3. 


V .GT. W 


call 


the result 


X 


(relational operators) 


4. 


.NOT. R 


call 


the result 


Y 


(highest logical operator) 


5. 


X .AND. Y 


call 


the result 


Z 


(second highest logical operator) 


6. 


Z .OR. S 


final 


operation 







Parentheses may be used to modify the order of evaluati 
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B. ARITHMETIC AND LOGICAL REPLACEMENT STATEMENTS 



General 

A replacement statement uses the replacement operator "=". A replacement statement has the 
form V=e meaning to replace the current value o£ the variable V with the current value of the 
expression e. The "=" means "is replaced by" rather than "is equivalent to". 

Ij. Ari thmeti c Repi acement Statements 

The general form of an arithmetic replacement statement is variable = expression where 

V is a scalar or a subscripted variable and e is an arithmetic expression. Execution of the statement 

causes the value of the variable to be replaced by the value of the expression. If the mode of the 

expression differs from the mode of the variable, the mode of the expression is converted to the mode 
of the variable before replacement takes place. 

Examples 

In the following, V, A, and B are real variables; I and J are integer variables; C is a com- 
plex variable; and D is a double precision variable. 

VAR=A*B/5 The current value of VAR is replaced by the result of A*B/5. 

A=I The value of I is converted to a real value and this result replaces the current value 

of A. 

1=1+1 The value of I is replaced by the value I+l. 

A=C The real part of the complex variable C replaces the value of A. 

C=A The value of A replaces the value of the real part of the complex variable C, the 

imaginary part of C is set to zero. 

C=ItJ I is raised to the power J and the result is converted to a real value which replaces 

the real part of the complex variable E. The imaginary part of C is set to zero. 

D=J The value of J is converted to double precision and stored in D. 

C=C3.4,2.0) The value of C is replaced by the complex constant (3.4,2.0). Note that CC=A,B) where 
A and B are variables is not allowed. 

3. Logical Replacement Statements 

The general form of a logical replacement statement is variable = expression where the 
variable is a scalar or subscripted logical variable, and the expression is a logical expression. 
Execution of the statement causes the variable to be replaced by a value of .TRUE, or .FALSE, depend- 
ing on whether the expression is true or false. 
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Examples 

In the following, H, K, and M are logical variables and I is an integer variable. 
H=.TRUE. The value of H is replaced by the logical constant .TRUE. . 

K=.NOT. M If M is .TRUE., the value of K is replaced by the logical constant .FALSE. . If 
M is .FALSE., the value of K is replaced by the logical constant .TRUE. . 

H=3 .EQ. I The value of I is converted to a real value. If the real constant 3. is equal to 
this result, the logical constant .TRUE, replaces the value of H. If 3. is not 
equal to this result, the logical constant .FALSE, replaces the value of H. 

NOTE: If a logical expression is used in an arithmetic assignment statement, the logical ex- 
pression has an integer value of 1 if it is .TRUE., and an integer value of zero if it is .FALSE. . 
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C. CONTROL STATEMENTS 

2j General 

Normally, FORTRAN statements are executed sequentially. That is, after one statement has been 
executed, the statement immediately following it is executed. Control statements allow the programmer 
to alter and control the normal sequence of execution of statements in a program. When it is neces- 
sary to alter this normal sequence, certain statements are labeled so that they may be referred to by 
control statements. A statement label can be any integer number. Any statement except declaration 
statements can be labeled. 

2. IF Statements 

Two types of IF statements are provided to make decisions; namely, logical IF and arithmetic 
IF. With a logical IF, decisions are based on a logical quantity being true or false. With an 
arithmetic IF, decisions are based on an arithmetic quantity being less than zero, zero, or greater 
than zero. 

Logical IF Statements 

The logical IF statement has the general form IF (expression) Statement where the ex- 
pression is any logical expression and the Statement is any executable statement except DO or another 
logical IF statement. 

If the expression is true, the statement is executed; execution then goes to the next state- 
ment except if the statement is a transfer statement. If the expression is false, the statement will 
not be executed; execution simply goes to the next statement. For example, in the following statements 
IF(X .GT. Y) A=B 
C=A*5 

if the value of X is greater than Y, the value of A will be replaced by the value of B; then C is set 
to A*5. If X is not greater than Y, A will not be set to the value of B. C will be replaced by A* 5 
whatever the value of A is . 

Arithmetic IF Statements 

The arithmetic IF statement has the general form IF (expression) nl, n2, n3 where the 
expression is any arithmetic expression (except complex) and nl, n2, and n3 are statement numbers. 
This statement causes control to be transferred to the statement numbered nl, if the value of the ex- 
pression is negative, to n2 if it is zero, and to n3 if it is positive. 

Examples 

IF (K-N) 10, 10, 20 If the expression K-N has a value which is less than or equal to 

zero, indicating that K is less than or equal to N, the next 
statement executed is statement 10; otherwise if K is greater 
than N, the next statement executed is statement 20. 
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IF (Y-ACI)) 10, 15, 15 



IF CXCI,J)*3f2) 12, 5, 30 



30 



12 



If Y is less than ACI), that is, the value of Y-A(I) is less 
than zero, the statement executed next is statement 10. If the 
value for Y-ACi) is zero or positive, the statement executed 
next is statement 15 . 

If the value of the expression CXCl,J)*3t2) is negative, state- 
ment 12 is executed next. If the value of the expression is 
zero, statement 5 is executed next. If the value of the ex- 
pression is positive, statement 30 is executed next. 



:=I+J 



=D-B 



= (F*B)/C 



3. 



DO Statements 



In a majority of numerical computation and information processing procedures, repetitive 
processes are used. Computer programs almost always involve a group of steps that need to be executed 
repeatedly. The iteration or DO statement controls repetitive processes. 

The statement DO 30 J=l, 10, 1 controls the repetition of all succeeding statements down 
through and including the statement labeled 30. Repetition of these statements is controlled by 
varying an index called J from an initial value of 1 to a terminal value of 10 in increments of one. 

The general form of a DO statement is DO n I=ml, m2, m3 
where n is the statement number of the last step in the repetition. 

I is the index variable whose value changes during the repetitive process. The name may be 

any non-subscripted variable, 
ml is the initial value for the index variable. 

m2 is the upper limit for the index variable; the index variable must exceed this value to 
terminate the repetition. 
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in3 is the increment in the index variable for successive repetitions. (m3 may be omitted, in 
which case the increment is assumed to be one.) 

NOTE: The indexing parameters ml, m2, and m3 may be any arithmetic expression. The increment 
m3 can be negative. 

The group of statements starting with the DO statement, through and including the statement 
labeled n, is called a DO loop. 

Example 1 

In this example, the 'increment' is negative. The value of I is printed beginning with 3, by 
steps of -1, ending with 0. 
>10 DO 30 1=3,0,-1 

>20 30 DISPLAY I 

>30 END 

>RUN 
3 
2 
1 


STOPPED AT: 20. 

> 

Example 2 

In this example, statement 100 is executed beginning with X=l, by steps of .2 and stops when 
X exceeds 1.5. 
>10 DO 100 X=l,1.5,.2 

>20 100 DISPLAY X 

>30 END 

>RUN 
1. 
1.2 
1.4 

STOPPED AT: 20. 

> 
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Throughout the DO loop, the index variable is available both in subscripts and as an ordinary- 
integer variable. For example, the statements 

DO 100, 1=3,20 
100 ACI)=I 

will set A(3) to a value of 3, AC4) to a value of 4, up to AC20)=20. 

NOTE: The indexing parameters of a DO statement il, ml, m2, m3) may not be changed by a 
statement inside a DO loop. 

Example 

The following DO loop is used to find the total of N numbers of the array X. 
SUM=0 

DO 25 1=1, N 
25 SUM=SUM+XCI) 

A CONTINUE statement may be used as the last statement for a DO loop. The above may be written 
as DO 25 1=1, N 

SUM=SUM+XCI) 
25 CONTINUE 

The CONTINUE statement causes no action; it merely serves as a dummy statement to refer to the 
end of the loop. 

Nested DO Loops 

A DO loop may include other DO loops provided that the range of each inner loop is contained 
completely within the range of the outer loop. That is, the two DO loops may not partially intersect 
one another. These loops are called nested loops. The following skeleton examples illustrate this: 
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Allowed 

DO 10 1 = 1, N 

DO 20 J=1,M 



20 CONTINUE 
-10 CONTINUE 

Not Allowed 
.DO 10 1=1, N 
•DO 20 J=1,M 



•-10 CONTINUE 



20 CONTINUE 



Allowed 

DO 15 K=1,N 



DO 5 A= 1,7. 5,. 5 



5 X=7+Y 



DO 7 BAJ=2,3 



DO 100 B=A*C,2 



•-100 CONTINUE 



7 X=22 



■15 CONTINUE 



Nested loops may end with the same terminal statement. The following statements sum up an 
array of three rows and four columns. 
TOTAL=0 . 
DO 15 1=1,3 
DO 15 J=l,4 
T0TAL=T0TAL+ACI,J) 
15 CONTINUE 

While it is permissible to transfer control out of a DO loop, transfer of control into a loop 
is not allowed. The following examples illustrate this: 
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Allowed 



Not Allowed 



DO 200 K=5,M,2 
C=K*P-Q 

IF(C .GT. 100) GO TO 50 



5 --- 



DO 10 1=1, N 



200 CONTINUE 



10 CONTINUE 



GO TO 5 



50 



One exception to this rule is when returning to a DO loop from an extended range. An extended 



range o£ a DO loop is statements outside the loop to which a statement within the DO loop transfers. 

They could be included as part of the loop but are written outside since they were to be executed 

again in the program. For example, in the following statements 

DO 100 K=l,10 

X=A*K 

Y=B*K+2 

GO TO 50 

60 C=X-Y 

-100 CONTINUE 

50 X=SQRT[X*2] 

Y=SQRT[Y*3] 
GO TO 60 

the statements 

50 X=SQRT[X*2] 

Y=SQRT[Y*3] 

are in the extended range of the DO loop. Thus, GO TO 60 is permitted since this returns control to 
the DO loop from an extended range of the loop. 

Implied DO Loops 

An implied DO loop may be used in an input, output, or a DATA statement. 
For example, the statement 



ACCEPT CA(I), 1=1, 10) 
reads in the data values for A(l) to ACIO) . 



The statement 



DISPLAY CB(I) ,1=1,10,2) 
prints the data values for A(l), AC3),..., up to AC9). 
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The statement 

DATA CACI),I=1,3)/1,2,3/ 
initializes the ACl) to 1, AC2) to 2, and AC3) to 3. 

4; GO TO Statements 

The GO TO statements permit the user to alter the sequence in which program statements are 
executed. The GO TO statements transfer control to the statement number specified in the GO TO state- 
ment. There are three types of GO TO statements: 

1. Unconditional GO TO statement. 

2. The Computed GO TO statement. 

3. The Assigned GO TO statement. 

Unconditional GO TO Statements 

The unconditional GO TO statement causes an unconditional transfer. The form of the statement 
is GO TO statement number 

Example 

GO TO 20 

The GO TO statement alone is unconditional. However, the GO TO statement may be written so 
that execution of it is conditional. For example the statement IF (X .EQ. Y) GO TO 30 
will cause transfer to the statement labeled 30 if the values of X and Y are equal. If X is not equal 
to Y, the next statement in sequence will be executed. 

The following statements illustrate the use of GO TO statements, kl, k2 --- kn is a list of 
integers. The following statements count the number of positive even and the number of negative even 
integers. Zeros are ignored. 

NPOSEVEN=0 

NNEGEVEN=0 

DO 10 1=1, N 

IF (KCI) .NE. .AND. (K(I) -K (I) /2*2) .EQ. 0) GO TO 20 
GO TO 10 
20 IF CKCI) .GT. 0) GO TO 30 

NNEGEVEN=NNEGEVEN+ 1 

GO TO 10 
30 NP0SEVEN=NP0SEVEN+1 
10 CONTINUE 
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Computed GO TO Statements 

The computed GO TO allows transfer to one of the several places in the program, depending on 
the integer value of a variable. The form of the statement is GO TO (nl ,n2 ,n3, . . . ) I 
where I is a non- subs crip ted integer variable and nl, n2,... are statement numbers to which the pro- 
gram will transfer depending on the value of I. If the value of I is 1, the program will transfer to 
the statement labeled nl; if the value of I is 2, the program will transfer to statement labeled n2, and 
so on. For example, GO TO C60,70,85)K will transfer control to statement 60, if the value 

of the expression K is 1, to 70 if the value is 2, and to 85 if the value is 3. 

If the value of K is not an integer, it will be truncated to an integer value. If K is nega- 
tive or greater than the number of statement numbers, an error will be indicated. 

Assigned GO TO Statements 

The assigned GO TO statement is useful for conditional branching. The form of the statement 
is GO TO variable 

This statement transfers control to the statement whose number was last assigned to the variable by 
an ASSIGN statement. The variable must appear in some previously executed ASSIGN statement. 
Examples 

GO TO L 

GO TO EX (3) 

The variable in an assigned GO TO statement is a control variable and has a statement number 
(a label) as its "value" and not a numerical quantity. 

5^ ASSIGN Statement 

This statement assigns a statement number to a variable for a subsequent assigned GO TO state- 
ment. The form of the statement is ASSIGN integer TO variable where integer is the statement 
niomber to which control will be transferred by the assigned GO TO statement. 

In the following 



ASSIGN 10 TO R 

IF C(A .LT. 100) .AND. CB .LT. 50)) ASSIGN 20 TO R 
A=450 
B=-200 
GO TO R 
10 X=CA+B)*2 
Y= (A-B)*X 
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GO TO 50 
20 X=3-B*A 
Y=3+B*A 
GO TO 50 

50 . 

statement label 10 is assigned initially to R. If the value of A is less than 100 and B less than 50, 
R is reassigned to statement 20. When the statement GO TO R is executed, control will be transferred 
to statement 10 or statement 20 depending on the value of R. 

Another Example 



ASSIGN 100 TO ACl) 
ASSIGN 200 TO AC5) 
ASSIGN 450 TO AC9) 
I=M*N-J 
GO TO A(I) 

100 . 
200 . 
450 . 

When the statement GO TO A(I) is executed, the program will go to statement 100, 200, or 450 
depending on whether the value of I is 1, 5, or 9 respectively. 

6. CONTINUE Statement 

The form is 

CONTINUE 

The CONTINUE statement is used as a dummy statement which may be placed anywhere in the pro- 
gram without affecting the sequence of execution. 

It may be used as the last statement of a DO loop. A CONTINUE statement must be used to avoid 
ending the DO loop with a GO TO, PAUSE, STOP, RETURN, Arithmetic IF, or another DO statement. 
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Example 1 

In the following statements which count the number of non-zero elements of the array A, the 
CONTINUE statement is used to skip the count if an element is zero. 

NCOUNT=0 

DO 15 1=1, N 

IF (A(I)-1)10,15,10 
10 NC0UNT=NC0UNT+1 
15 CONTINUE 

Example 2 

DO 40 1=1,20 
8 IF CX(I)-YCI))5,40,40 
5 X(I)=X(I)+1 
Y(I)=Y(I)-2 



GO TO 8 
40 CONTINUE 

In the above example, the CONTINUE statement is used as the last statement of the DO loop to 
avoid ending the loop with the statement GO TO 8. 

7. PAUSE Statement 

The form is 

PAUSE 

A PAUSE statement may be placed anywhere in a program. When this statement is executed, the 
running program will pause, return to the CCS command mode and print the message PAUSE. The computer 
will then type the line number at which the pause occurred and a >. Direct statements and CCS com- 
mands may then be executed. Execution of the program may be resumed at the point of interruption by 
typing CONTINUE. 

8; STOP Statement 

The form is 

STOP 

When a STOP statement is encountered, execution of the program is terminated. The user may 
enter statements for direct execution but the program cannot be continued. 



44 



9. END Statement 

The form is 

END 

CCS FORTRAN IV programs and subprograms must end with an END statement. The END statement 
may not be labeled with a statement number; that is, it is not possible to transfer to this statement 
from other parts of the program. 
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D. INPUT/OUTPUT STATEMENTS 



GENERAL 

Input and output statements are, as their name implies, the statements used to transfer (com- 
municate) data to and results from the computer. Input and output may be from the terminal or from 
data files stored on the disk. Terminal input and output always is symbolic and may be in free format 
(ACCEPT, DISPLAY) or in formatted (READ, WRITE, FORMAT) form. 

The data files may be in either symbolic or binary form. Data files always must be opened 
before use and closed after use (OPEN, CLOSE). Symbolic files always are input and output in for- 
matted form (READ, WRITE, FORMAT). Binary files always are input and output in a format free form 
(READ, WRITE). 

FREE FORMAT TERMINAL INPUT AND OUTPUT - ACCEPT, DISPLAY 

The simplest way to input and output data on the terminal is to use the free format input/out- 
put statements ACCEPT and DISPLAY which are unique to the Tymshare system. The form of the input com- 
mand is as follows: 
ACCEPT list 

Example 

ACCEPT A, B, (C (I), 1=1, 9) 
Input List 

The input list may include any legal FORTRAN IV variable and/ or array name. As seen in the 
example above, an implied DO loop may be used to input part of an array. For example, the implied DO 
loop (C(I), 1=1,9) will cause values for the variable C(l) through C(9) to be requested. An entire 
array which has been dimensioned previously in the program may be requested by its name. If, for 
example, B in the above statement had been dimensioned B(15), the statement ACCEPT B would request the 
array elements B(l) through B(15). 

Input In Response To An ACCEPT Statement 

When the ACCEPT statement is executed, the system will ring a bell and wait for the user to 
fill the variables and/or arrays specified. The user must then type in one value for each variable 
and/or array element specified in the variable list. The user may use a Carriage Return, a Line Feed, 
a comma, or a space bar to terminate each value input. The value supplied will assume the type (mode) 
of the variable or array in which it is stored. Unless otherwise specified in a TYPE statement, all 
variables and arrays beginning with I through N will store integer values; all other variables and 
arrays will store real values. 
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The DISPLAY command is used to output data to the terminal. The statement form of the DISPLAY 
command is as follows: 
DISPLAY list 

Example 

DISPLAY A, B, CCCI), I = 1,9), 5*6 + 5, SQRTlABS [-25] ] 

Variables, arrays, and expressions may be included in the output list. The DISPLAY command 
will cause the values stored in the variables and/or arrays specified to be printed on the terminal. 
Any expressions included in the list will be evaluated and their values returned. 

Literal Text In The I/O Lists 

An extension of the ACCEPT and DISPLAY statements allows the user to write literal text on the 
terminal either alone or in conjionction with variable values. The text desired may be enclosed in 
either single quotes, double quotes, or dollar signs as follows: 

ACCEPT 'TEXT' 
or "TEXT" 

DISPLAY $TEXT$ 

Examples 

DISPLAY $THIS PROGRAM COMPUTES THE . . .$ 
ACCEPT "THE ORIGINAL VALUE IS = $" , C 
DISPLAY 'THE MAXIMUM CAPACITY = ' , T, "TONS" 

FORMATTED INPUT AND OUTPUT - READ. WRITE, FORMAT 

Data also may be input and output in formatted form using the ASA Standard FORTRAN IV state- 
ments READ, WRITE, and FORMAT. The READ and WRITE statements are the executable input/output state- 
ments. The FORMAT statement is a non- executable reference statement which supplies certain information 
about the size and mode of the data values being read or written. The READ, WRITE, and FORMAT state- 
ments are used both for terminal input and output and for symbolic disk file input and output. The 
only difference between formatted terminal input/output and formatted symbolic disk file input/output 
is that symbolic disk files must be opened before use and closed after use. The general form of the 
input/output statements is : 

READ 

(file format 
number number 

The format number refers to the label of the FORMAT statement used to read or write the values of the 
variables and/ or arrays in the list. 



list 
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The file number indicates the file from which the data is to be read or on which it is to be 
stored. If a disk file is being used, the file number specified is the one used in the OPEN statement. 
The file number in this case may be 2, 3, or 4. If the terminal is being used, the file number will 
be or 1. The file nxjmbers and 1 are reserved for terminal input and output. The file number 
specifies terminal input; the file niimber 1 specifies terminal output. 

The input list may contain variables and arrays and implied DO loops. In addition to these, 
the output list also may contain expressions and constants (ii^cluding Hollerith constants "ABC", 
'ABC, $ABC$, 3HABC). 

Examples 

READ (3,7) A, I, CCCI), I = 1, 10) 

WRITE CI, 90) A, CBCR) R = 200,300,5), "5TG", SQRT[79] 

FORMAT Statement 

The FORMAT statement is a non-executable statement which specifies the size and type of the 
data values being read or written. The FORMAT statement consists of a set of field specifications 
which indicate how each value is to be read or written. The form of the FORMAT statement is: 

number FORMAT (field spec^, field spec2,...) 

Example 

90 FORMAT (15, 2X,F7 . 2 ,7 (3 (14) , 2 CF9. 3) ) , A5) 
Field Specifications 

The field specification supplies certain information about the size and mode of the value be- 
ing read or written. The information supplied by the field specification is used somewhat differently 
for input than for output. 

During Input . The size specification is crucial since it indicates the number of characters 
to be read from the file. The mode of the field being read is occasionally completely ignored (for 
example, an exponential value may be read using a decimal field specification) . The mode of the vari- 
able into which the value is being read overrides the mode specified in the FORMAT statement. 

During Output . The type of the field specified always determines the mode of the value output. 
For example, an F field will always output the value as a decimal; an I field will always output the 
value as an integer. The size specification indicates the number of spaces that will be used to print 
the number. The decimal places will be supplied as indicated in the FORMAT statement. 

NOTE: The field size is specified differently for the different types of fields. Careful at- 
tention should be paid to exactly how the field size is specified and what is included in the size 
specification. 
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Quick Reference Chart— Field Specifications 



Field Type 


Manner Specified 


Usage 


I 


rlw 


Integer field C123) 


F 


rFw.d 


External fixed point decimal Cl«23) 


E 


rEw.d 


Floating point C1.E09) 


J 


rJi.d 


Floating point C1.E09) 


D 


rDw.d 


Double precision C1.D09) 


G 


rGw.d 


Generalized (for E formats) 


L 


rLw 


Logical (T or F) 


A 


2]*Aw 


A 1 r\haniim<=>r'i <- rTnMP<?"\ 

J. UilclllLullCX V^UUINILO J 


H 


wHs 


Hollerith C3HEND) 


ti 


"s" 


("END") 




•s' 


C'END') 


$ 


$s$ 


C$END$) 


X 


wX 


Spacing - spaces w times 


T 


Tw 


Tab (spaces to column w) 


P 


£P 


Scaling 


/ 


/ 


Generates a Carriage Return 


Z 


Z 


Suppresses a Carriage Return 



Symbols 

w - field width (entire number of characters required) 

d - number of decimal digits 

i - nuHtoer of integer digits 

s - string of characters 

f - power of 10 

r - repeat count 

Numeric Fields 

I^ The integer field is represented by the letter I. The general format for an integer field 

specification is Iw where w represents the maximum length of the number which can be input or output. 
If the number is negative, one extra digit should be specified for the negative sign. 

Input 

If the entire field is filled on input; for example, 123 read with 13, the number being read 
must be an integer Cdecimal points, E's, and D's will not be allowed in the field). 
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However, if the number is terminated by a comma before the entire field is filled, free form 
input is assumed and any type of niimber Csuch as 1.5, 2E(D7) may be input. In this case, all charac- 
ters in the number will be counted. For example, -IE 07 has 6 characters so, if a comma terminator 
is used, a minimum field specification of 17 would be required. An 16 field would be filled com- 
pletely. 

Output 

Any number output using an I field specification will be output as an integer. Any decimal 
digits will be truncated. If the field width w is larger than the number, leading blanks will be sup- 
plied. If it is smaller, an error diagnostic will be given. 

F_ The external fixed point field (decimal) is represented by the letter F. The general form of 

the F field is Fw.d where w represents the entire width of the field which includes the plus or minus 
sign, the total number of characters in the field, the decimal point and any blanks, and the d indi- 
cates the number of decimal digits in the number. 

Input 

Any numeric constant (integer, real, double precision)! ju^y jjg read into a variable or array 
element using an F format. The w indicates the number of characters to be read, and the d indicates 
the number of decimal positions. If, for example, the niamber 12345 were read using an F5.2 field 
specification, the number would be stored in the variable as 123.45. If the number being read has a 
decimal point in it, the d specification will be ignored. Thus, the number 123.4 read with the format 
F5.2 would be stored as 123.4. 

Output 

Numbers output using an F field will be truncated to d decimal digits. If the width w is 
larger than required, the number will be filled in with leading blanks. For example, the number 
-1.59764 would be output as _- 1.597 using the field specification F7.3 and as _- 1.597640 using a 
field specification F10.6. If insufficient width (w) is specified, an error message will occur. 

E^ The floating point or exponential field appears in the following general form, Ew.d where w 

represents the entire width of the field including the plus or minus signs, the decimal point, all 
numeric characters, and the exponent; and d indicates the number of decimal digits . 

£ The double precision or D field is specified in exactly the same manner as the E field except 

that the exponent is written with a D instead of an E; that is, Dw.d. 

Input 

Any numeric constant may be input using the E and D format specifications . 1 The number of 
characters specified by the w will be read and d decimal digits assigned. For example, the number 
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+12345E+07 would be read, using the format E10.2, as 123.45E+07. If the number being read has a deci- 
mal point in it, the d specification will be ignored. Thus, -1.23D-05 would be read by D10.2 as 
-1.230D-05. 

Output 

Numbers output using an E or D format will be truncated to d decimal places. When the E field 
specification is used, the number will be output in E notation C1.5E+02). When the D field specifica- 
tion is used, the number will be output in D notation C1.5D+02). 

The J field also is used with floating point numbers. The general form of the J format, how- 
ever, is Ji.d where i represents the number of integer digits and d represents the number of decimal 
digits. 

Input 

Any numeric constant may be input using a J field specification. 1 When the field is input, i 
integers will be read followed by d decimal digits. If the number is input in exponential form, the 
exponential part of the number (E±07) is not specified in the field specification. For example, the 
number 12345E+04 would be read by J3.2 as 123.45 E+04. If the number being read has a decimal point, 
the decimal position specified by d will be ignored. For example, the number 123.45E-7 would be read 
by J4.3 as 123.450E-07. 

Output 

Numbers output with a J format will contain i integer digits, d decimal digits, and the ex- 
ponent. The decimal digits will be truncated if necessary. For example, the number 123.45 E+07 will 
be output as 123.4 E+07 by the format J4.1. 

G_ The G, or generalized field, is a combination of the F and E fields. The general form is Gw.d 

where w represents the total field length and d represents the number of decimal digits. 

Input 

Same as E and D field specification. Any numeric constant may be input; w characters will be 

read. 

Output 

The G field specification will output the number in decimal or exponential form depending upon 
its size. 



A complex number may be read also. Complex numbers require a separate field specification for the 
real and for the imaginary parts of the number. Thus, the number 2.3, 1E07 could be read using the 
format- statement FORMAT (F4. 1,E6.0) 
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The number (n) will be output in F format Cwithout an exponent) if it is within the range 
0.1 < n < 10'^, otherwise the number will be output in E format. For example, using the specification 
G12.3 the number 244. will be printed as 244. The number 2444. will be printed as .244E+04. 

Scaling - P 

A scaling factor may be used in a FORMAt statement to alter the normal form of the E, D, G, 
and F fields. The scaling notation is in the form fP where f indicates a power of ten by which scal- 
ing is to occur. Thus, IP would indicate multiplication by ten; 2P multiplication by 100 (10 ); -IP 
would indicate multiplication by .1,(10"^). 

Scaling may be. set and reset anywhere in the FOK^T statement. Once scaling is set it will 
be used throughout the program until it is reset. OP (10 to the power, or 1) may be used to 
set the scaling back to normal. 

F0R^4AT (1P,I4,F7.9,0P,F5.7,-2P,E7.6) 

Input 

When a number is input as a decimal in any of the above formats, the actual value of the num- 
ber will be altered. If the number is input in exponential (lEO?) or double precision (1D07) form, 
it will not be affected. 

Value Input Scaling Factor Value Stored 

109.3 IP 1093. 

1.2E05 IP 1.2E05 

111.3 -2P 1.113 

Output 

During output scaling is determined by the field specification used. If an F field specifica- 
tion is used, the actual value of the number will be altered. When E and D field specifications are 
used, the decimal point and the exponent will be realigned, but the actual value of the number will 
not be altered. When a G field specification is used, the decimal point and exponent will be re- 
aligned if the number is output in exponential form. If the number is output in decimal form, it will 
not be affected by the scaling. 

Format Normal Output Scaled Format Scaled Output 

F4.0 15. 1PF4.0 150. 

E12.5 .15000E 02 1PE12.5 1.50Q0E 01 

D12.5 .15000D 02 -1PD12.5 .01500D 03 

G3.1 1.5 -1PG3.1 1.5 
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Non-Numeric Fields 

k 

The logical field L is used with the logical values TC.TRUE.) and FC. FALSE.). It appears in 
the form Lw where w represents the number of characters in the logical field. Whenever a logical 
value is read, the first character read must be either a T or an F or an error message occurs. If a 
T is encountered, the logical variable will be assigned a value of .TRUE.; if an F is encountered, the 
logical variable will be assigned a value of .FALSE. . When a logical value is written, a T or an F 
always will be output. If the output field size is greater than one, leading blanks will be supplied 
automatically to fill the field. 

A ■ 

The alphanumeric or A field appears in the general form Aw where w is the number of characters 
in the field. Any group of characters may be stored in or read from a variable using the A format. 
For example, character set J.J.JONES, 23. could be stored in a variable which is read or written 
using the field specification A14. 

NOTE: Do not read or write numeric values using an A format. Use one of the numeric formats. 
H (Hollerith Field) 

The Hollerith field which appears in the general form wH causes the niomber of characters 
specified by w, which immediately follow the H, to be written on a file or the terminal during output. 
During input the number of characters specified Cw) will be skipped. For example, the format 
FORMAT C8HTYMSHARE) would print TYMSHARE during output and skip 8 characters during input. 

II I d; 
> > 4* 

A Hollerith field also may be specified by enclosing it in double or single quotes or dollar 
signs. The general form is "s", 's', or $s$ where s indicates a string of characters. The string 
will be printed on the terminal or output to a file during output. During input the number of char- 
acters specified in the string will be skipped. For example, FORMAT ("THE ", 'BIG $DOG$) would 
during input, skip the first 11 characters, but during output would print THE BIG DOG. 

Spacing 
X 

The X specification, which appears in the form wX, causes w blanks to be printed during output 
or w columns to be skipped during input. For example, 3X will leave three blank spaces during output 
or not read 3 characters during input. 
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I 

The T specification appears in the form Tw. Using the first column as a reference point, it 
spaces to the column specified by w. Thus, T40 will start reading or writing the next field in column 
40. The T specification may be used to backspace during output but not on input. For example, the 
statements 

WRITE (1,1) A, B 

FORMAT CT50,I6,T10,F5.2) 

would output the value stored in A starting in column 50 and then backspace to column 10 and output 
the value stored in B. 

Generating And Inhibiting A Carriage Return - /,Z 

A slash (/) is used in a FORMAT statement to generate a Carriage Return, which causes both 
reading or writing to skip to the next record. Multiple slashes may be used to skip more than one 
record when reading a file, or to generate blank records when writing a file. This allows the user 
to include the fields specifications for more than one record in a single format statement. 

If A = 111, B = 2.135E06, and C = 5.794, the statement 
WRITE Cl,100) A,B,C, "ABCF" 
100 FORMAT C13,3X,E10.3/,F5.3,5X,A4) 
will output the following: 

111 .213E+07 Crecord 1) 

5.794 ABCF (record 2) 

NOTE: When a slash is used between two field specifications, a comma need not be used to 
separate the specifications; for example, FORMAT CI5,F6.3/G8.0) 

A Z may be used at the end of a FORMAT statement to suppress the Carriage Return that is 
normally generated each time the FORMAT statement is scanned. This allows writing on or reading from 
a single record using more than one FORMAT statement. (This feature also is very useful in document- 
ing formatted terminal input.) 

Example 

The statements below 
WRITE (1,100) 
100 FORMAT ('VELOCITY=' ,Z) 

READ (0,200) V 
200 FORMAT (F6.2) 
will print 

VELOCITY = 
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and then wait on the same line for the user to input the value of V. CThe Carriage Return that would 
normally be generated after 100 FORMAT was scanned was suppressed by the Z.) 

Repeating A Field Specification 

If the same field specification is to be used a number of times, a repeat count may be used. 
The repeat count appears immediately preceding the field specification and specifies the number of 
times the field specification is to be read or written. For example, the format statement 

FORMAT C3I4,7F9.6) 

will read or write the first three values as four digit integers and the last seven values as nine 
place fixed point values with six decimal places. 

Parentheses may be used to indicate the repetition of a series of field specifications. The 
format statement 

FORMAT (2(213), 3(14,15)) 
will read or write the first four values as three digit integers (2(213)); it will then repeat the 
specifications (four digit integer, five digit integer) three times. 

Repeating FORMAT Statements 

A FORMAT statement will be rescanned automatically whenever the number of items specified in 
the input or output list exceeds the number of fields specified in the FORMAT statement. Each time 
the FORMAT statement is scanned a Carriage Return will be generated and reading or writing will be 
continued on the next line (record) . 

Data Field/Field Specification Compatibility 

The term data field refers to a single value or string of characters. Every data field must 
be read or written using the corresponding field specification in the FORMAT statement and must be 
stored in an array element or variable. 

An input data field may appear in one of three forms: 

1. Field filled--no decimal point. In this case the data appears as a solid string of char- 
acters not including decimal points. (NOTE: Blanks are treated as zeros in numeric 
fields.) The number of characters read in the field is the same as the field width speci- 
fied in the field specification. A data field which fills the field specification entirely 
is assigned its value according to the field specification. 
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Data Field 


Field Specification 


Value Assigned 


12345 


F5.2 


123.45 


123E04 


E6.1 


12.3E04 


123 


13 


123 


ABCDE 


A5 


ABCDE 


12345 


J1.4 


1.234 


TRY 


L3 


T (true) 



2. Field filled--with decimal points. In this case the data also appears as a solid string 
of characters except that a decimal point may be included in a numeric value. The decimal 
position specified in the input field overrides any decimal position specification given 
in the field specification. The number of characters specified in the field specification 
is the number of characters that will be read. (The decimal point is counted as a char- 
acter.) 

Data Field Field Specification Value Assigned 

1.2345 F6.2 1.23 

123. E04 E7.2 123. E04 

.123 J1.4 0.123 

3. Field not filled--terminated with comma. A comma may be used to terminate the input be- 
fore the entire field specified has been filled. Whenever a comma is encountered in a 
data field, reading is discontinued and blanks are supplied to fill the field. If a blank 
field is terminated by a comma, a value of 0.0 is assigned. A field terminated by a comma 
will be read in free-format form; that is, the number will be stored exactly as it appears; 
the field specification will be ignored entirely. (Remember- -the number of characters in 
the field must be fewer than the number specified as the field width.) Using the free-^ 
format comma field terminator relieves the user of the time-consuming task of filling all 
of the fields, thus minimizing the time required to create the file and also the space 
required to store the file. One exception exists: Alphanumeric (A) fields are not ter- 
minated by a comma since a comma is considered a valid character in an A field. However, 
an alphanumeric field may be terminated by a Carriage Return or a Line Feed. (NOTE: A 
Carriage Return or Line Feed also terminates a record.) 

Data Field Field Specification Value Assigned 

33.1, F3.3 33.1 

44.2, 15 44.2 (if read into a real vari- 
44, F7.2 44, 

, F5.3 0.0 
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If the data field and the field specification are not compatible; that is, not the same size 
and in some cases not the same type Cmode)! the data may be read or written improperly or an error 
message returned. If for example, the input data field is shorter than the one specified in the field 
specification using filled field data, only part of the data value will be read. The rest of the data 
value will be read as part of the second value. Conversely, if the input data field is longer than 
the field width specified, only the first portion of the data value (w characters) will be read into 
the value; the rest of the field will be treated as part of the next field. 

Example 

A should be 12.34; B should be 5.12. When the field width was counted the decimal point was 
forgotten. 

The file appears as: 12.345.12 and is read by: 
READ (3,1) A,B 
1 FORMAT(F4.2,F3.2) 

As: 

A = 12.3 B = 45. 

If the output data field is shorter than indicated in the field specification, blanks are 
supplied to fill the field. If the output data field is larger than the field width specified, deci- 
mal digits in the data value will be truncated whenever possible to reduce the data value to the field 
width specified. If the value cannot be truncated; for example, if the data field is 12345 and the 
field specification is 14, program execution will be interrupted and an error message returned. 

Data Records 

A data file is composed of data records. A data record consists of any number of concatenated 
data fields which are terminated by a Carriage Return or a Line Feed. A data record therefore nor- 
mally appears on a physical line. A data record may be continued on the next line by including an 
exclamation point (!) just before the Carriage Return or Line Feed. When an exclamation point is en- 
countered, reading is transferred immediately to the beginning of the next line ignoring the rest of 
the current line;, namely, the Carriage Return or Line Feed. By using the exclamation point, a data 
record of any length can be written. 2 



"If the field is entirely filled, an I field specification may be used only with integer data. For 
a complete list of such restrictions, see the section on tieid specifications. 

On output, a record of any length may be created by suppressing the Carriage Return with a Z in the 
format.. 
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Records might appear as follows: 

Record 1 12345678,910,111.11167249.7 

, ^ 1795.3,7,11, ,65411,3! 
Record 2 213579:6 

Data Record - Format Compatibility 

Normally, every data record is read using a separate FORMAT statement or by rescanning an 
entire FORMAT statenwnt with each data value Cdata field) being read by its corresponding field speci- 
fication. The Carriage Return Cor Line Feed) at; the end of each record corresponds to the Carriage 
Return supplied each time a FORMAT statement is scanned. 

Example 

The data file: 

1234.56789AAC 
7965.3477,8RTF 
9773.996RMT 
will be read by: 

READ (3,9) N1,A1,B1,R1 
9 FORMAT (12, F5.2, F3.1, A3) 
READ (3,9) N2,A2,B2,R2 
READ (3,10) X,Y,Z,L 
10 FORMAT (13, 4X, II, A2, LI) 

as : 

Nl =12 Al = 34.56 Bl = 78.9 Rl = AAC 
N2 = 79 A2 = 65.34 B2 = 77.0 R2 = 8RT 
X = 977 Y=6 Z=RM L=T (true) 

It is possible however, to read or write a number of data records using a single FORMAT state- 
ment by specifying a slash (/) in the FORMAT Statement at the end of the field specifications for each 
record. Remember, a slash (/) generates a Carriage Return which will cause a skip to the beginning of 
the next record. Multiple slashes may be used in the FORMAT Statement to skip records on input or 
generate blank records on output. 

Example 

The data file 12345,6789 
12345678 

will be read by READ (1,7) A,B,I,D,E 

7 FORMAT (F4.2,F5.2,I2/F5.2,J3.4) 



A slash may be used in place of the comma to separate the field specifications. 
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as A = 12.34 B = 5.00 I = 67 D = 12.34 E = 5.678 

It is possible also to read or write on a single record using more than one FORMAT statement 
by including a Z at the end of the FORMAT. The Z suppresses the Carriage Return that is normally- 
issued each time a FORMAT statement is scanned. 

The entire record or file need not be read. For example, the following data file: 
JONES,J.J. 29546793 
SMITH,R.V. 113476995 
will be read by: 

READ (3,10) (A(I),N(I) 1=1, M) 
10 FORMAT (A10,I5) 

as : 

ACl) = JONES, J. J. MCI) = 2954 

AC2) = SMITH,R.V. MC2) = 11347 

If the format indicates that more data values are to be read from a record than actually exist 
on that record, zero values will be supplied at the end to fill the record. 

Example 

The four record data file 
lllllllj 
2222222^ 
3333333^ 
4444444p 
will be read by 

READ (3,100) A,B,C,D 

100 FORMAT (4F5.2) All values read from the same record since no / indicated; 

zeros supplied to fill the field. 

A = 111.11 B = 0.11 C = 0.0 D = 0.0 

Further Examples 

The format statements 

READ (3,100) A,B,C,D 
100 FORMAT (4(F5.2/)) 
and READ (3,100) A,B,C,D 

100 FORMAT (F5.2) 
will both read one value from each record as follows: 

A = 111.11 B = 222.22 C = 333.33 D = 444.44 
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If a slash were included in the second statement above, 

READ C3,100) A,B,C,D 
Ipp FORMAT CF5.2/) 

every other record would be read, A would be A = 111.11, B would be B = 333.33, an end of file would 
be encountered, and the error message UNEXPECTED END OF FILE printed. Double slashes in the first 
statement; for example, FORMAT (4CF5.2//)) would have produced the same result. 

The file 

16,, 5643,, p 
will be read by the format statement 
READ C3,60) A,B,C,D 
60 FORMAT CF5.2) 

as 

A = 0.16 B = 0.0 C = 56.43 D = 0.0 
If a blank field is terminated by a comma, a value of zero is assigned. 
Constructing A Symbolic Data File 

A symbolic data file may be constructed in a number of ways: 

1. It may be read from paper tape. 

2. It may be constructed using the EDITOR commands READ, WRITE and APPEND. 

3. It may be constructed using the EXECUTIVE command COPY. 

4. It may be created by a FORTRAN IV program. 

The data file may contain any number of data records, each terminated by a Carriage Return or 
a Line Feed. Each data record may contain any number of data fields. The data fields may be filled 
completely or terminated by commas. The data fields must be compatible with the field specifications 
used to read them. The data records must be compatible with the FORMAT statements used. Enough data 
must exist on the file to fill all of the variables and/or array elements specified according to the 
formats specified or an End of Data error diagnostic will be given and execution of the program 
halted. 

Terminal Input In Response To A READ Command 

Data input from the terminal in response to a READ command is treated as if it were a symbolic 
data file. When the READ command is executed, the system will ring a bell and wait for the user to 
input one piece of data for each variable or array element in the list. Each data value should then 
be input using one of the three data field forms. The data must follow the field specifications given 
in the FORMAT statement. 
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DISK FILE INPUT/OUTPUT - OPEN. CLOSE 

Disk file input/output may be either symbolic or binary. Symbolic disk file input/output is 
read and written in the formatted form explained in the preceding section. Binary input/output is 
read and written in a format free form explained later in this section. Unlike terminal input/output, 
disk files always must be opened before use and closed after use. 

OPEN And CLOSE Statements 



To read or write on a disk file, the file first must be opened for input or output and its 
form, symbolic or binary, specified. The form of the OPEN statement is as follows: 



r , , INPUT SYMBOLIC 

file ,file 



OPEN " ,/ /, or , or 

[number' name ^^^p^^' ^^^^^^ 



File numbers 2, 3, and 4 are used for disk file input and output. Any file name may be used, 
although it is suggested that the file name be short (four or five characters) . A maximum of three 
files are allowed open at one time. 

After the last input or output operation has been performed with a disk file, the disk file 
should be closed. The form of the CLOSE statement is: 
CLOSE (file number) 

After a file is closed it must be reopened again before it can be reused. Whenever a file is 
closed and reopened, reading or writing starts automatically at the beginning of the file. There is 
no way to continue reading or writing on a file once it has been closed without rereading or rewriting 
to the point at which the user wishes to continue. 

NOTE: File numbers and 1 are reserved for terminal I/O and are never used in an OPEN state- 
ment. Note, however, that when a CONMVNDS FROM file is used, it is closed using the command CLOSECO). 
(See Commands From Files, page 00.) 

Binary/ Disk File Input/Output 

Binary files also may be read and written by a FORTRAN IV program. These files have the dis- 
advantage of not being readable in the EXECUTIVE or in EDITOR since they are written in a binary code, 
but they have the advantage of requiring less storage space than a symbolic file. Binary files must be 



"Another user's data file may be opened for input if 
the file name as follows: 

OPEN l^tf," (Account Number User Name) /©FILE/" , INPUT, 
Example 

OPEN [3 , " (A3R) /@DATA/" , INPUT , SYMBOLIC] 



it has been declared public or contains an @ in 

SYMBOLIC " 
or BINAR'l 
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opened for binary input or output in the OPEN statement. They are read and written using a format- 
free form of the READ and WRITE statements as follows: 



f file 



or ( u / list. 
WRITE V^^^^ 



Example 

OPEN [3, /BIN/, INPUT, BINARY] 
READ C3) A, B, C, (RCI) 1=1,9) 
CLOSE (3) 
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E. DECLARATION STATEMENTS 

It is essential in CCS FORTRAN IV, as in other FORTRAN IV compilers, that specific areas in- 
ternal to the computer be reserved and that certain data relevant to program operation be supplied by 
the program. This information and data is furnished to the compiler in the form of "non-executable" 
statements called Declaration Statements. 

Two of the declaration statements, END and FORMAT, have already been introduced. The FORMAT 
declaration is the only declaration statement that ever appears with a label. Nine declaration state- 
ments will be introduced in this section: Comment, DIMENSION, COMMON, Unlabeled DATA, and the Type 
Declarations (INTEGER, REAL, COMPLEX, LOGICAL, DOUBLE PRECISION). Four other declaration statements 
will be introduced in the following sections: SEGMENT, FUNCTION, SUBPROGRAM, and RETURN. 

COMMENT DECLARATION 

A comment consists of a string of characters commonly used to supply information (documenta- 
tion) on the program. Comments are ignored entirely when the program is executed. A comment may 
appear anywhere in the program except as the last statement (END must be the last statement) . A com- 
ment is indicated by placing either a C: or an * at the beginning of the line as follows: 
C:THIS IS A COMMENT 

*COMMENTS COMMONLY ARE USED FOR PROGRAM DOCUMENTATION 
DATA STATEMENTS 

Occasionally it is desirable to have certain data values reside within the program rather than 
use input or assignment statements. The DATA statement allows the user to store his data in the pro- 
gram. A DATA statement has two general forms as follows: 
DATA Variables/Value^/, Variable2/Value2/ 
DATA Variables,Variable2,Variable2/Values ,Value2 jValue^/ 

For example, both of the following DATA statements would set A = 1.1, B = 2E07, C = T (true), 
and I = 110. 

DATA A/1. 1/,B/2E07/,C/T/, 1/110/ 
DATA A,B,C,I/1.1,2E07,T,110/ 

Only constants and logical values (T and F) may be used as values in a DATA statement. An 
implied DO loop may be used with the second form of the DATA statement to assign values to an array 
as follows: 

DATA (Z(I), 1=1,5)71,2,3,4,5/ 
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If the same data value is to be repeated, it may be specified by placing the number of repeti- 
tions desired, an asterisk (*) , and the data value inside of the slashes as follows. The DATA state- 
ment 

DATA CZCl),I=l,5)/5*2/ 

will assign a value of two to all elements of the array AQ) through ACS). If the number of values 
does not match the niimber of variables, the value list will be rescanned from left to right until all 
of the variables are filled. The DATA statement 
DATA CZCI),I=1,5)/1,2,3/ 

will assign a value of one to ZCl) and ZC4), a value of two to ZC2) and ZC5), and a value of three to 
Z(3). 

Data statements may appear with or without a statement label as: 
DATA A/5/ unlabeled 
2 DATA B/7/ labeled 

Labeled and unlabeled data statements are treated quite differently and it is important that 
the user understand this difference. Unlabeled data statements are used to initialize variables, and 
are executed only once; namely, just before the first executable statement in the program. Labeled 
data statements technically are not declaration statements since they are executable statements. They 
are executed according to their position in the program, they are executed every time that they are 
encountered, and may be used to initialize and reinitialize variable values. 

Example 
>ENTER 10 CIO) 



10. 




1 CONTINUE 


20. 




DATA A/1/ 


30. 




2 DATA B/2/ 


40. 




DISPLAY A,B 


50. 




A=3 


60. 




B=4 


70, 




GG TO 1 


80. 




END 


90. 


CDC) 




>RUN 






1. 


2. 




3. 


2. 




3. 


2. 




3. 


2. 




3. 


2. 




3. 


2. 




3. 


2. 




3. 


2. 




3. 
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DIMENSION DECLARATION 

All arrays used in FORTRAN IV must be dimensioned to reserve sufficient space in the computer 
memory for all elements of the array. Dimension statements identify all array variables in the pro- 
gram by name and indicate how many memory locations are to be assigned to each array variable. The 
general form for a one dimensional array is as follows: 

DIMENSION array name (fii^st element: last element) 

Any legal variable name may be used as an array name. However, if a variable is used as an 
array name, it should not be used as a variable name. The boundaries of the array appear in paren- 
theses. The lower boundary or first element is presented first; the upper boundary or last element 
of the array is presented last. The array elements must be specified by integer constants in the 
main program. In a subprogram they may be specified by variables. The elements may be positive or 
negative integers. For example, the declaration 
DIMENSION AC-4:6) 

would reserve space for the eleven element array AC-4) through A(6) . If the first element of an array 
is one, only the last element need be listed. 

For example, the declaration 
DIMENSION MINC25) 
is equivalent to the declaration 
DIMENSION MINC1:25) 

Multidimensional arrays use an extended version of the DIMENSION declaration as follows: 

P,j„gj^gjQj, array J. first . last first . last 

name ^element ^ 'element element2 * element2* " 

For example, the dimension declaration 
DIMENSION R(0:6,-4:4) 
would set up an array with the range RC0,-4) through RC6,4). 

Space is reserved for multidimensional arrays as though they were linear strings. The ele- 
ments are stored in the string column by column varying the first dimension first, the second second, 
etc. For example, the array ACM,N) would be stored as follows: 

^1,I'^2,1»^3,^*•••^1, 1*^1,2*^2,2 '^3, 2' ••*\,2"--\,N 

An array also may be dimensioned in a Type Declaration statement using implied dimensioning. 
Implied dimensioning; for example, COMPLEX AC10,5), will be covered in the section on Type Declaration. 

If an array is dimensioned in the main program, there is no need to redimension it in a sub- 
program if it is passed to the subprogram as a parameter. If a COMMON declaration is used to pass the 
array however, the array must be redimensioned in the subprogram, 
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COMMON DECLARATION 

COMMON declarations are used to assign data to an area of "common storage" which may be re- 
ferred to by the main program and all subprograms Csubroutines , functions, and segments) in which a 
COMMON declaration appears. Unless a variable is specifically declared to be COMMON it will exist 
(have a particular value) only in the main program or subprogram in which it was given the value. The 
COMMON declaration has the following form: 
COMMON variable list 

Example 

COMMON A,B,I 

The variable list consists of the names of single variables or arrays which are assigned to 
common storage in the order in which they are listed. The COMMON declaration is position oriented; 
that is, the variable name itself is not important (except that it specifies the mode of the data), 
it is the variable position in the list that must be considered. For this reason a common value need 
not have the same variable name in the subprogram that it does in the main program. The variables 
must, however, have the same mode (integer, real, etc.).-'- 

Example 

COMMON A,B,C Main Program 

COMMON X,Y,Z Subprogram 

In the above situation A has the same value as X, B the same as Y, and C the same as Z. 

NOTE: If an array is passed to a subprogram using just the array name, the array must be re- 
dimensioned in the subprogram. 

DIMENSION A(10,10) . 

COMMON A -/ P^°Sram 

COMMON C V 

DIMENSION C(10,10) J Subprogram 

If the array dimensions are included in the COMMON statement, the array is not redimensioned 
in the subprogram. 

DIMENSION A(10,10) 
COMMON A(10,10) 

COMMON C(10,10) Subprogram 



^ Main Program 



One exception exists to this rule. REAL and INTEGER modes may be interchanged. This is possible be- 
cause REAL and INTEGER values both occupy two words of storage. Since the COMMON statement is posi- 
tion (word length) oriented, a variable which is declared REAL in the main program and INTEGER in the 
subprogram will still be reading the correct number of words and therefore maintain the correct posi- 
tion orientation. 
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TYPE DECLARATIONS 

Type Declaration is used to override the implied mode of a variable or function name. Unless 
the mode of the variable or function name has been changed explicitly in a type statement, all vari- 
ables and fimctions beginning with the letters I, J, K, L, M, or N are typed as integers and all other 
variables and function names are typed as real. The general form of the Type Declaration Statement is 
as follows: 

Type of Variables variable list 
The type specified may be INTEGER, REAL, DOUBLE PRECISION, COMPLEX, or LOGICAL. 

Examples 

INTEGER A, POUND, GAL, R 
REAL I,MIN,MAX 

DOUBLE PRECISION ART,HAR,STD,MEAN,X 
COMPLEX B,J,T,COM 
LOGICAL FIR,SEC,TR,N 

The variables in the list following the Type Declaration Statement are all declared to be of 
that mode; thus. A, POUND, GAL, R store integer values; B,J,T,COM store complex values, etc. 

A program may have as many Type Declaration Statements as are required. All of the Type Dec- 
larations are executed before the first "executable" step in the program; therefore, a type declara- 
tion may not be used to change the mode of a value while the program is being run. The mode of a 
value may be changed during program execution by assigning the value to a new variable Ca value assumes 
the mode of the variable in which it is stored), or by using one of the intrinsic functions specific- 
ally designed to convert a variable from one mode to another. 

Dimensioning With Type Declaration Statements 

The Type Declaration Statement also may be used as an implied dimension statement. Variables 
that are arrays may be subscripted in the Type Declaration to indicate the desired dimensions of the 
arrays, thus eliminating the need for a DIMENSION statement. For example, the following DIMENSION and 
Type Declaration statements 

REAL MIN,MAX 

DIMENSION MINC16,16) ,MAX(20,20) 
may be replaced by a single Type statement using implied dimensioning as follows: 

REAL MIN(16,16) ,MAXC20,20) 
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F. SUBPROGRAMS - FUNCTIONS AND SUBROUTINES 

A subprogram consists of a group o£ one or more statements which are stored together outside 
the main program and which may be called by a name assigned to the group. There are two general 
categories of subprograms; the function subprogram and the subroutine subprogram. 

A function subprogram is designed to return a single result and is called merely by the ap- 
pearance of the function name with the appropriate arguments. Three categories of function subpro- 
grams are covered -below: the Statement Function, Library Fianctions, and the Function Subprogram. 

A subroutine subprogram is not designed to return a specific result. A subroutine must be 
called explicitly by a CALL statement and may or may not return values. 

Two separate steps are necessary if a subprogram is to be used: 

1. The function or subroutine must be defined unless it is a library function 
(which is defined internally). 

2. The function or subroutine must be called. 
STATEMENT FUNCTIONS 

Statement functions are used when a single expression is to be performed repeatedly in the 
main program or one of the subprograms . 

Defining A Statement Function 

statement functions are defined with a single statement as follows. 

function r dummy "I ^ expression, 
name ^arguments J 

Example 

FUNCT[A,B] = A**2 + B*2 - C 

Any legal variable name may be used as a function name. The function will be of the same mode 
as this variable. The dummy arguments may be simple variables or array names Cthey may not be con- 
stants, subscripted variables, nor expressions). The dummy arguments are replaced at the time of exe- 
cution by the actual argtiments, which are supplied when the function is called. Every function must 
have at least one dummy argument. The dummy arguments are not considered to have mode. Mode is as- 
sumed only after the actual argument has replaced the dummy argument. 

The expression defines the computations which are to be performed whenever the function is 
called. In addition to the variables appearing as dummy arguments, the expression also may contain 
variables from the program or subprogram in which the function is defined and called. Note however 
that statement functions are not recursive; that is, the function name may never appear in the 
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expression. For example, the expression FUN[A,B] = A+B+FUN[3. 1 ,4. 1] is not allowed. 

Calling A Statement Function 

Statement functions may be called only in the program or subprogram Cfunction, subroutine or 
segment) in which they are defined. The function is called by the appearance of the function name 
and the argxament list (in brackets) . A function may be used anywhere that an arithmetic expression 
may be used, such as in a replacement statement CA = FUNT[3. 1,5.9] + 9*X+2) , or in a DISPLAY statement 
(DISPLAY FUNT[5,6]). The appearance of the function with the actual arguments in brackets causes the 
computations specified in the expression to be performed. The resulting value then replaces the func- 
tion reference in the arithmetic expression, DISPLAY statement, etc. The value returned will assume 
the mode of the function. NOTE: Unless specified in a Type statement, all function names beginning 
with I through N will return integer values, all others will return real values. The actual arguments 
used in calling a function may be constants, subscripted or non-subscripted variables, array names, or 
arithmetic expressions. 

Example 

>ENTER 10(10) 

10 C: MAIN PROGRAM 

20 SROOT [A,B] = (A**2+B**2) ** (1/2) 
30 A = 3. 

40 ANS= SROOT [A, 4] 

50 DISPLAY 'SQUARE ROOT = ' , ANS 

60 END 

70 DC 

>RUN 

The function was defined in statement 20. The function was called in statement 40. The fxmc- 
tion was not given a specific type, therefore the real mode was assumed. 

LIBRARY FUNCTIONS 

Many functions which are commonly used are stored permanently on the system as library func- 
tions. The library includes a variety of functions which may be used for computing trigonometric 
functions (angles must be specified in radians) , comparing variable lists, converting modes of variables 
within an expression, computing logarithmic functions, and so forth. 

The library functions may be used simply by naming the function and placing the argument list 
in brackets as follows: 

function name {arguments] 
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Example 

SIN [2.3145] 

The library functions are treated as expressions and may be used in replacement statements 
(A = SIN [1.596], in DISPLAY statements (DISPLAY EXP [1.957]), etc. The following list includes all 
of the CCS FORTRAN IV functions in the library. The list specifies the form of the function, the 
mathematical expression represented, the number of arguments required, the mode of the arguments, and 
the mode of the function. 



Function 


Statement Form 


Expression 


No. 
Arg's 


Mode 
Argument 


of 

Function 


Arctangent 


ATAN[arg] 
DATAN[arg] 
ATAN2[arg^,arg2] 
DATAN2[arg2,arg2] 


arctan (a) 
arctan CX,Y) 


1 
1 
2 
2 


Real 
Double 
Real 
Double 


Real 
Double 
Real 
Double 


Cosine 


COS[arg] 

DCOS[arg] 

CCOsfarg] 


cos Ca) 




Real 

Double 

Complex 


Real 

Double 

Complex 


iiy L\J\j xj.\^ xii\^ 

Hyperbolic Sine 
Hyperbolic Tangent 


SINH[arg] 
TANH[arg] 


cosh ( a.^ 
sinh Ca) 
tanh Ca) 




Real 
Real 
Real 


Real 
Real 
Real 


Sine 


SIN[arg] 

DSIN[arg] 

CSIN[arg] 


sin Ca) 


I 


Real 

Double 

Complex 


Real 

Double 

Complex 


Exponential 


EXP[arg] 

DEXP[arg] 

CEXP[arg] 


ea 




Real 

Double 

Complex 


Real 

Double 

Complex 


Natural Logarithm 


ALOG[arg] 
DLOG[arg] 
CLOG[arg] 


In Ca) 




Real 

Double 

Complex 


Real 

Double 

Complex 


Conmion Logarithm 


ALOG10[arg] 
DLOG10[arg] 


logioCa) 




Real 
Double 


Real 
Double 


Square Root 


SQRT[arg] 

DSQRT[arg] 

CSQRT[arg] 






Real 

Double 

Complex 


Real 

Double 

Complex 


Absolute Value 


ABS[arg] 
IABS[arg] 
DABS[arg] 
CABsfarg] 


|a| 

/Ca^ + b^) 
for a + bi 




Real 
Integer 
Double 
Complex 


Real 
Integer 
Double 
Real 


Maximum Value 


AMAXO 

AMAXl 

MAXO 

MAXl 

DMAXl 


MaxCaj^,a2. .) 


2+ 
2+ 
2+ 
2+ 
2+ 


Integer 

Real 

Integer 

Real 

Double 


Real 

Real 

Integer 

Integer 

Double 
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Function 


Statement Form 


Expression 


No. 

Arg ' s 


Mode of 
Argument Function 


Minimum Value 


AMINO 

AMINl 

MINO 

MINI 

DMINl 


MinCa^^ ,a2 . . ) 


2+ 
2+ 
2+ 
2 + 
2+ 


Integer Real 
Real Real 
Integer Integer 
Real Integer 
Double Double 


Remaindering 
(Modular arithmetic) 


AM0D[arg^,arg2] 

DMOD[arg3^arg2] 

M0D[argparg2] 


a l^lliUU ^2 J 


2 
2 
2 


Real Real 
Double Double 
Integer Integer 


Positive Difference 


DlM[arg^,arg2] 
IDIM[arg^,arg2] 


aj^-Min (aj^,a2) 


2 
2 


Real Real 
Integer Integer 


Truncation 


AINT 

INT 

IDINT 




1 
i 
1 


Real Real 
Real Integer 
Double Integer 


Float 


FLOAT [arg] 


Conversion 

XIlLCgCX LU IxUai. 


1 


Integer Real 


Fix 


IFIX[arg] 


Conversion 


1 


Real Integer 


Transfer of Sign 


SIGN[arg^,arg2] 

ISIGN[arg^,arg2] 

DSIGN[arg^,arg2] 


Sign a2 
times 1 a^ | 


2 
2 
2 


Real Real 
Integer Integer 
Double Double 


Most Significant Part 
of Double Precision 
Argument 


SNGL[arg] 




1 


Double Real 


Real Part of Complex 

A T* CTl ITTlf^n t" 


REAL [arg] 


a for a + bi 


1 


Complex Real 


Imaginary Part of 
Complex Argument 


AIMAG[arg] 


b for a + bi 


1 


Complex Real 


Single Precision 
Argument in Double 
Precision Form 


DBLE [arg] 




1 


Real Double 


Two Real Arg's in 
Complex Form 


CMPLX[arg^,arg2] 




2 


Real Complex 


Obtain Conjugate of 
a Complex Argument 


CONJG[arg] 




1 


Complex Complex 



FUNCTION SUBPROGRAMS 

Most useful functions cannot be defined by a single expression and therefore must be stored in 
a function subprogram. A function subprogram may contain any set of computations desired, but is de- 
signed to return one specific value, although other values may be returned. The function subprogram 
possesses a certain degree of independence and often may be debugged and run separately from the main 
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program. A function subprogram may call any subprogram Cfunction, subroutine, or segment) including 
itself, and thus may be used recursively. 

Defining A Function Subprogram 

A function subprogram always must have as its first statement a FUNCTION declaration and as 

its last statement an END declaration. Also, it must contain at least one RETURN statement, at least 

one dummy argument, and one value must be assigned to a variable with the same name as the function. 

The FUNCTION declaration tells the system that it is now entering a function subprogram, sup- 
plies the name of the function subprogram, and lists the dummy arguments. The general form of the 
declaration is as follows: 

FUNCTION function name [dirnimy arguments] 

RETURN 

END 

Any variable or array name may be used as a dummy argument. These arguments are merely place 
holders and will be replaced by the actual arguments specified when the function is called. At least 
one dummy argument must be specified, although more may be specified if desired. (Diommy arguments are 
covered more fully at the end of this section.) 

Within the function subprogram, the value that is to be returned must be stored in a variable 
that has the same name as the function. The value stored in this variable will be returned to the 
place at which the function was called when a RETURN statement is executed. 

Example Of A Function Subprogram 
FUNCTION HYP [A,B] 

IF (CA .EQ. 0) .OR. (B .EQ. 0)) GO TO 100 

HYP = SQRT[At2 + B+2] 

RETURN 

100 DISPLAY "ONE OF THE SIDES IS 0" 

RETURN 

END 

Type Specification Of A Function Subprogram 

The mode of the variable used as the function name determines the mode of the value returned. 
The implied mode will be assumed unless a specific type is specified using a Type Declaration State- 
ment in the calling program or in the FUNCTION statement. If the mode of the function subprogram is 
to be specified in the FUNCTION declaration, the following general form is used: 
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Type of Variable FUNCTION function name [argiament list] 
The mode specified may be INTEGER, REAL, COMPLEX, LOGICAL, or DOUBLE PRECISION. 

Example 

REAL FUNCTION INT [A,B,I,R] 

INT = A+B*I/R 

RETURN 

END 

In the above example the value returned, INT, will be in the real mode. Had the type not been 
specified the integer mode would have been assumed. 

Calling A Function Subprogram 

A function subprogram is called automatically whenever its name appears followed by the actual 
arguments required. When the function is called, control will be transferred to the statement follow- 
ing the FUNCTION declaration. The function subprogram will be executed step by step until a RETURN 
statement is executed. Then the value presently stored in the variable with the same name as the 
function will be returned and execution of the main program continued. 

A function call is treated as an arithmetic expression and may appear anywhere that an arith- 
metic expression is legal. The actual arguments may be specified as constants, predefined subscripted 
or non- subs crip ted variables, array names, or arithmetic expressions. 

For example, the function subprogram HYP defined earlier might be called in a replacement 
statement. 

R = 5.1 

HY = HYP [R,3.2] 

NOTE: All variable values are local to the subprogram in which they are defined. If the same variable 
value is to be used in a main and subprogram it must be "shared". Values, other than the "function 
value" which is automatically returned to the main program, may be "shared" by the main and subpro- 
grams in one of two ways. They may be passed either as parameters in the argument list or declared in 
COMMON using a COMMON declaration in both the main and subprogram. 

Example 

The function computes the mean of a series of numbers. The number of values in the series is 
given in the main program. 
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10 C: MAIN PROGRAM 

20 DIMENSION RC200) 

30 ACCEPT R Entire array accepted. 

40 RM = MEAN [R,200] Function subprogram called. 

50 DISPLAY $MEAN = $,RM 

60 END 

70 C: FUNCTION DEFINED 

80 REAL FUNCTION MEAN [S,N] 

90 Z = 

100 DO 100 I = 1,N 

110 Z = Z+SCI) 

120 100 CONTINUE 

130 MEAN = Z/N 

140 RETURN 

150 END 

Following is an example of a recursive fxjnction subprogram (a subprogram which calls itself). 

The function computes the factorial of any number N. 

50 R = FACTORIAL [7] 

60 END 

70 FUNCTION FACTORIAL [N] 

80 IF CN.EQ.O) GO TO 5 

90 FACTORIAL = N*FACTORIAL [N-1] 

100 RETURN 

110 5 FACTORIAL = 1 

120 END 

SUBROUTINE SUBPROGRAMS 

A subroutine subprogram differs from a function subprogram in that a subroutine subprogram 
is not designed to return one specific value to the calling program. A subroutine generally consists 
of a set or group of commands which are repeated frequently in the program. By storing the commands 
in a subroutine, the command set need only be listed once and may be called any number of times from 
anywhere in the program or subprogram. 

Defining A Subroutine Subprogram 

A subroutine subprogram always must have a SUBROUTINE declaration as its first statement, an 
END as the last statement, and at least one RETURN statement. The general form of the SUBROUTINE dec- 
laration is as follows: 
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SUBROUTINE subroutine name [dummy argument list] 



RETURN 
END 

Using arguments in a subroutine subprogram is optional, there need not be any arguments. 
Including arguments is merely a convenient way of referring to variable values in the calling program 
or subprogram. All of the variable values in the subprogram and main program are local variable 
values; that is, they apply only to the subprogram (or main program) in which they were defined and 
therefore do not exist in any other subprogram unless they are transferred as parameters in the argu- 
ment list or declared in COMMON using a COMMON declaration. If diommy arguments are used they must be 
non-subscripted variables or array names. As with dummy arguments in Function Subprograms, no mode 
is assumed unless it is specifically declared using a type statement in the subroutine, in which case 
the mode of the actual argument will be compared with the mode of the dummy argument and an error mes- 
sage returned if they are not the same. 

Example 

SUBROUTINE EXPRESSION [X,Y,Z] 

Z = (X**2+Y**2) 

RETURN 

END 

In this case the value Z was transferred back to the main program as part of the parameter 
(argument) list. It might instead have been shared with the main program using a COMMON statement. 

Calling A Subroutine Subprogram - The CALL Statement 

A subroutine subprogram must be called explicitly with the command CALL. It may not be re- 
ferred to as a function merely by the appearance of its name. The CALL statement has the following 
general form: 

CALL subroutine name [actual argument list] 

When the CALL statement is executed, the dummy arguments are assigned the values of the 
actual arguments and control is transferred to the statement following the SUBROUTINE declaration. 
The subroutine statements then will be executed in the order indicated until a RETURN statement is 
encountered which will transfer control out of the subroutine to the statement following the CALL 
statement. 
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A subroutine subprogram may be called from the main program or from any subprogram. Sub- 
routines are recursive. 

For example, the CALL statement for the previously defined subroutine might appear as 

follows: 

CALL EXPRESSION [T- 1.5, 7. 9, ANSWER] 

The entire program might appear as follows. 

10 C: MAIN PROGRAM 

20 READ CO, 50) T 

25 CALL EXPRESSION [T-1. 5, 7. 9, ANSWER] 

30 WRITE (1,50) ANSWER 

35 50 FORMAT CF7.2) 

40 END 

44 C: SUBROUTINE SUBPROGRAM 

45 SUBROUTINE EXPRESSION [X,Y,Z] 
55 Z = (X**2+Y**2) 

60 RETURN 

65 END 

NOTE: A subroutine subprogram does not return a specific value. Therefore, the subroutine 
name does not have mode. Do not confuse subroutine and function subprograms. 

Example 



5 


C: MAIN PROGRAM 


10 


REAL LOW 


15 


COMMON LOW, HI 


20 


ACCEPT A,B 


25 


CALL ABSOLUTE [A,B] 


30 


SQ = SQRT [HI - LOW] 


35 


END 


40 


C: SUBROUTINE SUBPROGRAM 


45 


SUBROUTINE ABSOLUTE [X,Y] 


50 


COMMON HI, H2 


55 


IF CX.GT.Y) GO TO 10 


60 


HI = X 


65 


H2 = Y 


70 


RETURN 
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75 10 HI = Y 
80 H2 = X 

85 RETURN 
90 END 

ARGUMENTS 

Two types of arguments are foimd in subprograms, the dummy arguments specified in the defi- 
nition of the subprogram, and the actual arguments specified when the subprogram is called. Subprogram 
arguments are a way of passing information Cvariable values, array values etc.) between the calling 
program and the subprogram. 

Dummy argiMents are merely place holders. No space is reserved for them in memory, they 
have no mode, and, in the case of arrays, they are never dimensioned. Only non-subscripted variables 
and array names may be used as dummy arguments. When a subprogram is called, the dummy arguments will 
assiame the value, mode, and dimension of the actual arguments listed in the calling statement accord- 
ing to their order in the list. Because the replacement process is position oriented, the number of 
arguments specified must be the same in both the actual and the dummy argument list. However, the 
name of the variables and/or arrays used to store the values may be changed. 

Although the dummy arguments are not normally considered to have a mode, mode may be assigned 
to them by including a Type declaration statement within the subprogram. When the dummy arguments are 
assigned mode, the mode of the actual argument is checked and must agree with the declared mode of the 
dummy argument or an error message is returned halting program execution. 

The actual arguments which are given when the subprogram is called may be variables, array 
elements, array names iir\ which case the entire array is shared), constants, or arithmetic expressions. 
A constant or arithmetic expression is used only to initialize the dummy variables. In this case 
information may be passed only to the subprogram, any values assigned in the subprogram to the corre- 
sponding dummy variable will not be transferred back to the calling program. If the actual argument 
is specified as a variable, array element, or array name, information may be passed both to and from 
the subprogram. In this case if the value of the corresponding dummy variable is changed in the sub- 
program, the value stored in the actual argument will be changed. 

RETURN 

Every Function and Subroutine subprogram must contain at least one RETURN statement. The 
RETURN statement is a control command which tells the system to leave the subprogram and return to the 
calling program. 
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G. PROGRAM SEGMENTATION 



Total program size may be increased from 200 statements to about 600 statements using program 
segmentation. Program segmentation requires that the total program be broken into a main program and 
segments which may be called from the main program or from another segment. As many as eight segments 
may be called provided that the entire program Call segments combined) does not exceed 600 statements 
and a single segment does not exceed 200 statements. While the program is rvmning, only one segment 
is in memory at one time. The segments not being used are kept on a high speed storage drum. When a 
segment is called it is brought into the computer memory and overlays the segment from which it was 
called. The old segment is stored on the drum and may be recalled if it is needed. 

DEFINING A PROGRAM SEGMENT 

A program segment is defined by a SEGMENT declaration statement as follows: 
SEGMENT segment name 

END 

A segment declaration statement may appear anywhere except in the first or main section of the 
program. Each segment is considered to be a separate program and may contain one or more subprograms. 
The main program and each segment must be terminated by an END statement. NOTE: If the end of a seg- 
ment and the end of a subprogram coincide, only one END statement is required. 

A COMMON statement may be used in the segment to transfer variable values to and from the main 
program and the segments. 

CALLING A PROGRAM SEGMENT 

A segment may be called in two ways as follows: 
CALL segment name 
or 

GO TO segment name 

When a CALL command is used to call a segment, control is returned to the statement following 
the CALL statement after the segment has been executed. 



NOTE: Subprograms may not be shared by the segments. To use the same subroutine in two segments, 
the subprogram must be defined in both segments. The COPY command may be used to copy the range of 
statements desired. 
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When a GO TO command is used to call a segment, control is not returned to the main program 
or segment from which the segment was activated. 

When a GO TO transfer is used, all variable values from the calling segment not declared in 
COMMON will be released when the transfer is made to the new segment to allow more storage of vari- 
ables local to that segment. 

FINISH 

The execution of a segment may be terminated before reaching the END statement by using a 
FINISH statement in the segment. Any number of FINISH statements may be placed in a segment. 
NOTE: A FINISH statement in the main program will halt execution of the entire program. 
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>LIST 



10. COMMON A,B>C>AR*SA 

20. ACCEPT "TYPE SIDES IN DESCENDING ORDER"* A, B, C 

30. AR=AREACA*B3 

40. DISPLAY "AREA OF LARGEST SIDE= "»AR 

50. GO TO VOLWE 

60. END 



70. FUNCTION AREACS1>S2] 

80. AREA = S1*S2 

90. RETURN 

100. END 



110. SEGMENT VOLUME 

120. COMMON A* B*C»AR,SA 

130. IF<AR.EQ.O) FINISH 

140. VOL = AR*C 

150. DISPLAY "VOLUME = "*VOL 

160. CALL SUR 

170. DISPLAY "SURFACE AREA = "*SA 

180. DISPLAY "END" 

190. END 



200. SEGMENT SUR 

210. COMMON A*B*C#AR*SA 

220. SA=2*(A*B+A*C+C*B) 

230. END 

>RUN 

TYPE SIDES IN DESCENDING 0RDER4, 3,2 
AREA OF LARGEST SIDE= 12. 

VOLUME = 24. 



SURFACE AREA = 52. 
END 

STOPPED AT: 180. 
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H. COMMANDS FROM FILES 

A COMMANDS file is a file on which the user may store CCS commands and data which would nor- 
mally be supplied from the terminal. For example, a COMMANDS file might appear as follows: 
QUIT^ 

COPY /DATAl/ TO TELp 
COPY /DATA2/ TO TELp 
CONTINUE p 
CONTINUE^ 

1.1, 3.5, 9.4, 3HENDp 
7.996, 34579, GRB^j 

This file would, when called, return control to the EXECUTIVE, copy the files /DATAl/ and /DATA2/ to 
the terminal, return to FORTRAN IV, and continue execution of the program. 

The FORTRAN IV statement COMMANDS /file name/, which is unique to TYMSHARE, allows the user 
to open a COMMANDS FROM file during execution of his program. When the COMMANDS file is opened the 
user will be in the direct (CCS) command mode. The first line of the COMMANDS FROM file must be a 
direct command or a CCS command. The command CONTINUE will cause the system to resume execution of 
the program immediately following the COMMANDS statement. 

Whenever a COMMANDS file is used to start or continue execution of a program Cusing the com- 
mands RUN or CONTINUE), the input statements ACCEPT and READCO,...) will take their input automatically 
from the COMMANDS file rather than from the terminal. A COMMANDS FROM file automatically assumes a 
file number of zero. 

The FORTRAN statement CLOSE (0) closes a COMMANDS file. After the COMMANDS file has been 
closed, the input again will be, taken from the terminal. NOTE: If the user wishes to use terminal 
input while a COMMANDS FROM file is open, he must. open one of the disk files for terminal use as 
follows: 

OPEN [3, "TELETYPE", INPUT] 
COMMANDS FROM files may be created by using the EDITOR commands READ, WRITE, and APPEND, or with the 
EXECUTIVE command COPY TEL TO /file name/. 

CAUTION: Do not confuse the FORTRAN Language statement COMMANDS /file name/ with the CCS command 
COMMANDS /file name/. 
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SECTION 4 
CCS FORTRAN IV COMMANDS 



CCS commands are not part of the FORTRAN IV language, but are features unique to the Tymshare 
System. The CCS commands make FORTRAN IV programs easy and efficient to create, manipulate, and 
operate on the Tymshare System. 

The user is in the command mode when CCS types a ">". He then may use any of the CCS commands 
described in this section. 

LINE NUMBERS 

Each statement of a CCS FORTRAN IV program has a line number. The line number is ignored by 
the compiler and is independent of any FORTRAN statement label. Line numbers range from .001 to 
999.999, inclusive. All FORTRAN statements are compiled and listed according to the order of their 
line niMbers. A segment of a program can be referred to by its range of line numbers. The form 
either can be 

line niimber 
or 

line niomber : line number 

A. ENTERING A PROGRAM 

Entering Statements From The Terminal 

ENTER 

FORTRAN statements may be entered from the terminal in one of three ways. In all forms of the 
command, the word "ENTER" is optional. 

(1) ENTER With A Line Number 
The command: 

>ENTER line number FORTRAN statement 
allows the user to enter a single FORTRAN statement into his program. The statement must terminate 
with a Carriage Return. A Line Feed may be used to continue a line. If the line number specified in 
the command is already the line number of a statement in the program, the statement typed in replaces 
the original statement. If the line number specified in the command is not already given to a state- 
ment in the program, the statement typed in is inserted into the program according to the line number 
order. 



81 



Example 

For the program: 

10. DIMENSION XCIO), ¥(10) 
20. INTEGER Z 
30. ZCI) = XCI) + Y(I) 
40. END 
The commands : 

>25 DO 10 1=1,10 
>ENTER 35 10 CONTINUE 
>20 INTEGER Z(10) 
would produce the program: 

10. DIMENSION X(10), YClO) 

20. INTEGER ZClO) 

25. DO 10 1=1,10 

30. ZCI) = XCI) + YCI) 

35. 10 CONTINUE 
40. END 

(2) ENTER With A Line Number Range 

A range of line numbers can be entered using ENTER as follows: 
>ENTER line number: line number 
FORTRAN statements 
DC 

This form of the ENTER command allows the user to enter one or more FORTRAN statements into his pro- 
gram. This form of entering statements terminates with a Control D CD*^) ♦ When statement input is 
terminated, any statements already in the program with the same line numbers as the statements just 
entered are replaced by the new statements. The statements typed in will be assigned line numbers 
evenly within the specified range. When this form of the ENTER command is used, CCS prompts an @ 
sign at the beginning of every line. FORTRAN statements typed in may be separated either by a Car- 
riage Return or semicolon. When a semicolon is used, more than one statement can be entered on a 
line. A Line Feed is used to continue a line. 

Examp le 

For the program: 

10. DIMENSION XCIO), YClO) 
20. INTEGER Z 
30. ZCI) = XCI) + YCI) 
40 . END 



The commands : 
>20:30 

§ INTEGER ZCIO) 
@D0 10 1=1,10 

@ZCI) = XCI) + YCI); 10 CONTINUE DC 
produce the program: 

10. DIMENSION XCIO), YClO) 

20. INTEGER Z(10) 

23. DO 10 1=1,10 

26. ZCI) = XCI) + YCI) 

29. 10 CONTINUE 

40 . END 

(3) ENTER With Prompted Line Numbers 

To enter a range of FORTRAN statements, the user may specify a starting line number, an 
increment, and a terminating line number in the ENTER command. The form of the command is: 
ENTER line number Cincrement) line number 

The following example illustrates this command: 

>10 (1)100 CReniember, the word ENTER is optional.) 

>10 A=5p 

>11 B=2*Ap 

>12 X=SQRT[A+B^ 

>13 Y+SQRT[A-B]> 
>14 DC-j ^ 

> 

The computer responds to the command by first deleting lines 10 to 100. Then as statements 
are entered from the keyboard, they are numbered beginning with line 10 and incremented by 1 until 
line 100 is reached. After each statement is entered, the computer prints on the terminal the line 
number to be assigned to the next statement entered. If the user wishes to terminate statement entry 
before line 100 is reached, as shown above after line 13, he merely types DC. 

An alternate form of the ENTER command using line number prompting, includes only a starting 
line number and an increment, as follows: 
>1(1) 

In this form, no final line number is specified. This form of the ENTER command is to be used 
only for initial program entry. When the command is used in this form, statements will be accepted 
until the user types a D^. In the example above, line numbering will begin at 10 and increment by 1 
until statement entry is terminated by a D^. 
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In any o£ the three forms of the ENTER command, a syntactically incorrect FORTRAN statement 
will cause a bell to ring; the statement in error will be typed out immediately with an upward arrow 
pointing to the first character which was not acceptable. The user may then correct the statement. 

Example 

The command 

>10 A=B+*D 
causes the system to type out 
A=B+*D 

T 

and if the user then types 
A=B+C*D 

this statement will be accepted as line 10 in the program. 

In all forms of the ENTER command, the word ENTER is optional. Thus the user can create a CCS 
FORTRAN program merely by typing in each source statement preceded by its intended line number. Lines 
can be changed, or new lines inserted simply by using line numbers which are identical to or fall be- 
tween line numbers previously used in the program. 

b. Entering Statements From A Symbolic File 

The command 

>COPY /file name/ TO line number: line number 
causes the FORTRAN statements from the symbolic file specified by the file name to be entered into the 
program just as with the ENTER command. The file must not have line numbers. Syntactically faulty 
statements are typed out with the upward arrow error indicator and are then discarded. 

c. Entering Statements From Paper Tape 

If the user's program is on paper tape, either of two commands may be used to input from paper 

tape. 

If each statement of the program on paper tape has a line number, use the command 
>LOAD TELETYPE 

The system will type 

EACH LINE MUST HAVE A LINE NUMBER. 
BEGIN INPUT. 

Then turn on the paper tape reader. After the tape is read, type a Control D (DC) . 

On the other hand, if the statements do not have line numbers punched on the paper tape, use 
the command 
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>COPY TELETYPE TO line number: line number 

In this case, the system will type 

NO LINE MAY HAVE A LINE NUMBER. 
BEGIN INPUT. 

Now turn on the paper tape reader. After the tape is read, type a Control D (D'^) . Statements 
are assigned line numbers within the range specified. 

These two commands should be used only with paper tape input because faulty statements are not 
printed out until Control D is pressed, which prevents immediate keyboard correction. 
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B. RENUMBERING A PROGRAM 

RENUMBER 

The laser may change the line numbers in his program using any of the forms of the RENUMBER 
command presented below. 

>RENUMBER^ reassigns line numbers to all statements in the program uniformly in the 
range 10 to 100. 

In the following example, the program was entered using line prompting starting from 1 at 
crements of 2. The RENUMBER command reassigns line ntmibers 10 to 100 to the program. 
>l(2)p 

1. C: PROGRAM TO COMPUTE THE AREA OF A TRIANGLE 
3. C: 

5. ACCEPT"ENTER VALUES OF A,B,C",A,B,C 
7. S=CA+B+C)/2. 

9. AREA=SQRT[S*CS-A)*CS-B)*(S-C)] 
11. WRITE CI, 200) A,B,C, AREA 
13. 200 FORMAT C4F16. 8) 
15. STOP 
17. END 

> RENUMBER p 



>LISTp 

10. C: PROGRAM TO COMPUTE THE AREA OF A TRIANGLE 
20. C: 

30. ACCEPT"ENTER VALUES OF A,B,C",A,B,C 

4P. S=CA+B+C)/2. 

50 . AREA=SQRT [S* (S-A) * (S-B) * CS-C) ] 

60. WRITEC1,200) A,B,C, AREA 

70. 200 FORMAT C4F 16. 8) 
80. STOP 
90. END 

> 

>RENUMBER line range^ 
reassigns line numberis uniformly to statements in the line range indicated. 
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Example 

>RENUMBER 1.5:47 

>RENUMBER Old Line Range AS New Line Range ^ 
reassigns line numbers to statements in the old line range imiformly within the new line range. 

Example 

>RENUMBER 1 : 20. 3 AS 10 : 50 

The new line number range must not include any line number belonging to a statement which was 
not included in the old line number range. This rule is to prevent destroying the rest of the program. 

This means that the RENUMBER command cannot be used to erase program lines nor to change the 
order of statements. For example, in the following: 
>1 

>1.5 — - 
>2.3 - — 

>7 

>RENUMBER 1:2.3 AS 1:6 
is correct, but 

>RENUMBER 1:2.3 AS 1:10 

is not allowed since this would erase statement 7. When the same line number is assigned to two state- 
ments by a RENUMBER command, the error message LINE RANGE MISMATCH will be printed. 
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C. LISTING A PROGRAM 

All or part of a program can be listed in formatted form for easy reading, or in unformatted 
form for quicker printing.' 

LIST 

The command 
>LISTp 

types the entire program in a formatted form with line numbers, statement labels, and statements 
aligned vertically. 

Example 

>LISTp 

1. C: PROGRAM TO COMPUTE THE AREA OF A TRIANGLE 

3. C: 

5. ACCEPT"ENTER VALUES OF A,B,C",A,B,C 

7. S=CA+B+C)/2. 

9. AREA=SQRT[S*CS-A)*CS-B)*CS-C)] 
11. WRITEC1,200) A,B,C, AREA 

13. 200 FORMAT C4F 16. 8) 

15. STOP 
17. END 

Example 

>LIST 11:15^ 
11. WRITEC1,200) A,B,C, AREA 

13. 200 FORMAT C4F 16. 8) 

15. STOP 

> 

The command 

>LIST line number ^ 

types out the statement with the specified line number Cif one exists) together with its line number. 

The command 

>LIST line number: line number^ 
lists the part of the program specified by the line number range (if any exists) together with the 
line numbers. 

The user may address any line or range of lines to be listed in a number of ways. He may in- 
dicate the position of a line relative to another line. 
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For example, in the following 
10. A=2.7 
15. B=7.1 
20. ACCEPT C 
25. D=A*C+B 
30. DISPLAY D 
35. END 

the statement ACCEPT C could be specified as line 20, or as line 25-1, or as line 10+2. The latter 
two methods are exanqjles of relative addressing. 

If, for example, the user decides to list lines 15 to 30 in the above program, any of the fol- 
lowing forms could be used: 
>LIST 15:30p 
>LIST 10+l:10+4p 
>LIST 10+l:35-l;p 
>LIST 15:35-lp 

The command 
>LIST. J) 

lists the current line. The period following LIST indicates the current line. 

Relative addressing can be used in all commands that allow line addressing. 

FAST 

An unformatted quick listing can be obtained with the command 
>FAST line range p 

The line range is optional; if it is omitted, the entire program is listed. The following is a list- 
ing of the previous example with the command FAST. 

>FAST^ 

1. C: PROGRAM TO COMPUTE THE AREA OF A TRIANGLE 

3. C: 

5. ACCEPT"ENTER VALUES OF A,B,C",A,B,C 

7. S=(A+B+C)/2 

9. AREA=SQRT[S*(S-A)*(S-B)*(S-C)] 

11. WRITEC1,2Q0) A,B,C, AREA 

13. 200 FORMAT C4F 16. 8) 

15. STOP 

17. END 

> 

All of the line addressing methods described under LIST can be used with FAST. 
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D. DELETING A PROGRAM OR STATEMENTS FROM A PROGRAM 

To delete either all or part of a program, the DELETE and CLEAR commands are used. 

DELETE 

This command will delete parts of a program. 
>DELETE line number ^ 
deletes the line specified. Thus, >DELETE 35 deletes line 35 from the program. 

>DELETE line number: line number^ 
deletes the range of lines specified. Thus, >DELETE 20:40 deletes lines 20 to 40 inclusively. 

Relative line addressing can be used with the DELETE command. For example, 
>DELETE . Deletes the current line. 

>DELETE .-3:100 Deletes from three lines above the current line up to and including 
line 100. 

>DELETE 10+1:10+4 Deletes from the first line after line 10 up to and including the 
fourth line after line 10. 

CLEAR 

This command will erase the entire program. The form of the command is 
>CLEAR ^ 
The system will reply with 
ERASE PROGRAM? 

to which the user will answer by typing Y for YES or N for NO. If the answer is Y, the program will 
be erased and the message OK. typed. If the answer is N, the command will be aborted. In either 
case, the user receives a command prompt (>) . 
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E. EXECUTING A PROGRAM 

To begin execution o£ a program after it is entered, use the conunand 
>EXECUTE^ 

or 

>RUNp 

Program execution will start with the first executable statement. 

When an input statement is executed, a bell will ring. This means that data is to be typed in. 

If an error occurs during execution, the first statement which is in error will be typed with 
an error diagnostic and execution terminated. The user should correct the error and re-run the pro- 
gram. 
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F. STORING AND RETRIEVING A PROGRAM 

SAVE 

After a program has been created, the user may wish to save it on the disk file so he may use 
it again. The command 

>SAVE /file name/ 

writes the program with the line numbers on the file specified. The file name is any name that the 
user chooses to identify his program. Once a program is saved on a file, the user may retrieve it in 
FORTRAN for execution or modification; or he may look at it in EDITOR (EDITOR is part of the Tymshare 
system which manipulates text) . 

If the user wishes to save a copy of his program on paper tape, he may type 
>SAVE 

then turn on the paper tape reader and type a Carriage Return. The program is punched on paper tape 
and listed on the terminal. 

LOAD 

Once a program has been saved on a file the user may retrieve it any time. To retrieve the 
program use the command 

>LOAD /file name/ 

An ALT MODE at this point will abort loading (and hence, preserve the current program). A 
Carriage Return causes the file to be loaded. If a Carriage Return is pressed, the system will type 
O.K. and start loading the program specified. When loading is completed, a > is typed. 

During loading of a program, faulty statements are printed with an up arrow and discarded. 
These statements should be corrected before the program is executed. 

After a program is loaded, it is ready for execution or modification. If the user has modi- 
fied his program and wishes to keep the modifications, he should save the program again either under 
the same file name or under a different file name. 

An important use of the LOAD command is to merge files. For example, if the user is creating 
a program on the terminal, the command >LOAD /A/ will merge by line number the contents of file /A/ 
with the current program. If the same line number occurs in both the current program and in file /A/, 
the line in file /A/ will be the one retained. If the line in file /A/ is syntactically incorrect, 
the line in the current program will be retained. 

Files can be merged by the LOAD command as follows. The file /PROG/ contains the following 

program: 
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10. ACCEPT X,Y 
20. Z=SQRT[X**2+Y] 

30. 'II»iXA*ATAN2[X,Y] 
40. END 

And the file /CORR/ contains the following: 

35. DISPLAY X,Y,Z 

20. Z=SQRT[X**2+Y**2] 
When both files are loaded by executing two LOAD commands as follows, 

>LOAD /PROG/ 

O.K. 

>LOAD /CORR/ 
O.K. 

> 

the following program results: 
10. ACCEPT X,Y 
20. Z=SQRT[X**2+Y**2] 
30. THETA=ATAN2 [X,Y] 
35. DISPLAY X,Y,Z 
40. END 

The two files have been merged by line numbers. Note that the individual files need not be ordered 
by line numbers as is the case with the file /CORR/. 

Since files are merged by the LOAD command, the user can prepare insertions and corrections to 
a program off line on paper tape and merge the changes with the program to be altered. Thus, the cost 
of program development is reduced by reducing terminal time. 

COPY 

This command allows the user to copy a part of a program. The general form of the command is 
>COPY source TO destination 
where source and destination can be either files or line ranges or the terminal. The source is always 
retained. 

Exan^les 

>COPY 10:25 TO 75:100^ Deletes lines 75 to 100 and then inserts the contents of lines 

10 to 25 uniformly renumbered in the range 75 to 100. Lines 10 
to 75 are undisturbed. 
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>COPY 1:100 TO /KZZP/ 



>COPY /KZZP/ TO 200:300;> 



>COPY /FILE A/ TO /FILE B/- 



>COPY TELETYPE TO 40:50 



Copies lines 1 to 100 to file /KZZP/ without line numbers. 
Lines 1 to 100 are undisturbed. 

Deletes lines 200 to 300, then inserts the contents of the file 
/KZZP/ uniformly numbered in the range 200 to 300. The file 
/KZZP/ must not have line numbers. 

Copies the contents of /FILE A/ to /FILE B/. The file /FILE A/ 
may be any type, a symbolic file with line numbers, a binary 
file without line numbers, or a DUMP file. 

Deletes lines 40 to 50 and inserts the contents of what is typed 
on the terminal after this command uniformly numbered in the 
range 40 to 50. Statements typed on the terminal must not have 
line numbers. 



MOVE 

The general form of the command is 
>MOVE source TO destination 
where source is a line range and destination is either a line range or a file name. The MOVE command 
is similar to the COPY command except that the source is deleted after it is moved. 



Examples 



>MOVE 10:25 TO 75:100 ^ 



>MOVE 1:100 TO /KZZP/p 



Deletes lines 75 to 100, inserts the contents of lines 10 to 25 
uniformly numbered in the range 75 to 100, then deletes lines 
10 to 25. 

Moves lines 1 to 100 to the file /KZZP/ without line numbers, 
then deletes lines 1 to 100. 



THE EXECUTIVE COMMANDS DUMP AND RUN 

Two commands in the EXECUTIVE are described here briefly since they may be used to save and 
retrieve FORTRAN program. 



QUIT 



The user may return to the EXECUTIVE from FORTRAN with the command 
>QUIT^ 



The dash typed by the con^uter means that the user is returned to the EXECUTIVE, 
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DUMP 

The EXECUTIVE command 

-DUMP /file name/^ 
NEW FILE j> 
or OLD FILE^ 

saves the entire contents of the user's program which includes the text of the program, the compiled 
program, and all data values on a DUMP file. He may reload the file at a later time to continue where 
he left off. 

RUN 

The EXECUTIVE command 

-RUN /file name/ p 

reloads the FORTRAN program that is saved under a DUMP command. The user is returned to FORTRAN and 
put in the command prompt mode which allows him to do whatever he wishes with his program. 

The use of the EXECUTIVE commands DUMP and RUN is the fastest way to save and retrieve FORTRAN 
programs. It is recommended that the user with large programs use these commands to save time in do- 
ing day-to-day debugging. NOTE: If the user wants to look at the contents of the file in EDITOR he 
must save it with the FORTRAN command SAVE (which saves the symbolic file). 

SHRINK 

SHRINK is used to avoid excessively large DUMP files. SHRINK eliminates unnecessary informa- 
tion such as segment pages, array storage, and all data values to minimize the amount of disk space 
used by the program under the EXECUTIVE DUMP command. The sequence of commands used after the user 
has created a program and wishes to save it as a DUMP file is: 

>SHRINK ^ 

>QUIT ^ 

-DUMP /PROGl/^ 
NEW FILE p 

intervening work 
or LOGOUT 

-RUN /PROGI/4J 
FORTRAN 

> 
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LOCK 

The LOCK command inhibits all commands but 
>RUN or >EXECUTE, >CONTINUE, >QUIT; 
and may be used when a program is to be made available as a running program but not to be listed or 
saved by another user. 

The form is 
>LOCK^ 
The computer types 

IS PROGRAM ALREADY SAVED? 

to which the user will answer by typing Y for YES or N for NO. If he types Y, the computer will type 
O.K. and the program is locked. If he types N, the command will be aborted. If other characters are 
typed, computer will type WHAT? to which he should type Y or N. A command prompt will then be given. 

If the user wants to keep a copy of the program he should save it before the LOCK command is 
used since he will not be able to use the SAVE coimnand after LOCK. 

If a locked program is also declared proprietary with the EXECUTIVE command DECLARE, the 
user is assured complete protection of his program. No one will be able to use it or list it. 

COMMANDS 

The user may have a set of CCS commands saved on a file. When he uses the command 
>COMMANDS /file name/-) 
all the commands in the file will be executed. 

For example, the file /COM/ has the following instructions: 
LIST 
RUN 
QUIT 
FILES 
CONTINUE 
When the user types 

>COMMANDS /COM/^ 

his current program will be listed and executed. QUIT returns him to the EXECUTIVE. The files he has 
in his directory will be listed by the EXECUTIVE FILES command. Then CONTINUE puts him back into 
FORTRAN IV. 
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If this is a standard procedure that he wants to go through with each program, all he has to 
do is to give the commands 

>LOAD /PROGRAM/ p 
O.K. 

> COMMANDS /COM/^ 
each time he enters the system or each time he runs a program. 



97 



G. PROGRAM CONTROL AND DEBUGGING 

For debugging and modifying programs, CCS provides the user with the commands REFERENCES and 
DEFINITIONS. With these commands, the user can refer to the parts of the program he is interested in. 

^ To Locate Variables And Label References And Definitions 

REFERENCES 

The command 

>REFERENCES Label or Variable Line Range ^ 
will type out all statements together with their line numbers which make reference to the specified 
variable name or statement label. Declarative occurrences of the identifier are ignored (such as 
DIMENSION statements. Type Declaration statement, subroutine declaration, etc.). The line range is 
optional. If a line range is used, all non- declarative statements within the line range will be typed 
out. If the. line range is omitted, all non-declarative statements containing the label or variable 
specified in the entire program will be typed out. 

For example, >REFERENCES 30 will cause the statement labelled 30 to be typed out, and all 
statements in the entire program with GO TO 30 to be typed out. 

>REFERENCES Z 10:50 will cause all statements with Z Cexcept declaration statements) within 
the line range 10 to 50 to be typed out. 

DEFINITIONS 

The command 

>DEFINITIONS Label or Variable Line Range^ 
will type out all statements together with their line numbers which define the specified variable or 
statement label. Line range is optional. If it is omitted, all declarative statements containing the 
label or variable specified are typed out. 

For example, 

>DEFINiTIONS Z will print all declarative references to the variable Z. 
>DEFINITIONS Z 10:50 will print declarative references to the variable Z within the line 
range 10 to 50. 

Following are some examples of how the user can make use of the commands REFERENCES and DEFI- 
NITIONS for program debugging and modification. 

1. Suppose during the execution of a program, the user gets the error message CONFLICTING 
TYPE DECLARATION with the statement 3.7 REAL X,Y,Z printed. He may type the command 
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>DEFINITIONS Z ^ 
which will type out the following statements: 
1.4 SUBROUTINE S[L,M,Z] 
3.2 INTEGER T,L,Z 
3.7 REAL X,Y,Z 

This tells him immediately that he had declared Z to be both INTEGER and REAL. He then could modify 
either statement 3.2 or 3.7. Suppose he has decided that Z should be real and has kept the statement 
3.7 REAL X,Y,Z. In running the program again, he gets the message CONFLICTING TYPE DECLARATION with 
the statement 3.7 REAL X,Y,Z printed. He may type the command 

>REFERENCES S^ 
which will type out the following statements: 
9.1 CALL S [5,6,3.9] 

12.2 CALL S [5,7,8] 

14.5 CALL S [11,2,SQRT[17.4]] 
The user immediately finds out that in statement 12.2, he had used the integer 8 in the subroutine 
CALL when it should have been a real number. He then may change 8 to 8. and re-run the program. 

2. Suppose the user has a subroutine in his program called GEM and wishes to add a parameter 
to all references to the subroutine. He may do this simply by using the command 

>REFERENCES GEM 9 

All subroutine calls to GEM will be printed. He then may modify these lines. 

3. Suppose the user encounters the error message SUBSCRIPT OUT OF BOUNDS and the statement 
30. X=AC35) printed. He may type 

>DEFINITIONS A 

which types out the statement 10. DIMENSION AC20) . He then may modify the statement to reserve more 
space for the array A. 

b^; P rog ram Interruption 

ALT MODE/ESCAPE 

A program may be interrupted at any point during execution by pressing the ALT MODE/ESCAPE 
key. The computer will type the message INTERRUPTED BEFORE STEP , and return the user to the com- 
mand prompt mode with a > sign. Execution of the program can be resumed from the point of interrup- 
tion by typing the command 
>CONTINUE 
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Breakpoints 

BREAK 

As many as ten breakpoints can be set easily in the program with the BREAK command as follows: 
BREAK 5 8 10:14 p 

This command sets a breakpoint at lines 5 and 8 and one breakpoint at every line between 10 and 14 in- 
clusive. NOTE: Type a space between each breakpoint desired. During execution, when a line is en- 
countered where a breakpoint has been set, the following occurs: 

1. Execution is interrupted just before executing the statement at the line where the 
breakpoint was set. 

2. The niomber of the line where the breakpoint was set is printed on the terminal 
immediately followed by > on the same line, e.g., 10. 7>. 

3. The user is now in the statement prompt mode Cas indicated by a line number and a >) . 
Execution can be resumed from the breakpoint by typing CONTINUE^ . To list all active breakpoints, 
type BREAK^ . All active breakpoints can be cleared with the >RESETp command. 

dj. Immediate Execution 

Once program execution has been interrupted by the ALT MODE key or by encountering a break- 
point, statements can be entered for immediate execution. 

To enter a statement, the user merely types @ followed by the statement to be executed, fol- 
lowed by a Carriage Return. When the Carriage Return is depressed, the statement just typed will be 
compiled immediately, and, if the statement was syntactically correct, executed. In any case, the 
statement then will be discarded. Any FORTRAN IV statement except a declarative statement may be 
executed in the immediate mode. CCS commands such as LIST and DELETE may be executed also. Program 
execution now may be resumed from the point of interruption by typing CONTINUE^p . 

e. Step Execution 

After program execution has been interrupted, the user may single step through his program, 
executing one statement at a time. This is accomplished by typing NEXTp for each statement to be 
executed. 

The following program demonstrates some of the debugging features. 
LISTi) 



100. Y=.5 

110. Z=.2 

120. LAST=5 

130. 10 DO 15 C0UNT=1,LAST 

140. Y=Y*SIN[Y]+Y 

150. Z=Z*SIN[Z]+Z 
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160. X=X*SIN[X]+X 

170. 15 CONTINUE 

180. A=SQRT[Z**2+Y**2+X**2] 

190. IFCA.GT.l) GO TO 30 

200. 20 LAST=LAST*2 

210. GO TO 10 

220. 30 DISPLAY LAST, A 

230. END 



>BREAK 180 200^ 
> EXECUTE <> 

VARIABLE HAS NOT BEEN ASSIGNED A VALUE 

160. X=X*SIN[X]+X 

160. >@X=5^ 

160. >CONTINUEi? 

180. xaDISPLAY X,Y,Zp 

0.557564177 0.9253519603 0.7916238 
180. >CONTINUE«> 

5 1.339336499 

> 

In this program the user has set two breakpoints, one at Statement 200 and another at State- 
ment 180. He then began execution and encountered an error. The direct statement @X=5 placed the 
value 5 in the variable X. Execution was resumed at Statement 160 by the command CONTINUE. A break- 
point was encountered at Statement 180 and control returned to the user. The user now asked for 
X,Y,Z to be displayed by the direct statement SDISPLAY X,Y,Z. Execution was resumed again by the com- 
mand CONTINUE. The program then typed the answers 5 and 1.339336499 and stopped. 

f^^ Easy Program Checking 

Two commands allow the user to verify quickly the executability of his program without actually 
running it. 

CHECK 

This command tests for the executability of a program by searching for structural errors such 
as no END statement in a subroutine or in the main program, labeling or declaration errors, calls to 
undefined routines, etc. Any error discovered will be printed on the terminal. 

Example 

>CHECK^ 
100 A=3 

DUPLICATE LABEL 
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If there are no structural errors, the following will be printed: 
>CHECKp 
O.K. 

> 

The user now may execute his program by typing RUN. 
INITIALIZE 

This command, in addition to performing all the functions of the CHECK command, allocates data 
storage and then begins execution of the program, "breaking" at the first executable statement. The 
statement with a ">" will be indicated. At this point, the user may enter FORTRAN statements for 
immediate execution or resume program execution by typing CONTINUE 

Example 

>INITIALIZEi> 
10. > 
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H. - EDITING 

One of the most powerful features of CCS FORTRAN IV is the ability to edit while in FORTRAN. 
Editing features allow editing of both program statements and CCS commands. All of the editing char- 
acters available in the text editing language EDITOR can be used in FORTRAN. 



The following is a summary of the editing control characters: 



EDITING CONTROL CHARACTERS 


Control 
Character 


Symbol 
Printed 


Action 


A or 


-<- 


For Deleting 

Deletes the previous character typed. Repeated use deletes 
several characters. 






Deletes the entire line being typed. 




\ 


Deletes the preceding word in the line. 


p'' and a 
character 


% 


Deletes up to but not including the character typed after it. 


„c , 
X and a 

character 


% 


Deletes up to and including the character typed after it. 


c 

K 




Deletes the next character in the old line, prints the character 
it deletes. 


b 


& 
"» 


Deletes the next character in the old line. 


M or ^ 




is the same as a Carriage Return; deletes the rest of the old 
lines and ends the edit. 






For Copying 

Copies and prints the rest of the old line and ends the edit. 






Copies without printing the rest of the old line and ends the edit. 






Copies and prints the rest of the old line; edit continues at the 
end of the new line. 






Copies without printing the rest of the old line; edit continues 
with the new line acting as old line. 






Prints the rest of the old line plus the new line; edit continues. 


T^ 




Same as R^ except that it aligns old and new lines. 


C^ 




Copies the next character in the old line. 


0*^ and a 
character 




Copies up to but not including the character typed after it. 


Z and a 
character 




Copies up to and including the character typed after it. 
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EDITING CONTROL CHARACTERS (Continued) 


UUl 1 UT U 1 

Character 


oyiriDO 1 
Printed 


Acti on 


u 




Copies characters from the old line up to the next tab stop in 
the new line. 




< 
> 


For Inserting 

Inserts text into the old line; first E^ prints <, second E^ 
prints >. 






Others 

Types spaces up to the next tab stop. 


,c 

J or -1 




J is the same as Line Feed; new line continues. 


— c ■ 

N 




Backs up in old and in new line. 



cu Commands Or Statement Input Editing 

1, During normal input of CCS commands or FORTRAN statements, A^, w'^, and are available. 
A Control A or back arrow deletes the previous character typed. 

W Deletes the preceding word in the line. 
Q Deletes the entire line being typed. 

20 IF (X.GRA^-e-T.Y) THEN W^\gO TO 100 
>LIST 

20 IF CX.GT.Y) GO TO 100 

2. If the user types a Control Y (Y^) at the end of a command or a statement, what he just 
typed in becomes the old line for editing purposes. All editing characters now can be used. 



Example 



Example 1 

>LIFT 1:10 y'^^ 
The LIST command is executed. 

z'^ILIST 1:10 ^ 

Example 2 

>15 A=SQRT[X+Y]Y^^ 

c c 
>Z +15 A=SQRT[X+Z]^ He types a Z + which copies the old line up to +, then types Z and ] 



104 



Editing When Syntax Errors Occur 



If a syntax error occurs during the ENTER command, the erroneous line with an up arrow will 
be printed. The line in error is the old line for editing purposes. 

Example 

>10 A=3-*SQRT[3. 14159]^ 
WRONG NUMBER OF OPERANDS: : OPERATORS 
A=3-*SQRT[3. 14159] 
Z^-A=3-E^<Ae'^>d'^*SQRT[3. 14159] 



>LIST 10 
10. A=3-A*SQRT[3. 14159] 

> 

If a syntax error occurs during LOAD or COPY /file/ TO line range, the erroneous line is 
printed, but input continues to be accepted from the file. Thus, bad lines are lost and cannot be 
edited from the terminal. They must be retyped. 

c. Addressing The Lines To Be Edited 

To address a line to be edited, type 
>EDIT line number ^ 

To address a range of lines to be edited, type 
>EDIT line niamber:line number^ 

The EDIT command types out the addressed line Cthe old line) and allows editing to be done on 
that line. EDIT also allows editing of a range of lines when a line range is given. Each line is 
typed out sequentially and made available for editing. When the last line in the line range has been 
edited, control is returned to the user. 

Example 

In the following program, 

10. ACCEPT A,B,C 

20. S=A-B/2 

30 . AREA=SQRT [S* (S+A) *S* (S+B) *S* (S-C) ] 

40. WRITE(1,200) A,B,C,AREA 

50. 200 FORMAT C4F 10. 2) 

60. STOP 

70. END 



The user typed a Z followed by which copied 
the line up to He then typed an E^ which 

printed a <; he inserted the variable A, then typed 
another EC to end the insertion. To copy the rest 
of the line, he typed a D^. 

The corrected line is listed. 
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the user edits lines 20 through 50 as follows: 

>EDIT 20:50 p 
S=A-B/2 £> 
S=CA+B+C)/2^ 

AREA=SQRT [S* CS+A) *S* CS+B) *S* (S-C) ] ^ 
AREA=SQRT [S* (S-A) * CS-B) *?<^{%%%%T)C {S-Q ]p 



WRITE(1,200) A,B,C,AREA^ 
DC ^ 

200 FORMAT (4F 10. 2) 

Z'^l 200 FORMAT (4F16. 8)^ 



The new program is listed as below: 
>LIST 



10. ACCEPT A,B,C 

20. S=CA+B+C)/2 

30. AREA=SQRT[S*CS-A)*CS-B)*CS-C)] 

40. WRITE CI, 200) A,B,C,AREA 
50. 200 FORMAT C4F16. 8) 

60 . STOP 

70 . END 



> 

The command 

>MODIFY line number^^, 

or 

>MODIFY line number: line number^ 
is the same as the EDIT command except that the lines to be edited will not be typed out. 

d^^ Editing Input Data 

Data typed into a running program may be edited with the following editing characters: 
A^ Deletes the preceding character and allows it to be retyped. 

Deletes the current field and allows it to be retyped. 
q'^ Restarts the entire input statement from the beginning. 



Computer prints line 20 and a Carriage Return. 

User types in new line 20 and a Carriage Return 

Computer prints line 30 and a ^ . 

User types new line up to the * after (S-B) , then 

a P'^ followed by a left parenthesis, which deletes 

the old line up to but not including the (. He 

then types in D'' to copy the rest of the line. 

Computer prints the next line. 

User decides that this statement is correct and 

types a D^. 

Computer prints the next statement. 
User types Z^l which copies the old statement up to 
the character 1; he then types in the rest of the 
line and a ^ . 
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Example 

>10(10) 

10. DIMENSION ACS) 
20. ACCEPT A 
30. DISPLAY A 
40. END 



>RUN 

2354AC^,34.6AC^78,23.4,56.7WC\46.8,10.1 
235. 34.78 23.4 46.8 10. 1 

STOPPED AT: 30. 
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SECTION 5 
USER AIDS 

A. ABBREVIATED COMMANDS 

To save time, any command may be abbreviated to as few characters as is necessary to identify 
the command uniquely. Thus the command FAST may be abbreviated simply as 

since there are no other commands that start with the letter F. But to list a program, at least 
>LI^ 

is required since the commands LOAD and LOCK also start with the letter L. 

If insufficient identification is given to a command, the message INSUFFICIENT COMMAND 
IDENTIFIER will be given. 

B. COMMAND MODELS 

If the user is not certain how to use a particular command, he may type a question mark while 
entering the command. A model of the command is then typed out. 

For example, 
>COPY?53 

COPY <LINE RANGE OR FILE NAME> TO <LINE RANGE OR FILE NAME> 
> 

>DEL? ^ 

DELETE <LINE RANGE> 



To obtain a listing of command models for all the CCS commands type a question mark as follows: 
>?<> 

INITIALIZE CPROGRAM) 

EXECUTE 

RUN 

CONTINUE (EXECUTION) 

CLEAR (PROGRAM) 

MAP (STORAGE ALLOCATION) 

NEXT (STATEMENT) 

QUIT (THIS SUBSYSTEM) 

CHECK (PROGRAM STRUCTURE) 

@ <IMMEDIATE FORTRAN STATEMENT> 

SAVE [ <FILE NAME> ] 

LOAD [<FILE NAME>] 

[ENTER] <LINE NUMBER> <SOURCE STATEMENT> <C.R.> 

[ENTER] <LINE>(<INCREMENT>) [<LINE>] <STATEMENTS> <CONTROL D> 

[ENTER] <LINE>:<LINE> <STATEMENTS> <CONTROL D> 

COMMANDS (FROM) <FILE NAME> 
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EDIT <LINE RANGE> 
MODIFY <LINE RANGE> 
DELETE <LINE RANGE> 
LIST [<LINE RANGE>] 
FAST [<LINE RANGE>] 

COPY <LINE RANGE OR FILE NAME> TO <LINE RANGE OR FILE NAME> 

MOVE <LINE RANGE> TQ <LINE RANGE OR FILE NAME> 

RENUMBER [ <OLD LINE RANGE> [ AS <NEW LINE RANGE> ] ] 

REFERENCES <IDENTIFIER> [<LINE RANGE>] 

DEFINITIONS [ <IDENTIFIER> [<LINE RANGE>] 

BREAK <MAX OF 10 LINE RANGES> 

RESET (BREAKPOINTS) 

LOCK 

SHRINK 



CAUTION: A space is required between a command and the rest of the command. For example, 
LOAD/FILE/ will cause an error diagnostic. The user must type LOAD /FILE/. 
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C. A REVIEW OF CCS PROMPTS 

One of the following may be prompted by CCS. 

1. A > is a command prompt. When CCS gives a >, it is ready to accept a CCS command. 

2. When a line number is printed followed by a >, it means execution of program is continu- 
able at that line. For example, when CCS prints 

10. > 

the user may enter any CCS command including an @ followed by a statement for direct execution, or he 
may type CONTINUE to resume execution of the program at line 10. 

3. When CCS prints an e sign, it is ready to accept a FORTRAN statement. The statement will 
be included in the line range indicated. This occurs when the user uses the ENTER command with a line 
number range. 

4. When CCS prints a line number, it is ready to accept a FORTRAN statement. This occurs 
when the user uses the ENTER command in a line prompt mode. 
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SECTION 6 
SAMPLE PROBLEMS 

This section contains programs written in FORTRAN IV and executed on the Tymshare system. The 
problems are presented in an increasing degree of complexity. 
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MONTHLY PAYMENT PROGRAM 



1. DEFINE THE PROBLEM 

A. Input 

1. Original debt (P) 

2. Annual interest CI) 

3. Nxamber of monthly payments to be made CN) 

B. Compute 

Monthly payment (M) 

M = P-UI^l)^ 

N 

Ci+1) -1 

C. Output 

Monthly payment, M. 
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2. 



FLOWCHART 



C START ) 



INPUT: ORIG. 
iDEfeT CP) 
ANNUAL INTEREST (1) 
NO. OF MONTHS (N) 



COMPUTE : 
MONTHLY 
INTEREST 
I = 1/12 



COMPUTE : 
Q = I + 1 
M = P*I*QfN/CCQ^N)-l) 



I 



OUTPUT: 
MONTHLY 
PAYMENT M 

I 

© 
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3. FORTRAN CODE AND SAMPLE EXECUTION 



>5C5> 

5. REAL I*N*M 

10 • 7 DISPLAY "♦**** 

15. ACCEPT "PRINCIPALa $"#P#" INTEREST* 'M # "NO. MONTHS* 
20. 1=1/12 
25. Q«I+t 

30. M«P*I*Q»N/(<QtN>-i> 

35. DISPLAY "MONTHLY PAYMENT « $"#M 

40. GO TO 7 

45. END 

50. 

>RUN 

PRINCIPAL* $1200 
INTEREST* .06 
NO. MONTHS* 12 

MONTHLY PAYMENT « $ 103.2797157 
*** 

PRINCIPAL* $1200 
INTEREST* .06 
NO. MONTHS* 6 

MONTHLY PAYMENT * $ 203.514547 

PRINCIPAL* $1200 
INTEREST* .06 
NO. MONTHS* 18 

MONTHLY PAYMENT « $ 69.87807671 

PRINCIPAL* $ 
INTERRUPTED BEFORE I 

15. ACCEPT "PRINCIPAL* $"#P#" INTEREST* "* I #"N0. MONTHS* 

15. >SAVE /PR/ 
NEW FILE 



15. >CLEAR 
ERASE PROGRAM? Y 
OK. 
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DOUBLE DECLINING BALANCE DEPRECIATION 



1. DEFINE THE PROBLEM 

The problem is to compute the double declining balance depreciation on any given asset over 
any specified niimber of years using formatted I/O. 

A. Input : 

1. Cost of the asset (C) . 

2. Estimated useful lifetime (U) . 

B. Compute : 

1. Depreciation D = 2 — 

2. Book values C = C - D 

C. Output : 

For the entire range of years . 

1. Year U) 

2. Amount of depreciation CD) 

3. Book value CC) 
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1. FLOWCHART 



^ START ^ 

t 



INPUT 
C, U 



TYPE THE 
HEADING 







) — 


DO 




COMPUTE 




I = 1, U 





D = 2 * C/U 
C = C - D 







^ STOP ^ 



OUTPUT 
I, D, C 



CONTINUE 
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3. 



FORTRAN CODE AND SAMPLE EXECUTION 



-FORTRAN 



>ENTER IOC 10) 



10. C: DOUBLE DECLINING BALANCE DEPRECIATION PROGRAM 
20. WRITE (1*2) 

30. 2 FORMAT (//"COST OF ASSET= $'SZ) 
40. READ (0*3) C 
50. 3 FORMAT CF12.2) 
60. WRITE (1*4) 

70. 4 FORMAT ('ESTIMATED USEFUL LIFETIME^ '^Z) 
80. READ (0*3) U 
90. WRITE (1*5) 

100. 5 F0RMAT(/4HYEAR*8X* 12HDEPRECIATI0N*8X, lOHBOOK VALUE) 
110. DO 100 I=1*U 
120. D=2*C/U 
130. C=C-D 

140. WRITE (1*7) I*D*C 

150. 7 FORMAT ( I 4* 8X* "$"* F 1 0. 2* 8X* "$"* F8. 2 ) 
160. 100 CONTINUE 
170. END 



>SAVE /DDB/ 
NEW FILE 



>RUN 



COST OF ASSET= $3500.00 
ESTIMATED USEFUL LIFETIME= 7. 



YEAR DEPRECIATION BOOK VALUE 

1 $ 1000.00 $ 2500.00 

2 $ 714.28 $ 1785.71 

3 $ 510.20 S 1275.51 

4 $ 364.43 $ 911.07 

5 $ 260.30 $ 650.77 

6 $ 185.93 $ 464.83 

7 $ 132.81 $ 332.02 



STOPPED ATS 160. 
>QUIT 
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LEAST-SQUARE LINE 



1. DEFINE THE PROBLEM 

Fit a least-square line o£ the form Y = A + BX to the following set of data X,Y where X 
is the independent and Y the dependent variable. 



X 


1 


3 


4 


6 


8 


9 


11 


14 


Y 


1 


2 


4 


4 


5 


7 


8 


9 



Solution : 

The worH involves computing first the sums shown below. 



X 


Y 


X2 


X Y 


1 


1 


1 


1 


3 


2 


9 


6 


4 


4 


16 


16 


6 


4 ■ 


36 


24 




S 


64 


40 


9 


7 


81 


63 


11 


8 


121 


88 


14 


9 


196 


126 


EX = 56 


EY = 40 


EX2 = 524 


EXY = 364 



Then compute the regression coefficients A and B 

A - C^Y) (gx2) - (EX) (EXY) 

2 2 
N EX - (EX) 

o N EXY - C2X) (EY) 
B = ~ 

N EX - (EX) 
Then Y = A + BX 



118 



A. Input : 

1. Number of data points (N) 

2. The X values (XCl) - XC8)) 

3. The Y values CYU) - YC8)) 

B. Compute : 

Equations given above: 

1. Sums needed 

2. Coefficients A and B 

C. Output : 

The equation of the least-square line 
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2. FORTRAN CODE AND SAMPLE EXECUTION 



PLEASE LOG IN: 
ACCOUNT: A3 
password: 
USER NAME: MM 
PROJ code: 

READY 4/30 18:30 
-FORTRAN 



FORTRAN IV A21.00 



>10(5) 

10. *C0MMENT :THIS PROGRAM FITS A LEAST SQUARE 

LINE OF FORM: Y=A+BX TO A SET OF DATA (X*Y) 

WHERE IS THE INDEPENDENT VARIABLE. 
15. * 

20. DIMENSION X(8)*Y(8)*XX(8)#XY<8) 

25. accept "NUMBER OF DATA POINTS= 'SN 

30. accept "THE X VALUES ARE "*X 

35. ACCEPT "THE Y VALUES ARE "*Y 

40. *C0MMENTt LOOP TO CALCULATE XX AND XY 

45. * 

50. DO 100 I=1*N 
55. XX<I)«X(I>t2 
60. XY(I>=X(I)*Y<I) 
65. 100 CONTINUE 
70. TX=0. 
75. TY=0. 
80. TXXsO. 
85. TXY=0. 

90. *C0MMENT: loop to CALCULATE TOTALS 
95. * 

100. DO 200 I«1*N 
105. TX=TX+X<I> 
110. TY=TY+YCI> 
115. TXX=TXX-i-XX(I) 
120. TXY=TXY+Xy<I) 
125. 200 CONTINUE 

130. *COMMENT: TO^C^LCULATE COEFFICIENTS A,B 

135. A= C TY*TXX- TX*TXY ) / < N* TXX- TX*TX ) 

140. B=CN*TXY-TX*TY)/CN*TXX-TX*TX) 

145. WRITECU3) A*B 

150. 3 FOHMATC"THE LEAST SQUARE LINE IS: Y="* F6 . 3> "+"> F6 . 3* "X" ) 

160. END 
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I 



>RUN 

NUMBER OF DATA POINTS= 8 

THE X VALUES ARE 1 > 3* 4* 6* 8* 9, U * 14 

THE Y VALUES ARE 1* 2>4#4* 5* 7* 8, 9 

THE LEAST SQUARE LINE ISi Y« 0.545+ 0.636X 

STOPPED ATt 150. 
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PAYROLL CHECKS (FILE I/O) 



1. DEFINE THE PROBLEM 

A. Input : 

1. Read from file /PRF/ employee number (EMPCI)) and pay rate (RATECI)). 

2. Number of hours worked for each employee. 

B. Compute : 

1. Gross pay for each employee. 

C. Output : 

Gross pay file. 

File Formats 

/PRF/ - Input file 

number of employees, Emp^, Rate^^, Emp2, Rate2, .... 

5,99,2.00,77,2.50,88,3.75,55,5.50,66,3.35 

/GPF/ - Output file 
Emp^, Pay^ 

Emp , Pay^ 
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2. FORTRAN CODE AND SAMPLE EXECUTION 



•FORTRAN 



>LOAD /PAY/ 
OK. 

>LIST 



10. DIMENSION EMPC200)#RATE<200)#HRSC200)jPAY<200> 

12. CI OPEN FILE AND READ DATA 

20. OPEN C3#/PRF/#INPUT*SYMBOLIC3 

30. READ <3*7> N 

35. 7 FORMAT <14*Z) 

40. READ <3#8) <EMP< I )*RATE< I )# I^UN) 

45. 8 FORMAT <I3*F7.3#Z) 

50. CLOSE (3> 

54. * OPEN FILE* INPUT NUM.HRS.* COMPUTE PAY AND OUTPUT TO FILE 

55. WRITE < 1*4) 

56. 4 FORMAT <$ EMP.$* 6X*$HRS. $ ) 

60. OPEN C 4* /GPF/* OUTPUT* SYMBOLIC 3 

70. DO 800 Ial*M 

80. WRITE <l*9) EMP<I) 

90. 9 FORMAT <I4*8X*Z) 

100. READ CO* 6) HRSCI> 

110. 6 F0RMAT<F6.2> 

laO. JaEMPCI) 

130. PAY(J)aHRS<I>*RATE<I) 

140. 200 WRITE (4*5) EMPC 1 )*PAY( J) 

145. 5 FORMAT <I3*2X*F7.2) 

150. CLOSE <4) 

160. END 

>RUN 

EMP • HRS . 

99 40. 

77 44.4 

88 48. 

55 39. 

66 45. 



STOPPED ATt 150. 

>f DISPLAY PAYC66) 
150.75 

>QUIT 
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-COPY /GPF/ TO TEL 



99 80.00 

77 Ul.OO 

88 180.00 

55 214.50 

66 150.75 
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THE DIRECTION OF AN AIRCRAFT 



1. DEFINE THE PROBLEM 

A plane is flying from point A to point B. The direction is controlled by a radio direction 
signal from point B. The direction of the plane is corrected for the wind by the autopilot so that 
the plane is always pointed at B. The problem is to write a program to compute the true distance 
flown. 

Let the heading of the aircraft be constant over some small amount of time. The plane will 
then follow path D2, compute the new distance D3 to point B. Repeat this until point B is reached. 
All calculations are done in the complex plane for simplicity. 
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Input : 




1. 


Velocity of aircraft 


CVM) 


2. 


Heading of aircraft 


CM) 


3. 


Velocity of wind 


CWM) 


4. 


Heading of wind 


(AW) 


5. 


Distance 


(DM) 



B. Compute : 

1. Deviation from flight path QDRM) 

2. Total distance traveled (DT) 

C. Output : 

1. Deviation from flight path CDRM) 

2. Total distance traveled CDT) 
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2. 



FORTRAN CODE AND SAMPLE EXECUTION 



10 • 




COMPLEX ud»uo*v *v *D 


12 • 




D I SPiaAY ENTER TnE VELtUuIlX ANU rtl^HiJ i Nva Ur inc. HxnofVfrir i 


14* 




ACCEPT vM*AA 


20 • 




DISPLAY ENTER THE VELOt/Iii ANU DlfilLUllUN Ur inc. wlnu 


22 • 




ACCEPT wM*Aw 


2o • 




DISPLAY THE DIbfANuE Ur inc, lARoEl 


30 • 




ACCEPT L DM J 


32 * 




RADIAN-30U/ . 1 «11 oy ^ 


34 • 




Aw=AW+loU 


36 • 




DRM^O 


38 • 




DI=0 


40 • 




Aw=AW /EAD I AN 


42 • 




A A A A A t% t A IWt 

AA^AA/RADIAN 


44 • 




AA^AtT AA 

AA«Aw-AA 


46 • 




W I =^WK[*S INC AA J 


50 • 




WR-WM4cC OS L AA J 


52 • 




w«CMPLXCWR*Wl J 


54 • 




VasCMPLXC VM# J 


56 • 




Dl «0 


58 • 




D=CMPLXtDM*0 J 


62 • 




DT=0 


64 • 




T= . 1 


66 • 




TT«0 


70 • 




K=2 


72 • 




FLssO 


74 • 


3 


D2»CV+w >*T 


76 • 




D3=D-D2 


84 * 


5 


Dl-CABSCDS J 


86 » 




01=ABSCD1 J 


88* 




DRsEEALCDS J 


90 • 




IFCDI .LT .DRM) GO TO 50 


92* 


40 


DRM-D I 


96 • 


50 


TT«TT+T N 


rt fir 
98 • 




DIsAIMA6lD3 3 


100 • 




DMaCABSCD33 


102. 




D2MaCABSCD2 3 


104. 




DT»DT+D2M 


106* 




VMaCABSCW3 


108. 




VIi«VM*DR/DM 


110. 




VI=VM*DI/DM 


112* 




V=CMPLXCVR,vn 


114. 




DaD3 


il6. 




A«VM*T 


120. 




iF<Dl.GT.A) GO TO 3 


124. 


12 


T«T/10 


126. 




FL«FL+1 


128. 




IFCFL.NE.K) 60 TO 3 


132. 


10 


DISPLAY- LARGEST DEVIATION FROM THE FLIGHT PATH=">DHM 


136. 


585 


DiSPLAY-THE TOTAL DISTANCE TRAVELED WAS 'SDT 


138. 




END 
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>RUN 

EirrER THE VELOCITY AND HEADING OF THE AIRCRAFT 
600*45 

ENTER THE VELOCITY AND DIRECTION OF THE WIND 
200*135 

THE DISTANCE OF THE TARGET 
2000 

LARGEST DEVIATION FROM THE FLIGHT PATH= 256.3855477 
THE TOTAL DISTANCE TRAVELED WAS 2116.221673 



STOPPED ATt 136. 
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CHECKING ACCOUNT PROGRAM 



In this problem, we wish to compute the monthly service charge for a regular checking account, 

The amount of the service charge is based on the average monthly balance and the number of checks 

written. The charge may be computed from the following table. 

AVERAGE MONTHLY BALANCE 



NUMBER 




$200 


$300 


$400 


$500 


$600 


$700 


$800 


$900 


$1000 


$1100 


$1200 


$1300 $1400 $1500 


OF 


UNDER 


to 


to 


to 


to 


to 


to 


to 


to 


to 


to 


to 


to to to 




$200 




$399 


$499 


$599 


$699 




$899 




$1099 


$1199 


$1299 


$1399 $1499 $1599 





J .75 


% .47 


% .33 






$ 




$ 




% 


$ 


$ 


% % % 


1 


.82 


.54 


.40 






















9 
L 


.09 




Ml 


.00 
























.Do 


.Ot 


AO 




















•f 


1 




.UX 






















R 

9 


1 in 




.DO 


.Of 




















6 


1.17 


.89 


.75 


.61 




















7 


1.24 


.96 


.82 


.68 


.54 


















8 


1.31 


1.03 


.89 


.75 


.61 


















9 


1.38 


1.10 


.96 


.82 


.68 


.54 
















10 


1.45 


1.17 


1.03 


.89 


.75 


.61 
















11 


1.52 


1.24 


1.10 


.96 


.82 


.68 


.54 














12 


1.59 


1.31 


1.17 


1.03 


.89 


.75 


.61 














13 


1.66 


1.38 


1.24 


1.10 


.96 


.82 


.68 


.54 












14 


1.73 


1.45 


1.31 


1.17 


1.03 


.89 


.75 


.61 












15 


1.80 


1.52 


1.38 


1.24 


1.10 


.96 


.82 


.68 


.54 










16 


1.87 


1.59 


1.45 


1.31 


1.17 


1.03 


.89 


.75 


.61 










17 


1.94 


1.66 


1.52 


1.38 


1.24 


1.10 


.96 


.82 


.68 


.54 








18 


2.01 


1.73 


1.59 


1.45 


1.31 


1.17 


1.03 


.89 


.75 


.61 








19 


2.08 


1.80 


1.66 


1.52 


1.38 


1.24 


1.10 


.96 


.82 


.68 


.54 






20 


2.15 


1.87 


1.73 


1.59 


1.45 


1.31 


1.17 


1.03 


.89 


.75 


.61 






21 


2.22 


1.94 


1.80 


1.66 


1.52 


1.38 


1.24 


1.10 


.96 


.82 


.68 


.54 




22 


2.29 


2.01 


1.87 


1.73 


1.59 


1.45 


1.31 


1.17 


1.03 


.89 


.75 


.61 




23 


2.36 


2.08 


1.94 


1.80 


1.66 


1.52 


1.38 


1.24 


1.10 


.96 


.82 


.68 


.54 


24 


2.43 


2.15 


2.01 


1.87 


1.73 


1.59 


1.45 


1.31 


1.17 


1.03 


.89 


.75 


.61 


25 


2.50 


2.22 


2.08 


1.94 


1.80 


1.66 


1.52 


1.38 


1.24 


1.10 


.96 


.82 


.68 .54 .00 



A. Input : 

1. Number of checks written (NUMCHKS) 

2. Average monthly balance (AVGBAL) 

3. Current monthly balance (CURBAL) 

B. Compute ; 

1. This month's service charge (MATRIX (I, J)) 

2. New balance (CURBAL) 

C. ' Output : 

1. This month's service charge 

2. New balance. 
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2. 



FORTRAN CODE AND SAMPLE EXECUTION 



LIST 



10. C: CHECKING ACCOUNT 

20. C:i 

30. DISPLAY 'THIS IS A PROGRAM TO COMPUTE THE MONTHLY SERVI 

CE CHARGE FOR A REGULAR CHECKING ACCOUNT AT A COMMERCIAL B 
ANK.* i 



40. 


CSI 












50. 




REAL MATRIX* NUMCHKSI 










60. 




DIMENSION MATRIX (15# 26 >i 










70. 




WRITE< 1*50)1 










80. 


50 


FORMAT CURRENT BALANCE = * *Z)1 










90. 


75 


FORMAT </ AVERAGE BALANCE FOR THIS MONTH 


s 


• *Z)J 




100. 


80 


FORMAT</* NUMBER OF CHECKS THIS MONTH = * 


*Z)i 






1 10* 




ACCEPT CURBALl 










1 20* 




DO 100« Isl*15« 










i 30* 




IFCI.LT. 15) MATRIXCI,26) « 2*50 - <I*.14)i 










1 40. 




IF<I.E0. 15) MATRIXCI#26) = O.OOl 










I 50. 




IFCI.EQ. 1) MATRIX(I#26) « 2.50j 










160. 


100 


CONTINUE* 










I 70. 




DO 200* I«l*15i 










180. 




DO 200* Jai*25l 










190. 




MATRIX<I*J) s MATRIX<I*26)-(C26-J)*.07)J 










200. 


200 


CONTINUEI 










210. 




WRITE<1*75)I 










220. 




ACCEPT AVGBALl 










230. 




WRlTECl*80)i 










240. 




ACCEPT NUMCHKSI 










250. 




I s AVGBAL/lOOi 










260. 




J s NUMCHKS -t- li 










270. 




IFCMATRIX<I*J) .GE. 0.54) GO TO 4001 










280. 




MATRIXCI*J) s 0.001 










290. 




IF CCI .E0.2.0.AND. J.EQ.l .0) .OR. CI.EQ.3.0.AN0. 


J. 


EQ. 


3. 


0) 






.0R.<I.EQ.4.0.AND.J.e0.5.0)> MATRIX<I*J) s 





.47 


} 




300. 




IF ((I.Ed.3.0.AND.J.£0.2.0).OR.(I.EQ.4.0*AND. 
MATRIX <I*J) a 0.40i 


J. 


EQ. 


4. 


0) ) 


310. 




IF C(I.E0.3.0.AND. J.EO. 1 .0 ) • OR* ( I .EQ .4.0 • AND. 
MATRIX <I*J) s 0.331 


J. 


EQ. 


3. 


0)) 


320. 


400 


DISPLAY 'THIS MONTHS SERVICE CHARGE =**MATRIX(I 


* J) 


i 




330. 




CURBAL = CURBAL - MATRIX(I*J)1 










340. 




DISPLAY 'THE NEW CURRENT BALANCE #CURBALl 










3 50. 




ENDJ 











> EXECUTE 

THIS IS A PROGRAM TO COMPUTE THE MONTHLY SERVICE CHARGE FOR A REGULAR 
CHECKING ACCOUNT AT A COMMERCIAL BANK. 
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CURRENT BALANCE = 513.67 



AVERAGE BALANCE FOR THIS MONTH « 336.71 

NUMBER OF CHFXKS THIS MONTH = 13 
THIS MONTHS SERVICE CHARGE = 1.24 
THE NEW CURRENT BALANCE = 512.43 
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SPECTRAL EFFICIENCY OF A BLACK BODY 



The distribution of energy emitted from a black body, with respect to wave length is given by 
Planck's distribution law: 

ECA) dX = 4 Ce ^^^^^ -1)'^ dA (1) 
A^ 

where and C2 are the first and second Planck radiation constants. 

-15 

Ci = STThc = (4.9918 t 0.002) X 10 erg-cm. 
C2= hc/k = 1.43880 i 0.0007 cm-deg. 

The total density of radiant energy in unit volume is (1), integrated over all wave lengths. 



8 TT CkT)^ 
Total - j^3^3 



8 TT^ k"^ 

where a = 



u^du ^ 8 Tr^CkT)^ ^ ^ 
e^-1 15h3c3 



ISh^c^ 

(2) is known as the Stefan-Boltzmann law. 

If we consider radiant energy in terms of flux, the total radiation crossing a unit area in 
unit time in all directions in one hemisphere is usually written: 



K . ^ = aT^ 
total 



271 ^ k^ - 5 erg 
where a (Stefan-Boltzmann Constant) = = 5.672 x 10 ^ — 

ISh^c^ cm^deg'^ 

For practical calculations this reduces to: 

, , .4 E = total radiant flux in watts/cm^ 

total ~ T = Temperature in degrees K. 

If we integrate Planck's equation between the wavelength interval in which we are interested 
and compare this to the total radiation, we can obtain the relative efficiency of a black body emitter, 
for a selected wavelength interval, as a function of temperature. 

E = b f x^ fe 1-4388/Tx where x = wavelength in cm. 

partial J 
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For Epay^iai in watts/cm^ b = 2.39 X lO'^^V 

Relative Efficiency = E = E ^. ,/E , x 100 
' partial' total 

which reduces to 

E^^^ = 64.774 j X5 Ce^-^^^^^/^^-l)-^ dx/T^ 
^1 

The integral can be evaluated by Simpson's 1/3 rule as follows: 

'"^ 

E = jT f(x)dx = I [fCx^) + 4 f Cxi+h)+2 fCxj + 2h) 
Xi 

+ 4 f Cxi + 3h) + — + 2 f - 2h) 
+ 4 f CX2 -h) + f CX2)] 

, , X2 - Xi 

Where h = where n = number of terms in the series; n is even. 

We can evaluate the integral easily by separately summing the terms multiplied by 4 and the 
terms multiplied by 2 and then adding in the first and last terms. 

Since the Planck function is used in several places, it is convenient to use the DEFINE 
statement. 

If the user uses an n = 100, he will get results accurate to the five significant digits 
printed. In fact, n = 10 is quite accurate for temperatures over 5000° K. The reader may want to 
experiment with this. The following chart shows some sample results of such experimentation. 
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SPECTRAL EFFICIENCY PROGRAM 
Test for sensitivity to N over visible spectrum 





Xi = 


4.0 X 


10"^ 




= 7.0 X 


10-^ 








N = 


1000 




100 




10 






2 




TEMP. 




















500° 


1.5850- 


■12 


1.5850- 


■12 


1.6222- 


■12 


4. 


,3237- 


•12 


1000° 


1.8041- 


■04 


1.8041- 


■04 


1.8054- 


■04 


2. 


,3938- 


■04 


5000° 


3.1033 


01 


3.1033 


01 


3.1034 


01 


3. 


,1065 


01 


10000° 


3.2117 


01 


3.2117 


01 


3.2117 


01 


3. 


,2040 


01 


15000° 


1.8245 


01 


1.8245 


01 


1.8245 


01 


1. 


,8315 


01 


20000° 


1.0464 


01 


1.0464 


01 


1.0464 


01 


1. 


,0554 


01 



DEFINE THE PROBLEM 

A. Input ; 

1. Upper and lower wave length boundaries (XI, X2) . 

2. Temperature ranges in degrees Kelvin (Tl, T2) . 

3. The temperature interval CT3) . 

4. The number of integration intervals CN) . 

B. Compute : 

The spectral efficiency of a black body. 

C. Output : 

1. Temperature (T) 

2. Percent efficiency (E) 
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2. 



FORTRAN CODE AND SAMPLE EXECUTION 



>RUN 

mVELENGTHCCM) RANGES FROM Xl« 4.0E-5 TO Xa= 7.0E-5 
TEMPERATURE (DEGREES K) RANGES FROM! 

Tl=500 TO T2a200Q0 IN STEPS OF T3a500 
EVEN NUMBER OF INTEGRATION N«aO 

TEMPERATURE PERCENT EFFICIENCY 



500 


• 1587E-1 1 


1000 


V1804E-03 


i 500 


.5453E-01 


2000 


• 7629E+00 


asoo 


.3274E+01 


3000 


.7941E+01 


3500 


• 1402E-i'02 


4000 


•2043E+0a 


4500 


• 2626E'f02 


5000 


•3103E+02 


5500 


•3455E+02 


6000 


•3688E+02 


6500 


•3815E+02 


7000 


•3857E+02 


7500 


.3832E+02 


8000 


•3758E+02 


8500 


.3649E+02 


9000 


.351 5E+0S 


9500 


.3367E+02 


10000 


•321 lE+02 


10500 


.3052E+02 


1 1000 


.2893E+02 


11500 


•2738E+02 


12000 


.2587E+02 


12500 


.2442E+02 


13000 


.2304E+02 


13500 


.2174E+02 


14000 


•2050E+02 


14500 


. 1934E+02 


1 5000 


.1824E+02 


15500 


.1721E+02 


16000 


.1625E+02 


16500 


•1535E+02 


17000 


.1450E+02 


17500 


. 1371E+02 


18000 


.1297E+02 


18500 


.1228E+02 


19000 


.1163E+02 


19500 


.1103E+02 


20000 


.1046E+02 



PAUSE 
ISO. > 
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•FORTRAN 



>LOAD /PR/ 
OK. 

>LIST 



10. 3 ACCEPT WAVELENGTH C CM) RANGES FROM XI « "#X1*"T0 X2« **,X2 

20. ACCEPT "TEMPERATURE (DEGREES K) RANGES FROM Tt=*STl#"TO T2= 

"*T2#"IN STEPS OF T3»'ST3 

30. ACCEPT "EVEN NUMBER OF INTEGRATION N=*SN 

40. IF <MODCN#23.NE.O) N»N+1 

50. Ha<Xa-Xl)/N 

60. WRITE (1#9) 

70. 9 FORMAT (/8X,$TEMPERATURE$* 10X# SPERCENT EFFICIENCYS) 

80. Fi:A*B3«1.0/<At5*(EXPE1.4388/<B*A)3-l)) 

90. DO 100 T«Tl#Ta*T3 

100. E«<64.774>*(CH/3)*(FCXl#T3+4*<SUMCXl+H*X2-H*2*H*T3)"i-2*<SUM 

CXJ+2*H#X8-2*H*2*H#T3)+FCX2#T3))/Tt4 

110. WRITE (1,10) T,E 

120. 10 FORMAT ( 1 3X, I 5# 1 6X,E1 1 .4) 

130. 100 CONTINUE 

1 40 . PAUSE 

150. 60 TO 3 

160. END 



170. FUNCTION SUMCW,X,Y,T3 

180. S«0 

190. FCA#B3«1 .0/(At5*(EXPC1.4388/(B*A)3-l )) 

200. DO 200 R»W,X,Y 

210. S«FCR#T3+S 

220. 200 CONTINUE 

225. SUM«S 

230. RETURN 

240. END 
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THE BASE OF A SEMICIRCLE 

1. DEFINE THE PROBLEM 

Given n number of chords inscribed within a semicircle and their respective lengths, find the 
diameter of the semicircle. 

If we construct triangles with the radii of the semicircle and bisect each chord length with 
a construction line, C, we have n half angles, ^ whose sum is equal to 90°, or -n/l radians; from 
the sum of the angles equal to 180°. First, we approximate the angle ^: 

(j) = [(sideCl)/2) /sum of sides] 

We also know that the angle 6 is equal to 90°, from this we may approximate the diameter using 
the law of sines; 

Cdiameter/2)/sin <)> = CsideCl)/2) /sin ^ 

Knowing that the sine of 90° is equal to 1.0000 we can see: 

diameter = 2 [(side CI) /2) /sin c|)] 

Since our first calculation of (j) was only a rough approximation and thus the diameter was 
only a rough approximation, we now want to approximate the diameter to a greater degree of accuracy. 

We know that when the sum of the half angles is equal to 7r/2, or 1.57079, we have a very close 
approximation of the diameter of the semicircle. We may approximate the half angles by using the 
tangent and arctangent functions . 

We know that; 

sideCl)/2 

C (construction line) 



since. 



Therefore : 



and that; 



C = Cdiameter/2)2 - Cside(I))2 
Cside(I)/2) 



tan = 



= tan 



[ (diameter/2) 2- (side (1) ) ^] ^ 

-1 (side CI) /2) 

[(diajneter/2)2-CsideCI)/2)2]l/2 



From the above calculation, if the sum of the angles is greater than 7r/2, or 1.57079, the diam- 
eter must be made larger by some increment. If the sum of the angles is equal to it/2, the approximation 
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is exact to five decimal places, and if the sum of the angles is less than Tr/2, the approximation 
the diameter is too large and must be reduced by some increment. 

Following is a program to solve the problem by incremental approximations of the diameter. 




I ~— DIA 

A. Input : 

1. Niomber of sides CN) 

2. The sides SIDE(N) 

B. Compute ; 

Do computations explained above. 

C. Display : 

1. Length of the diameter. 
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2. FORTRAN CODE AND SAMPLE EXECUTION 

10. CS BASE OF A SEMI-CIRCLE 

20. c: given: the number of sides and their lengths 

40. ACCEPT " THE NUMBER OF SIDES = "»N 

70* CALL CALCCN] 

80* END 

90* SUBROUTINE CALCCN3 

100* DIMENSION SIDE(N> 

110* REAL INCR 

120* SUNSIO « 

130* CNTR « 

140. ALPHA « 

150* MAX « 

160* MIN • 1 

170. DO 30* Ia|*N 

180. ACCEPT " SIDE « "# SI DEC I) 

210* SUMSID • SUMSI0'*^SI0E<I> 

220. MAX « AMAX1CMAX«SIDE<I>3 

230* 30 CONTINUE 

240. ALPHA « CCMAX/2>/SUMSlO>*3«14l59 

250* DIAl > (NAX/2>/SlNCALPHA] 

260* DIA « OiA|4c2 

270. INCR • MAX/ 10 

280* 100 AW6LE « 

290* DO 40* I»1*N 

300* y « SIDECI) 

310. X « S0RTC:<DIA*«2)-CSIDE<I)**2)3 

320. 40 AN6LE » ANGLE * ATANCY/Xl 

330. IAN6LE « ANGLE^IOOOOO 

340. IF<I57079-IAN6LE)50*60*70 

350. SO DIA « DIA ♦ INCR 

360. 60 TO lOO 

370. 70 DIA ■ DIA - INCR 

380. INCR • INCR/tO 

390. OIA - DIA * INCR 

400. CNTR « CNTR * 1 

410. IF (CNTR-20>100*100*60 

420. 60 OISPLAY • THE BASE OF THE SEMI-CIRCLE «• *OIA 

430. RETURN 

440. END 
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>EXECUTE 

THE NUMBER OF SIDES » 3 
Slot « t 
SIDE s 2 
SIDE « 3 

THE BASE OF THE SEMI-CIRCLE ■ 4- 11 3103501 



>EAECUTE 

THE NUMBER OF SIDES « 2 
SIDE « 1 
SIDE a I 

THE BASE OF THE SEMI-CIRCLE • I •4142145 
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APPENDIX A 



EFFICIENT STORAGE ALLOCATION 

The computer has a specified amount of memory available for both program and data. By seg- 
menting your program; that is, storing only part of it in memory and the rest on the drum (another 
data storage device) , the total number of statements and amount of data in a program may be increased 
greatly. 

A maximum of four pages of memory is available for both the program (segment) and the data. 
The amount of storage needed for data must be estimated. The amount of storage required for each 
program segment may be found using the CCS command MAP. 

A page may not be shared by the program statements and the data. Therefore if part of a page 
is used for data storage, the entire page must be considered filled. 

The maximum size of the segment may be estimated as follows: Maximum Segment Size = 4 pages 
minus the data storage requirements. Data storage needs may be estimated using the following values: 

1 Page = Approximately 2000 Words 

1 Real Number = 2 Words 

1 Integer Number = 2 Words 

1 Complex Number = 4 Words 

1 Double Precision Number = 3 Words 
Thus, if a program contains the real array AC130,20), 5200 words or 3 pages would be required for 
data storage leaving one page for program segment storage. 

The actual segment size may be determined using the CCS command MAP. The MAP command prints 
a table as shown in the example below. 

>MAP 

TEXT PAGES =2.8 
UNUSED PAGES = 3 

SEGMENT 1 2 3 

STATEMENTS 150 175 127 130 

PAGES USED 0.6 0.8 0.5 0.9 

> 

If any of the segments exceeds the maximum size allowed, they must be reduced in size before 
the program will run. 
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APPENDIX B 



INTERNAL REPRESENTATION OF ASCII CODE 



Octal 



Octal 



!n tat ion 


Character 


Representation 


Character 


00 




AO 


@ 


01 


1 


41 


A 


02 


1 1 


A? 


B 


03 


ff 


43 


c 


04 


$ 


AA 


D 


05 


9- 

'D 


45 


E 


06 


e 
Q 




F 


07 


t 


47 


G 


10 


c 

I 

^ 

J 


50 


H 


11 


51 


I 


12 


* 


52 


J 


13 


+ 


53 


K 


1 A 




54 


L 


15 




55 


M 


16 




56 


N 


17 


/ 


57 





20 





60 


P 


21 


1 


61 


Q 


22 


2 


62 


R 


23 


3 


63 


S 


24 


4 


64 


T 


25 


5 


65 


U 


26 


6 


66 


V 


27 


7 


67 


w 


30 


8 


70 


X 


31 


9 


71 


Y 


32 




72 


z 


33 


> 


73 


[ 


34 


< 


74 


\ 


35 




75 


] 


36 


> 


76 


-^ 


37 


? 


77 
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SUMMARY - FORTRAN IV LANGUAGE STATEMENTS 

ABBREVIATIONS 

In. = line number 
var = variable 
exp. = expression 
vl, = variable list 
sn. = statement number 
val = value 

arg. list = argument list 

< > = anything enclosed in < > is optional. 

CONSTANTS . 

Integer e.g. 1, 400 

Real e.g. 1.9, 1,7E 07, 5E-55 

Complex (3,5.2) 

Logical .TRUE., .FALSE. 

Hollerith 4HTIME, "XXX", 'RED', $RTM$ 

ARITHMETIC OPERATORS (In Order Of Priority) 
t or ** Exponentiation 

Unary minus (negative) 
*, / Multiplication, Division 
+, - Addition, Subtraction 

RELATIONAL OPERATORS 

.EQ. Equal to 

.NE. Not equal to 

.LT. Less than 

.LE. Less than or equal to 

.GT. Greater than 

.GE. Greater than or equal to 

LOGICAL OPERATORS (In Order of Prionty) 
.NOT. 
.AND. 
.OR. 
.EOR. 
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ASSIGNMENT STATEMENTS 

variable = expression 

INPUT/OUTPUT STATEMENTS 



ACCEPT 



vl, "text", $text$, 'text' 



DISPLAY vl, "text", $text$, 'text', expression 

Formatted I/O (symbolic) 




file format 
no. no. 



list 



File Niombers 

= Terminal Input 

1 = Terminal Output 
2,3,4 = disk file I/O 



Open & Close Disk File 

OPEN 



file /file, 
no. 'name' ' 



INPUT SYMBOLIC 

or , or 
OUTPUT BINARY 



CLOSE 
Binary I/O 

READ 

or 
WRITE 



file 
no. 



r fii( 

\no. 



list 



Commands From Files 

COMMANDS /file name/ 

CONTROL STATEMENTS 

IF (logical exp) executable statement 
IF (arithmetic exp) sn^^, sn2, sn^ 

i'^i't final ^n-f-Prva1> 

DO sn var = , interval> 



GO TO<sn>or<var>or< 



segment^ 
name 

GO TO (sn^^, sn2, sn^...) var 



CALL< 



subroutine 



name 

ASSIGN integer TO var 
CONTINUE RETURN 
PAUSE FINISH 
STOP END 



r actual 1 ^segment ^ 
arg. list name 
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DECLARATION S TATEMENTS 
comments 



•[ 

DATA I var^/val^/,var /val /... 

I^var^, var2,...7val^, vdl^,.,./ 

DIMENSION ^^^^y / • "PP®^ 1°^®^ . "PP®^ 

name y limit/ limit^, limit2* '^^^^'^2*' 



INTEGER 
REAL 

DOUBLE PRECISION 

COMPLEX 

LOGICAL 



vl or label 



function name 

FUNCTION faction 
name 



SUBROUTINE subroutine ^ 



[[ dummy 1 
arg. list 
dummy 
arg. listj 

[dummy 
arg. listJ' 



SEGMENT segment 
name 

COMMON vl. 

Cfield spec) 
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SUMMARY - CCS FORTRAN IV COMMANDS 



Line number = .001 - 999.999 



<ENTER> 



In FORTRAN statement 

In: In (D^ terminator) 

In (interval) <ln> CD^ terminator) 



COPY 



/file name/ 
In: In 
TELETYPE 



TO 



/file name/ 
In: In 
TELETYPE 



MOVE In: In TO 



/file name/ 

In: In 
TELETYPE 



RENUMBER 



LIST 
FAST 

DELETE 



(starts at 100 in steps of 10) 

In: In 

In: In AS In: In 



In 

In: In 



(entire program) 



REFERENCES 
DEFINITIONS 



(current line) 

["label <ln:ln> 
Lor variable 



EDIT 
MODIFY 



In 

In: In 



BREAK In In: In 
RESET 

CHECK 

INITIALIZE 

SHRINK 

LOCK 

MAP 

@ FORTRAN statement 



(maximum of 10 lines or line ranges) 



RUN 

EXECUTE 

NEXT 

CLEAR 

QUIT 

CONTINUE 



SAVE /file name/ 
LOAD /file name/ 
COMMANDS /file name/ 
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